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CHAPTER  1 


INTRODUCTION 

1.1  Statement  of  the  Problem 

Equivalent  circuits  containing  frequency  dependent  elements,  or 
circuits  valid  only  in  narrow  bands  of  frequencies,  have  been  used  for 
years  by  the  electromagnetics  community  to  model  the  terminal  behavior 
and  receiving  properties  of  radiating  and  scattering  structures.  In 
contrast,  this  work  is  concerned  with  the  synthesis  of  broad-band 
equivalent  circuits  with  frequency  and  time  invariant  elements,  for 
electromagnetic  energy  collecting  structures  with  definable  ports. 

Our  goal  can  be  perhaps  best  explained  with  the  help  of  Figure 
1.1.  In  Figure  1.1(a)  is  shown  an  energy  collecting  structure  (a 
ylindrical  dipole  in  this  case)  excited  by  a  plane  wave  having  some 
time  history  f(t),  and  loaded  at  its  terminals  by  a  general  impedance 
Z^,  which  is  permitted  to  be  nonlinear  and/or  time-varying.  Our  task 
is  to  proceed  from  the  known  mathematical  description  of  the  structure 
and  a  given  specification  of  the  spatial  character  of  the  impinging 
field,  to  synthesize  an  equivalent  active  network  (Figure  1.1(b)) 
consisting  of  a  finite  number  of  lumped,  linear  circuit  elements,  whose 
behavior  at  the  terminals  is  a  good  approximation  to  that  of  the  original 
object  for  any  loading  conditions.  The  network  should  be  driven  by  an 
autonomic  voltage  source  with  a  waveform  f(t)  and,  with  the  source  shorted, 
it  should  replicate  the  driving-point  immittance  of  the  structure  with 
no  external  incident  field  present.  This  equivalent  network  should 
also  be  applicable  in  the  time-harmonic  mode  of  operation. 

The  question  which  naturally  arises  is,  "Why  would  one  want  to 
construct  an  equivalent  circuit  for  an  electromagnetic/radiating 
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(b) 


Figure  1.1  Illustration  of  the  synthesis  problem,  (a)  Energy-collecting 
structure  excited  by  a  plane  wave  with  a  time  history  f(t)  and  headed 
at  a  prot  by  a  load  Z^;  (b)  its  equivalent  circuit  representation. 
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problem?"  Beam  [1976a]  in  his  work  introducing  the  singularity  expansion 
based  equivalent  circuit  synthesis  states  that  these  circuits,  under 
certain  circumstances,  could  be  helpful  in  providing  some  or  all  of  the 
following: 

1.  physical  insight 

2.  computational  convenience 

3.  capability  of  using  established  circuit  transformation 
techniques 

4.  capability  to  combine  the  electromagnetic  analysis  with  physical 
circuit  elements,  transmission  lines,  etc. ,  which  are  constructed 
as  part  of  an  antenna  or  scatterer 

5.  capability  of  using  existing  computerized  circuit  analysis 
programs. 

In  a  more  recent  paper  Baum  and  Singaraju  [1980]  add  to  this  list  two 
new  areas  of  application: 

6.  physical  construction  of  equivalent  circuits  for  use  in  special 
types  of  electromagnetic  pulse  (EMP)  simulators 

7.  radar  target  detection  and  camouflage  techniques. 

The  desired  features  of  the  equivalent  networks  depend,  of  course, 
on  the  area  of  application.  In  general,  they  should  be  simple,  easy  to 
construct,  and  accurate  in  the  frequency  band  of  interest.  It  is 
desirable  that  the  networks  be  physically  realizable,  i.g.,  they  should 
not  contain  any  negative  or  complex-valued  elements.  It  is  self-evident 
that  this  requirement  is  essential  if  the  equivalent  circuit  is  to  be 
physically  constructed.  This  feature  is  not  essential,  but  it  is 
potentially  helpful  in  other  applications,  because  negative-valued 
elements  are  difficult  to  interpret  physically  and  they  can  cause  numerical 
problems  in  some  circuit  analysis  codes.  A  potential  application  of  the 
hardware-realized  circuits  is  in  the  high-voltage  environment  of  transient 
pulse  simulation.  This  environment  precludes  the  use  of  any  controlled 
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sources,  gyrators,  and  ideal  or  perfect  transformers,  because  these 
elements  can  only  be  realized  by  active  devices  which  usually  do  not 
tolerate  high  voltages  and  currents.  It  is  also  desirable  that  the  angle 
of  incidence  and  other  parameters  of  the  incident  field  enter  the 
circuit  in  the  simplest  possible  way.  For  example,  we  would  prefer 
that  a  change  in  the  angle  of  incidence  8  in  Figure  1.1(a)  would  change 
perhaps  some  element  or  source  values  but  not  the  topology  of  the  circuit. 
Furthermore,  neither  the  topology  nor  the  element  values  should  depend 
on  the  time  history  f(t).  It  would  be  also  convenient  for  the  autonomic 
source  to  be  connected  to  only  one  pair  of  terminals  in  implementation. 

Perhaps  the  most  convenient  and  powerful  mathematical  description 
of  an  electromagnetic  energy  collecting  structure  for  purposes  of 
equivalent  circuit  synthesis  is  the  Singularity  Expansion  Method  (SEM) 
introduced  by  Baum  [1971]  (see  also  [Baum,  1976b]  and  [Baum,  1976c]. 

SEM  is  a  generalization  of  well-known  techniques  of  linear  circuit  theory 
in  which  the  singularities  of  a  transfer  function  are  used  to  determine 
the  transient  response  by  the  Heaviside  expansion  theorem.  In  electro¬ 
magnetic  theory,  the  singularities  are  found  by  first  applying  a  two- 
sided  Laplace  transform  to  Maxwell's  equations  and  then  constructing 
an  integral  equation  for  the  scattered  field.  Complex  singularities 
{s^}  appear  as  poles  of  the  inverse  of  this  equation  and  these  are 
the  frequencies  for  which  non-trivial  solutions  of  the  corresponding 
homogeneous  integral  equation  exist.  These  non-trivial  solutions  are 
termed  the  natural  modes  of  the  structure.  From  these  modes  and  the 
integral  operator  the  so-called  normalization  coefficients  can  be  found. 
Poles,  natural  modes  and  normalization  coefficients  constitute  the 
fundamental  SEM  parameters.  If  the  poles  and  the  associated  modal 
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distributions  are  known,  one  can  compute  the  coupling  coefficients  which 
weight  a  given  singularity's  contribution  to  the  response  for  a  given 
incident  field.  From  these  quantities  and  the  normalization  coefficients 
the  residues  associated  with  the  poles  can  be  evaluated.  Since  a 
distributed-parameter  system  has  an  infinite  order  of  complexity,  it 
can  only  be  described  by  an  infinite  series  of  terms  associated  with 
poles  and  other  singularities  (see  Section  2.4).  In  practice,  however, 
for  band-limited  excitations,  the  response  of  such  system  is  dominated 
by  a  limited  number  of  resonances,  so  that  the  series  can  be  truncated 
without  introducing  appreciable  error. 

Following  this  introduction,  the  main  objective  of  this  work  can 
be  stated  as  follows:  To  develop  simple  and  reasonably  accurate 
procedures  for  the  synthesis  of  finite,  lumped-parameter,  physically 
realizable  circuits  for  electromagnetic  energy  collecting  structures 
from  their  SEM  description.  These  procedures  should  be  practically 
tested  on  some  representive  canonical  structures,  their  responses  for 
broadband  excitations  should  be  computed  and  compared  with  available 
data  or  with  responses  obtained  by  other  means. 


1.2  brief  Literature  Survey 

One  of  the  early  attempts  to  extend  circuit  theory  concepts  to 
distributed-parameter  system  was  made  by  Sohelkunoff  [1944],  who 
applied  the  Foster  [1924]  expansion  theorem  to  immittances  of  lossless 
structures  with  an  infinite  number  of  degrees  of  freedom.  He  observed 
that  this  infinite  series  representation  was  slowly  convergent  or  did 
not  converge  at  all.  He  then  used  a  Mittag-Lef f ler  theorem  from  the 
theory  of  functions  of  a  complex  variable  to  convert  it  into  a  more 
rapidly  convergent  form  (he  in  effect  introduced  what  we  call  in  Chapter 


2  Che  modified  pole  modules).  Subsequently,  Schelkunoff  modified  this 
series  to  allow  for  slight  dissipation.  He  did  not  present  any  equivalent 
circuits  for  this  case. 

In  a  later  work,  Schelkunoff  and  Friis  [1952]  discussed  the 
representation  of  driving-point  immittances  of  antennas  in  terms  of 
their  poles  and  zeros.  They  conjecture  there  that  the  impedance  of 
any  physical  circuit  may  be  expressed  as  a  ratio  of  two  products  of 
linear  factors  exhibiting  the  natural  oscillation  constants  of  the 
circuit  with  its  terminals  at  first  floating  and  then  short-circuited. 

For  a  distributed-parameter  system  the  number  of  factors  is  infinite 
and  the  convergence  of  this  product  form  is  governed  by  the  Weierstrass 
theorem  from  the  theory  of  functions  of  a  complex  variable.  Then  they 
state  that  the  static  capacitance  (inductnace)  of  a  dipole  (loop)  antenna, 
together  with  its  natural  oscillation  constants,  determines  the  antenna 
impedance  at  all  frequencies.  They  do  not  present  any  equivalent 
circuits  for  antennas,  except  for  parallel  and  series  resonant  circuits 
valid  in  the  vicinity  of  the  antenna  resonance. 

An  explicit  expression  for  the  driving-point  admittance  can  be 
written  only  for  a  few  structures,  such  as  the  spherical  or  spheroidal 
antenna  with  circumferential  slot.  This  expression  is  in  a  form  of  an 
infinite  series  of  terms  associated  with  the  eigenfunctions  of  the 
problem.  For  the  spherical  antenna  [Stratton  and  Chu,  1941]  the  nth 
term  involves  a  ratio  of  a  spherical  Hankel  function  of  order  n  and  its 
derivative.  As  shown  by  Chu  [1948],  this  ratio  can  be  expanded  in  a 
finite  continued  fraction  with  positive  real  coefficients,  which  can 
be  physically  realized  by  an  LC  ladder  network  terminated  into  a  single 
resistance.  A  complete  equivalent  circuit  consists  of  a  parallel  connection 
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of  an  infinite  number  of  such  "eigennetworks".  Unfortunately,  this 
representation  is  not  practical  because  it  converges  very  slowly. 
Sohelkunoff  [1952]  and,  more  recently,  Franoesohetti  [1976]  (see 
also  [Buaci  and  Franceschetti,  1974])  observed  that  each  of  the  eigen¬ 
networks  contributes  to  the  static  capacitance  C  of  the  antenna.  They 
were  able  to  eum  explicitly  the  infinite  series  for  C.  The  subtraction 
of  this  series  term-by-term  from  the  series  for  the  admittance  resulted 
in  a  more  rapidly  convergent  series.  Franceschetti  constructed  an 
equivalent  circuit  representation  for  this  modified  series  by  intro¬ 
ducing  gyrators  (with  frequency-dependent  gyration  resistance)  to 
each  eigennetwork.  Neither  Schelkunoff  nor  Franceschetti  discussed 
the  source  problem. 

An  important  development  for  the  equivalent  circuit  synthesis  for 
energy  collecting  structures  was  the  introduction  in  1971  of  the 
Singularity  Expansion  Method  [Baum,  1971].  Barnes  [1973]  constructed 
perhaps  the  first  SEM  driving-point  equivalent  circuit  for  a  thin  dipole 
antenna.  He  used  the  approximate  poles  and  residues  derived  analytically 
by  Lee  and  Leung  [1972].  In  his  circuit  each  pole  pair  was  represented 
by  a  series  RLC  module.  Hess  [1975]  has  used  SEM  analysis  to  develop 
a  transfer  function  model  for  determining  EMP  coupling  to  a  dipole  antenna. 
He  used  the  SEM  data  obtained  numerically  by  Tesche  [1973].  He  did 
not,  however,  synthesize  any  equivalent  circuits  so  no  loading  effects 
could  be  analyzed  by  this  method. 

Perhaps  a  more  systematic  approach  to  SEM  equivalent  circuit 
synthesis  was  propsoed  by  Baum  [1976a]  (some  of  this  material  appeared 
recently  in  [Baum,  1978]).  In  this  work  he  introduced  the  formalism  and 
terminology  which  stimulated  and  influenced  most  of  the  subsequent 
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research  in  this  area,  including  the  present  effort.  In  his  paper, 
however,  Baum  developed  only  formal  equivalent  circuits  without  addressing 
the  issues  of  their  practical  realizability.  In  a  latter  work  Baum  and 
Singarao'u  [1980]  utilized  a  sinusoidal  natural-mode  approximation  and  an 
analysis  similar  to  that  presented  by  Marin  [1974],  to  derive  driving- 
point  admittance  equivalent  circuits  for  a  thin-wire  dipole  antenna. 

They  did  not,  however,  present  any  active  equivalent  circuits  for  an 
illuminated  dipole. 

Sohaubert  [1979]  constructed  lumped-parameter  equivalent  circuits 
for  a  center-driven  dipole  and  a  Yagi  antenna  using  experimentally 
derived  SEM  description.  He  used  Prony's  algorithm  [Van  Blccricum  and 
Mittra,  1975]  (a  technique  for  parametrization  in  terms  of  complex 

0 

exponentials)  to  extract  SBM  data  from  the  antenna  terminal  voltage 
waveform  due  to  a  step-like  excitation  and  used  these  data  to  construct 
the  driving-point  impedance  function  in  the  form  of  a  ratio  of  two 
polynomials  in  s.  He  then  used  Brune's  [Brune,  1931]  synthesis  method 
to  derive  equivalent  circuits.  As  usual  in  Brune's  method,  perfect 
transformers  were  employed,  Schaubert  did  not  consider  the  source 
problem  but  limited  his  concern  to  the  determination  of  the  driving- 
point  impedance  part  of  the  equivalent  circuit. 

The  first  extensive  study  of  the  physical  realizability  of 
passive  equivalent  circuits  for  thin-wire  dipole  and  loop  antennas 
was  performed  by  Streable  and  Pearson  [1981],  They  utilized  the  Bott- 
Duffin  [Bott  and  Duffin,  1949]  synthesis  method  to  realize  input  admittances 
using  different  grouping  of  pole  terms  and  used  the  resistive  padding 
technique  [Guillemin,  1977]  for  an  approximate  synthesis  of  pole-pair 
circuits.  In  a  companion  paper,  Pearson  and  Wilton  [1981]  proved 
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that  the  eigenadmittances  (reciprocal  eigenvalues)  of  a  structure  are 
positive  real  (PR)  [ Brune ,  1931]  functions  of  frequency  s  and  speculated 
on  the  PRness  of  the  admittances  associated  with  individual  current 
eigenmode  contributions  to  total  current.  They  also  briefly  discussed 
the  source  synthesis  issue,  but  they  did  not  present  any  actual  active 
circuits. 

Another  approach  to  equivalent  circuit  synthesis  has  been  to  compute 
the  driving-point  or  transfer  function  of  the  structure  on  the  imaginary 
axis  in  the  s-plane  and  then  to  use  a  curve  fitting  technique  to 
develop  the  corresponding  rational  function  in  s,  which  can  be 
subsequently  submitted  to  some  standard  synthesis  procedure  [Weber  and 
Toulios ,  1974;  Sharpe  and  Rousszt  1979].  This  technique  was  used  by 
Sharpe  and  Roussi  [1979]  to  obtain  a  ladder  network  representation  for 
the  input  admittance  of  a  biconical  antenna  [ Tai ,  1981].  This  method, 
however,  does  not  offer  the  advantages  of  the  SEM  appraoch,  as  discussed 
above . 

Recently  Tesche  and  Giri  [1981]  presented  some  active  SEM  equivalent 
circuits  for  a  dipole  antenna.  Their  networks,  however,  are  intended 
for  computational  purposes  only  and  contain  negative-valued  elements. 

The  foregoing  discussion  can  be  summarized  with  the  observation 
that,  to  date  and  to  the  best  of  the  author's  knowledge,  no  active, 
physically  realizable  equivalent  circuits  valid  over  significant 
frequency  band-widths  have  been  synthesized  and  analyzed  for  electro¬ 
magnetic  energy  collecting  structures. 
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1.3  Scope  of  the  Present  Work 


The  outline  of  this  work  is  as  follows.  In  Chapter  2  the  SEM 
formalism  is  defined  and  the  SEM  expansion  for  the  current  density 
on  a  scattering  object  is  developed.  This  representation  is  derived 
formally  via  the  eigenfunction  expansion.  It  is  based  on  some  evidences 
for  which  a  rigorous  proof  has  not  been  furnished,  to  date.  Some 
questions  concerning  the  formal  development  fo  the  SEM  exapnsion  are 
elaborated  upon  in  Appendix  A.  The.  convergence  properties  of  the  two 
SEM  coupling  coefficient  forms  (termed  class  1  and  class  2),  as  a 
function  of  the  time  origin  location,  are  investigated.  Then  the  most 
convenient  choice  of  the  time  origin  location  for  the  purposes  of  SEM 
equivalent  circuit  synthesis  is  determined.  Also,  the  dominant-pole 
approximation  is  discussed  and  numerical  evidence  is  presented  showing 
that  it  can  be  successfully  used,  at  least  in  the  case  of  highly- 
resonant  structures. 

In  Chapter  3  is  presented  a  formal  development  of  SEM  equivalent 
circuits  for  energy-collecting  structures.  Both  the  short-circuit  and 
the  open-circuit  boundary  value  problems  are  discussed.  With  some 
exceptions,  the  material  of  this  chapter  is  a  recasting  of  the  rf suits 
obtained  previously  by  Beam  [1976a],  in  a  notation  consistent  with 
the  rest  of  this  work. 

In  Chapter  4  the  realizability  conditions  for  the  pole-pair 
immittances  and  voltage  transfer  functions  are  discussed  and  explicit 
realizations  for  them  are  given.  The  real-part  padding  technique  for 
the  synthesis  of  approximate  driving-point  circuits  is  also  addressed, 
and  an  explicit  expression  for  the  amount  of  padding  necessary,  in 
terms  of  the  SEM  parameters,  is  derived.  Most  of  the  material  of 
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this  chapter  is  original  to  the  present  work. 

Chapters  5  and  6  present  an  application  of  the  developed  synthesis 
procedures  to  the  thin-wire  loop  and  dipole  antennas,  respectively. 
Practical,  physically  realizable  active  equivalent  networks  are  presented 
and  analyzed  in  the  time  domain  with  the  help  of  a  circuit  analysis 
program.  The  results  are  compared  with  known  repsonses  or  with  responses 
obtained  by  the  classical  frequency  domain  technique.  In  addition,  in 
Chapter  6  a  sinusoidal  current  mode  approximation  is  used  to  develop  a 
particularly  simple,  and  yet  reasonably  accurate,  equivalent  circuit 
for  a  symmetrical  dipole,  or  a  cylindrical  post  over  a  ground  plane. 

An  application  of  the  equivalent  circuits  developed  in  Chapters  5  and 
6  to  the  analysis  of  antennas  with  nonlinear  loads  is  described  in 
Appendix  B  where  representative  equivalent  circuit  responses  for  a 
thin-wire  loop  and  a  dipole  antenna  loaded  with  a  semiconductor  diode 
are  included.  The  calculation  of  these  responses  with  other  methods, 
though  tractable,  is  quite  tedious. 

Chapter  7  is  devoted  to  the  equivalent  circuit  synthesis  for  a 
spherical  antenna.  The  dominant  pole-pair  circuits  are  synthesized 
by  the  real-part  padding  technique  and  the  quality  of  the  dominant 
pole  approximation  is  examined.  In  Appendix  C  the  static  capacitance 
of  the  spherical  antenna  is  computed  both  analytically  and  numerically. 

The  conclusions  from  this  work  are  summarized  in  Chapter  8. 
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CHAPTER  2 


THE  SEM  REPRESENTATION  FOR  EM  SCATTERING 

2.1  Introduction 

In  Sections  2. 2-2.4  of  this  chapter  the  SO!  formalism  is  defined 
and  the  SEM  expansion  for  the  current  density  on  a  scattering  object 
is  derived  formally  via  the  eigenfunction  expansion.  In  Section  2.5 
the  convergence  properties  of  the  two  SEM  coupling  coefficient  forms 
(termed  class  1  and  class  2)  are  investigated  as  a  function  of  the 
time  origin  location.  Then,  in  Section  2.6,  the  optimal  choice  of  the 
time  origin  location  for  the  purposes  of  SEM  equivalent  circuit  synthesis 
is  determined.  Finally,  in  Section  2.7  the  dominant  pole  approximation 
is  discussed  and  numerical  evidence  is  presented  that  it  can  be  success¬ 
fully  used,  at  least  in  the  case  of  highly-resonant  structures. 

2.2  Integral  Equation  Formulation 

The  integral  equation  for  the  current  density  J(r,s)  on  a  surface 
S  of  a  perfectly-conducting  object  immersed  in  a  lossless,  homogeneous 
medium  may  be  formulated  as 

<F(r,r '  ,s)  ;?(r',s)>  -  ^(r.s)  ,  reS  (2.1) 

where  Ii(r,s)  is  the  known  incident  or  source  field  and  r(r,r',s)  is  the 
dyadic  kernel  of  the  integral  equation.  The  symmetric  product 
notation 

<F(r);  G(r)>  i  jj  ?(r)  •  G(r)dS 
S 

is  used  to  denote  integration  over  the  surface  S  of  the  object. ^ 

*^If  two  or  more  ranges  of  integration  are  used  then  subscripts  can 
be  placed  on  <,>. 
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As  in  [Mann  and  Latham,  1972]  it  is  assumed  that  S  is  a  smooth  surface 
of  finite  extent  enclosing  a  simply-connected  region.  The  solution  of 
(2.1)  can  then  be  sought  in  L2(S),  the  Hilbert  space  of  all  square- 
integrable  functions  with  support  S.  The  frequency  variable  is 
s“0  +  jco  and  is  that  associated  with  the  Laplace  transform  pair*^ 


L{f(t) }  r.  f(s)  -  f(t)e  dt 


(2.2a) 


f(t)  =  L_1{f (s) }  -  2~y  f(s)estds 


(2.2b) 


where  tildes  (~)  over  quantities  distinguish  Laplace  domain  members  of 
transform  pairs  and  C„  is  the  Bromwich  contour. 

D 

This  work  deals  specifically  with  the  electric  field  integral 
equation  (EFIE)  characterization  for  a  conducting  object.  In  this 

2) 

case  the  kernel  of  eq.  (2.1)  becomes  the  symmetric  impedance  kernel 


Z(r,r',s)  »  ZQy 


fl +— T  G  (r,r '  ,s) 

i  Y  J  0 


where 


G  (r.r's,) 


-Y | r-r  ’ 
4tt  I  r-r  ’ 


(2.3) 


(2.4) 


'The  unilateral  Laplace  transform  is  used  since  only  causal  functions 
of  time  are  dealth  with.  This  is  accomplished  by  a  judicious  choice 
of  the  time  origin  of  the  problem  (see  Section  2.5). 

2) 

The  primed  and  unprimed  differential  operators  act,  respectively,  on 
the  primed  and  unprimed  variables.  Strictly  speaking  (2.1)  is  only 
meaningful  when  interpreted  in  terms  of(  some  sense  of  principal  value 
integration.  We  use  it  and  the  kernel  Z  in  a  symbolic  sense  for  efficiency 
of  notation.  Bringing  the  gradient  operator  in  (2.3)  outside  of  the 
integration  leaves  an  integrable  kernel,  and  it  is  upon  this  integro- 
differential  operator  that  we  may  build  a  meaningful  mathematical 
theory. 


i*  -  ' 


is  the  scalar  Green  function,  1  denotes  the  identity  dyadic,  y  -  s/c 
is  the  complex  propagation  constant,  and  Zq  is  the  intrinsic  impedance 
of  the  medium.  The  speed  of  light  is  denoted  by  c.  The  EFIE  incident 
field  function  is  the  tangential  component  of  the  incident  electric 
field  on  S 


I^(r , s)  -  E* (r,s)  •  (T- nn)  (2.5) 

where  n = n(r)  is  a  unit  vector  normal  at  r  to  the  object  surface  S. 

In  most  cases  of  practical  interest  the  incident  field  is  "factorable", 
i.e.,  it  can  be  written  in  the  form 

E*(r, s)  -  EQ(r)p(r,s)f (s)  (2.6) 

-*■  ->■  ^  -► 
where  E^Cr)  indicates  a  spatial  distribution  factor,  p(r,s)  is  a  propa- 

gation  factor  which  links  space  and  time  dependencies,  and  f(s)  is 

the  Laplace  transform  of  the  time  history  of  the  wave.  It  is  assumed 

here  that  f(s)  is  algebraic,  which  is  also  satisfied  in  most  practical 

cases.  For  example,  a  plane  wave  carrying  a  time  history  f(t)  may  be 

written  as 


E^r.s) 


iiYP 

o 


<r'ro>i<«> 


(2.7) 


where  Eq  is  a  constant  polarization  vector,  p  is  a  unit  vector  in  the 
direction  of  propagation,  and  rQ  is  a  vector  pointing  to  a  chosen 
phase  reference  point.  The  time  domain  counterpart  of  (2.7)  is 

?i(r,t)  «  Eo  6  [t  - 

-!of[t-p  •(«„)/«] 

-  E  f(t-  t  -p  •  r/c)  (2.8) 

o  o 


(r-rQ)/c  *f(t) 
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where  6(t)  is  the  Dirac  delta  function,  *  denotes  convolution  in  time, 
and 


t  -  -p  •  r  /c 
o  o 


(2.9) 


It  is  noted  that,  since  the  causality  assumption  is  imposed  on  f(t), 
i.e .j  f(t)  “0  for  t<0,  the  excitation  first  reaches  the  point  r^r^ 
at  the  time  t“0.  Hence,  r^  defines  the  time  origin  of  the  problem: 
the  time  is  counted  beginning  from  the  moment  when  the  wave  front 
passes  the  point  r  ■  rQ.  The  wave  front  reaches  the  coordinate  origin 

r  =•  (5  at  t  =  t  . 

o 

2.3  Formal  Solution  of  the  EFIE 

The  kernel  (2.3)  is  complex-symmetric  [Cochran,  1972]  and  so  one 
may  define  a  set  of  eigenvalue/eigenfunction  pairs  through 

<Jn(r',s);Y(r' ,r,s)>  -  <Z  (r,r '  ,s) ;  J(r',s)>  *  Xn(s)Jn(r,s)  .  (2.10) 

It  follows,  formally,  that  the  kernel  and  its  resolvent  can  be  expanded 


~  ^  „  J  (r,s)  J  (r',s) 

Z(r,r',s)  =  Ms>  X  - - - 


(2.11) 


<J  (r,s) ;J  (r,s)> 
n  n 


-lr— .  .  r  1  Jn(r,s)  Jn(7',s) 

(r,r  ,s)  =■  I  - -r—; - 

n  <J  (r,s) ;J  (r,s)> 


The  solution  to  the  integral  equation  (2.1)  is  thus  written 


(2.12) 


I 


J(r,s)  -  <Z  1  (?,r',s);I1(?',s)> 

_  i  <J  (r’,s);Ei(r’,s)>  ~ 

-  I  rh  -5-TT - -  ?n(r.s)  •  (2.13) 

n  ^n1,  <J  (r,s);J  (r,s)> 
n  n 

The  EEM  representation  (2.13)  is  formal  in  the  sense  that  in 

deriving  it  some  tacit  assumptions  have  been  made  which,  in  general, 

may  not  be  satisfied.  Among  other  things,  the  major  pitfall  in  the 

expansions  above  is  that  since  the  kernel  (2.3)  is  not  Hermit iar.-symmetric , 

there  is  nothing  to  guarantee  the  existence  of  a  complete  set  of  eigen- 

vectors  Jn(r,s),  as  defined  in  (2.10).  One  cannot  even  be  sure  that 

the  normalization  factors  in  the  denominators  above  do  not  vanish 

[Farm,  1973] .  Perhaps  the  most  common  situation  in  practice  in  which 

the  eigenfunction  representation  fails  is  when  the  eigenvalues  A  (sj. 

n 

have  branch  points.  At  a  branch  point  in  the  s-plane  two  or  more 
eigenvalues  coalesce  [Coehran,  1972],  so  that  the  inclusion  of  root 
vectors1^  in  the  expansion  may  be  necessary  [Ream,  1980].  In  some 
cases  even  for  degenerate  eigenvalues  a  set  of  independent  eigenvectors 


;Let  A  be  a  linear  operator  on  the  space  considered  and  let  i  # 0 
be  an  eigenvector  associated  with  an  eigenvalue  Ait  i.e.,  (A-lAi) 
$i,l*0,  where  I  denotes  the  identity  operator.  Consider  the 
subspace  of  all  vectors  that  are  annihilated  by  some  power  of 
(A-I<\i).  If  for  some  power  k  there  exists  a  vector  4>i,k  the 

subspace  such  that  (A-lA^)*1-!^ ^  /  0  and  (A-IAi)k<J>^ ^  ■  0,  and  such 
that  there  is  no  vector  in  the  subspace  that  is  not  annihilated  at 
least  by  (A-lA^)^,  then  <&ifk  is  a  root  vector  (generalized  eigenvector 
[Pease,  1965])  of  rank  k,  where  k  is  the  length  of  the  longest  chain 
associated  with  eigenvalue  A^.  The  chain  {4>i,l»4>i,2» •  •  •  »4>i,k) *  with 
the  elements  defined  by  <f>i,k-l  ■  (A-lAi)$itk, . . .  ,<Ji,l  *  (A-lAi)k~A<{i^t k, 
is  called  a  Jordan  chain  of  length  k.  The  union  of  all  root  vectors 
of  a  linear  operator  A  corresponding  to  all  eigenvalues  of  A  is  called 
the  root  system  of  A  [Ram,  1980]. 
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can  be  found  so  that  the  root  system  coincides  with  the  eigensystem. ^ 

2) 

This  is  true  for  normal  operators.  For  example,  it  was  shown  by 
Ream  [1973]  that  the  scalar  analog  of  the  operator  in  (2.1)  is  normal 
if  S  is  a  sphere.  If  the  operator  is  not  normal,  it  is  not  always 
diagonalizable  and  considerable  complexity  arises  in  the  computational 
scheme  [Dolph  et  al.,  19807. 

There  is  an  important  class  of  bodies  for  which  it  can  be  shown 

that  branch  points  do  not  occur.  These  are  structures  such  as  the 

sphere  or  circular  loop,  in  which  the  eigenvectors  are  completely 

determined  by  the  geometrical  symmetry  of  the  structure  and  hence  are 

frequency  independent.  This  follows  from  the  Rayleigh  quotient 

3) 

representation  of  the  eigenvalue, 

<J  (r);Z(r,r',s);J  (r')> 

Vs)  -  *  ■»  -S  V  -  •  (2.14) 

<J  (r);J  (r)> 
n  n 

The  kernel  in  (2.14)  can  be  expanded  in  a  Laurent  series  converging 
for  all  | s |  >0,  and  integrated  term-by-term  to  obtain  a  series  with 
the  same  region  of  convergence  and  analyticity.  This  argument  is  not 
valid,  however,  if  the  body  is  not  completely  symmetric  because  the 
eigenmodes  Jn  are  then  generally  functions  of  s.  To  date,  analytic 

^In  the  language  of  finite-dimensional  spaces  this  means  that  the 
geometric  multiplicity  of  an  eigenvalue  is  equal  to  its  algebraic 
multiplicity  [Stakgold,  1967],  We  say  that  an  eigenvalue  X*  has 
algebraic  multiplicity  k  if  the  term  (X  -  X^)  is  repeated  k  times 
in  the  characteristic  equation  of  the  operator.  Geometric  multiplicity 
of  Xj  denotes  the  number  of  independent  eigenvalues  corresponding  to  it. 

2)  3d  3 

An  operator  A  is  called  normal  if  A  A  ■  AA  ,  where  A  is  the  adjoint 

operator. 

3) 

The  triple  product  notation  is  introduced  for  a  triad  of  elements, 
the  middle  element  of  which  is  dyadic:  <  ;  ;>■<;<;». 
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properties  of  the  eigenvalues  of  only  rotationally-symmetric  objects, 
such  as  the  sphere  and  the  circular  loop,  have  been  thoroughly  examined. 


The  singularity  expansion  for  the  surface  current  density  on  a 
scattering  object  can  be  obtained  from  the  eigenfunction  expansion  by 
performing  the  residue  series  expansion  of  (2.12)  or  (2.13)  in  terms 
of  the  complex  natural  resonances  (the  poles)  of  the  object,  defined 
as  follows: 

(s  4:  X  (s  .)  *  ol  .  (2.15 

(  ni  n  ni 

It  is  noted  that  for  s-s  .,  in  view  of  (2.10)  and  (2.15),  the  following 

ni 

is  true: 


<Z(r,r’,sn.);  Jn(r',sn.)> 


(2.16) 


— ►  — ►  -►  -V 

Hence,  J  (r,s  ,)=J  .  (r)  are  the  natural  modes  of  the  object  associated 
n  ni  ni 

with  the  poles  s^.  In  practice  the  homogeneous  integral  equation 
(2.16)  is  replaced  by  a  matrix  equation  by  some  approximate  technique 
such  as  the  method  of  moments  (MoM)  [Harrington,  1968],  and  the  poles 
are  sought  as  zeros  of  the  determinant  of  the  coefficient  matrix. 

Before  expanding  the  reciprocal  eigenvalue  factor  which  appears  in 


(2.12)  and  (2.13)  in  terms  of  its  singularities  in  the  complex 
frequency  plane,  it  is  expedient  to  subtract  from  ^^(s)  the  possible 
pole  at  the  origin  and  define 


h  (s)  ■  ■* - 

n  Xn<9) 


(2.17) 


where  r  denotes  the  residue  of  the  pole  at  s»0.  Assuming  first 
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Hi?*,!-- 


order  poles  one  can  postulate  uaing  the  Mittag-Lef f ler 
expansion  [Jeffreys  and  < Jeffreys ,  1956]  the  representation 
{ Pearson  and  Wilton,  1980;  Pearson,  1981] 


h  (s)  »  h 
n 


( o )  +  £  r  .  -  +  -  +  b  ( 

n  “  m  s-s  .  s  .  n 

i  ni  ni 


(2.18) 


where  b  (s)  denote  the  branch-integral  contributions  which  mav  occur 
n 

and  r  are  the  residues  which  can  be  computed  by  using  eq .  (2.14) 
as  well  as  the  orthogonality  property  of  the  eigenfunctions  to  yield 


r  =*  Pj—  \  (s  )~|  =3  ,<J  (r) ; J  .  (r)> 

ni  [ds  n  ni  j  ni  ni  ’  ni 

with  the  normalization  constants  ft  given  by 
6nl  *  fj  ?<?,;•, 


(2.19) 


(2.20) 


Marin  a>id  Latham  [1972]  have  shown  that  for  perfectly  < , nducting 
bodies  of  finite  extent  imbedded  in  a  lossless  medium *Z  1  (r.r's)  is 
analytic  in  s  except  for  pole  singularities.  This  does  not  mean,  however, 
that  there  cannot  be  branch- integral  contributions  to  (2.18)  which 
cancel  out  in  the  summation  (2.14),  i.e.. 


„  ~  <J  (r,s);E  (r)p(r,s)> 

y  b  (S)  —2 - ° - 

n  <J  (r,s);J  (r,s)> 

n  n 


J  (r,s)  =  0  . 
n 


(2.21) 


It  is  demostrated  in  Appendix  A  that  such  phenomena  can,  indeed,  occur. 
Only  for  bodies  such  as  the  sphere  or  the  circular  loop  antenna,  in 
which  geometrical  symmetry  completely  determines  the  eigenfunctions 

x^It  was  conjectured  by  Beam  [1971]  that  perfectly  conducting,  finite- 
extent  oodies  in  homogeneous,  lossless  media  possess  only  simple  (first 
order)  poles.  A  proof  of  this  conjecture  has  been  put  forth  recently 
by  Sanoer  and  Varvatsis  [1980].  The  simple  pole  restriction  exludes 
also  essential  singularities  in  a  finite  complex  plane  since  they  can 
be  regarded  as  infinite  order  poles. 
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(and  hence  they  are  frequency  independent),  has  it  been  possible  to 
show  (Section  2.3)  that  branch  points  cannot  occur,  and  poles  are  the 
only  singularities  of  the  reciprocal  eigenvalue  in  the  finite  complex 
frequency  plane,  so  that  bn(s)*0.  On  the  other  hand,  any  losses  in 
the  medium  always  contribute  a  branch  cut  in  the  complex  frequency 
plane.  Although  in  the  lossless  medium  case  the  branch-integral 
contributions  ultimately  cancel  leaving  the  singularity  expansion  of 
current  unaffected,  they  are  likely  to  have  bearing  on  SEM  formulations 
where  eigenmode-by-eigenmode  expansion  enters  in,  including  the  SEM 
equivalent  circuit  synthesis  discussed  in  this  work.  Pearson  and 
Wilton  [1981]  have  shown  that  the  reciprocal  eigenvalues  in  (2.13) 
are  positive  real  (PR)^  functions  of  s;  they  have  also  speculated 
that  the  same  may  hold  for  the  terminal  admittances  associated  with 
individual  current  eigenmode  contributions  to  total  current.  The  PRness 
is  a  necessary  condition  for  the  realizability  of  driving-point  functions. 
If  the  branch-integral  constituent  in  (2.18)  is  non-zero,  its  ommission 
may  lead  to  non-realizable  immittances.  This  difficulty  can  potentially 
be  circumvented  by  grouping  terminal  eigenadmittances  which  share 
common  branch  points,  so  that  the  branch-integral  constituents  add 
to  zero  and  PRness  is  preserved.  Since  in  most  practical  cases  the 
branch- integral  terms  cannot  be  explicitly  identified,  the  t>n(s) 
factors  in  (2.18)  will  be  suppressed  throughout  the  rest  of  this  work. 

The  expansion  (2. 1C)  is  valid  provided  hn(s)  is  bounded  on  a  set 
of  circular  contours  centered  at  s*0  and  passing  between  the  poles 
[Jeffreys  and  Jeffreys ,  1956].  If  hn(s)  is  not  only  bounded,  but 
actually  decays  on  this  set  of  contours,  one  can  show  that  [Pearson,  1981] 

^See  Section  4.2  for  the  definition  of  a  positive  real  function. 


(2.22) 


and 


(2.23) 


It  is  noted  that  even  when  bn(s)*0  in  (2.18),  the  poles  and  the 

residues  alone  do  not  determine  h  (s)  and  some  additional  information 

n 

is  necessary  to  determine  the  constant  term  h  (0).  Often  it  can  be 

n 

argued  from  the  high  freuqency  and/or  low  frequency  considerations 

that  h  (0)  *  0.  In  that  case  it  follows 
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Similarly,  using  (2.17)  and  (2.23)  yields 


(2.24) 


(2.25) 


Following  the  terminology  introduced  by  Beam  [1976a],  the  terms  in 
brackets  in  (2.24)  and  (2.25)  are  referred  to  as,  respectively,  the 
modified  and  unmodifeid  pole  modules.  Both  forms  have  been  used  in 
the  SEM  equivalent  circuit  synthesis  to  date.  Throughout  most  of 
this  work,  the  modified  pole  modules  are  used  since  they  impose 
less  stringent  conditions  on  the  behavior  of  Xn(s).  It  is  of  interest 
to  point  out,  that  the  unmodified  pole  modules  expansion  (2.25)  was 
proved  by  Umashankar  and  Wilton  [1974]  to  be  valid  for  the  circular 
loop  in  a  lossless  medium.  However,  in  the  actual  computations  the 
truncated  series  (2.25)  converged  very  slowly  at  low  frequencies  and 
they  were  forced  to  use,  in  effect,  the  modified  from  (2,24). 
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The  substitution  of  the  expansion  (2.24)  into  (2.13)  yields 


J(?,s)  -  4  B  fi  (s)j 


s  oo  oo 


l 1 


B_.n. 


.  s  .  (s  -  s  . )  ni  ni 
n  i  ni  ni 


with  the  coupling  coefficients  rj  .  (s) 

ni 


oo(*)?(s) 
(s)Jni(r)f(s) 
defined  as 


(2.26) 


rj(1)  -  <J  .  (r) ; E  (r)p(r,s  .)> 
ni  ni  ’  o  r  ni 


—  the  class  1  coupling  coefficient,  or 


n^(s)  “  <5ni(r);Eo(r)p(r,s)> 


(2.27) 


(2.28) 


—  the  class  2  coupling  coefficient1^  [Baum,  1972;  Bccum,  1976b].  It 

is  noted  that  the  class  1  form  is  simply  a  specialization  of  the  class  2 

coupling  coefficient:  *  fj^(s  .).  Also,  both  forms  are  identical 

ni  ni  ni 

for  nonpropagating  excitation  fields,  when  p(r,s)  *  1.  For  the  plane 
wave  incident  field  the  two  classes  of  coupling  coefficients  take  the 
form 
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nni  "e  <Jni(r);Eoe 
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(2.30) 


The  coupling  coefficients  defined  here  differ  slightly  in  form  from 
the  ones  originally  introduced  by  Baum  [1976b].  In  Baum's  work  the  time 
origin  coincides  with  the  space  origin  but  the  SEM  series  can  be 
"turned  on"  in  negative  time,  if  necessary,  by  choosing  the  value  of  the 
turn-on  time  parameter.  In  the  present  formulation,  which  follows 
naturally  from  the  assumed  form  of  the  incident  field  (27),  the  time 
origin  rQ  is  chosen  independently  from  the  space  origin  and  the  SEM 
series  is  always  turned  on  at  t-0. 


In  the  next  sections  the  convergence  properties  of  both  coupling 
coefficient  forms  are  discussed  and  what  is  believed  to  be  the 
best  choice  of  the  coupling  coefficient  for  the  purposes  of  SEM  circuit 
synthesis  is  determined.  Meanwhile  the  generic  notation  hni(s)  is 
used  to  embrace  both  forms. 

Several  comments  are  in  order  concerning  eq.  (2.26).  First,  the 

term  associated  with  the  pole  at  the  origin  is  shown  explicitly  and 

it  is  assumed  that  this  pole  belongs  to  the  zeroth  eigenmode.  It  can 

be  argued  from  the  form  of  the  integral  equation  that  this  mode  must 

be  divergenceless  (solenoidal) .  If  the  structure  is  not  capable  of 

supporting  a  solenoidal  current  mode,  it  can  not  have  a  pole  at  the 

origin.  Second,  it  is  noted  that  a  concellation  of  a  pole  can  occur 

if  it  coincides  with  a  zero  of  the  corresponding  coupling  coefficient. 

A  limiting  procedure  should  be  used  in  such  case  to  determine  the 

contribution  of  such  a  pole-zero  pair.  Finally,  it  should  be  noted  that 

perhaps  the  most  improtant  feature  of  the  representation  given  in 

(2.26)  is  that  the  set  of  complex  natural  frequencies  (s  },  the  set  of 

n 

complex  natural  modes  (Jn^(r)},  and  the  set  of  complex  normalization 
coefficients  {8^1  are  dependent  on  the  object  parameters  only  but  are 
indept  ient  on  the  excitation.  The  effect  of  the  exciting  field  is 
contained  entirely  within  the  set  of  coupling  coefficients  (fjn;l(s)}. 

2.5  Convergence  Properties  of  Class  1  and  Class  2  Coupling  Coefficients 

In  this  section  are  investigated  the  convergence  properties  as 
Re{sni}-*-°°  of  the  singularity  expansion  (2.26)  with  both 
coupling  coefficient  forms,  as  a  function  of  the  choice  of  the  location 
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of  the  time  origin  rQ.  The  tools  for  this  analysis  were  provided  in 
a  recent  paper  by  Baxan  and  Pearson  [1981].  It  is  helpful  to  follow 
them  in  defining  the  following  characteristic  dimensions  for  use 
subsequently:  the  largest  linear  dimension  of  the  scatterer 
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the  "clearance  distance 


(2.31) 
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and  the  "contact  distance" 
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The  latter  two  distances  define,  respectively,  the  signed  distance  from 
the  coordinate  origin  to  the  point  at  which  the  wave  front  clears  the 
object,  and  the  signed  distance  from  the  coordinate  origin  to  the 
point  at  which  the  wave  front  first  contacts  the  object  (see  Figure  2.1). 

To  investigate  the  convergence  properties  of  the  SEM  expansion 
(2.26)  the  asymptotic  properties  of  the  terras  B^  ?i  (s)  need  to  be 
examined.  The  behavior  of  B  ^  for  Re{sn^}-^-°°  was  argued  by  Baum  and 
Pearson  [1981]  to  be 
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The  behavior  of  the  class  1  coupling  coefficient  is  governed  by 
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Figure  2.1  Portrayal  of  a  scatterer  in  a  plane-wave  field.  Indicated 
arc  the  maximum  extent  of  the  object  Lc.  the  clearance  distance  Lx 
relative  to  the  coordinate  origin,  and  the  contact  distance  L?. 


This  term  decays  provided 
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This  condition  is  most  stringent  when  the  direction  of  propagation  of 
the  incident  plane  wave  is  along  the  largest  dimension  of  the  object. 
In  this  case  (2.37)  reduces  to 


P 


(2.38) 


Hence,  if  the  time  origin  is  placed  on  the  surface  of  the  scattering 
object  the  convergence  of  the  class  1  SEM  expansion  term  is  assured. 

Similar  analysis  of  the  class  2  coupling  coefficient  form  shows 


that 


B  ~(2)  .  sniLo/c 

6ni  nni  <S)  *  6 


(2.39) 


Since  Lq  is  always  positive,  one  can  conclude  from  the  last  expression 
that  the  class  2  form  is  always  convergent,  independet  of  the  choice 
of  the  time  origin. 


2.6  Optimal  Choice  of  the  Time  Origin 

In  this  section  are  discussed  the  numerical  properties  of  the  two 
classes  of  coupling  coefficients  in  the  case  of  plane  wave  incident 
field.  Also,  the  most  convenient  choice  of  the  time  origin  rQ  for  SEM 
circuit  synthesis  if  determined. 

Although  both  class  1  and  class  2  coupling  coefficients  give  valid 
representation  for  the  surface  current  density  when  used  in  (2.26)  , 
they  exhibit  quite  different  properties  from  the  practical  computation 
point  of  view.  This  is  evident  when  one  performs  the  Laplance  inversion 
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of  (2.26)  to  obtain  the  time  domain  response 


J(r,t)  -  L-1{J(r,s)}  . 


(2.40) 


In  this  process  terms  like 
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must  be  evaluated.  Using  the  class  1  coupling  coefficient  (2.29)  in 
(2.41)  gives 
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whereas  the  class  2  coupling  coefficient  (2.30)  gives 
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Thus,  in  the  time  domin  the  class  1  and  class  2  coupling  coefficients 
take,  respectively,  the  form 


^It  should  be  kept  in  mind  that  any  singularities  of  f(s)  must  be 
taken  into  account  in  this  process. 


(2.44) 
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It  should  be  noted  that  the  Integrand  in  (2.45)  gives  a  nonzero 
contribution  onlv  when 


or,  equivalently,  when 


p  *  r/c  <  t  -t 


p  *  (r  -  r  )/c  <  t 
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These  expressions  can  be  interpreted  with  the  help  of  Figure  2.1  as 
follows.  The  pole  terms  in  (2.40)  are  "turned  on"  at  t  =  0, 
as  sl>.Ii  as  the  wave  front  passes  the  time  origin  r  =  rQ.  However,  t *  - y 
are  computed  differently  depending  on  whether  class  1  or  class  2  cour>l ing 
coetticients  are  used.  The  class  1  coupling  coefficients  are  constants 
resulting  from  the  integration  over  the  entire  object,  as  indicated  in 
(2.44).  The  class  2  coefficients  (2.45),  however,  are  obtained  by 
integrating  only  over  that  part  of  the  object  illuminated  at  the  given 
moment  by  the  incident  plane  wave.  Thus,  they  are  time  dependent 
until  t*tc>  when  the  edge  of  the  incident  wave  has  cleared  the  object. 
After  that  time  they  take  on  constant  values  identical  with  those  of 
the  class  1  coefficients.  It  is  evident  from  the  above  discussion, 
that  the  class  2  coupling  coefficients  form  is  much  more  expensive 
computationally  than  the  class  1  form.  It  is  also  unsuitable  for  the 
purposes  of  SEM  circuit  synthesis  because  of  its  s-dependence  in  the 
frequency  domain.  It  has,  however,  some  advantages,  as  discussed  below. 
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It  should  be  noted  that  the  series  (2.40)  must  sum  to  zero  in  the 

period  of  time  before  the  leading  edge  of  the  incident  wave  arrives 

at  the  observation  point  (port  location).  However,  the  terms  of  the 

series  contain  the  factors  exp|sn^[t-p  •  (r-ro)/c]j  which  take  on  large 

values  in  early  times,  because  Re{sni) <  0.  Since  these  large  numbers 

must  sum  to  zero,  the  result  is  sensitive  to  errors  in  the  pole  values 

[Marin,  1972;  Pearson,  1979].  The  class  1  form  requires  also  that  a 

possibly  full  collection  of  poles  belonging  to  a  given  eigenmode  be 

included  in  the  time  domain  counterpart  of  (2.26)  [Pearson,  1979]. 

To  illustrate  this  point,  in  Figure  2.2  are  shown  current  responses 

at  <{>  »  180°  of  a  loop  antenna  of  radius  b  illuminated  by  a  step  function 

plane  wave  polarized  as  indicated  in  the  inset,  computed  in  three  ways. 

The  first  two  curves  were  obtained  by  using  class  1  coupling  coefficients, 

but  different  collections  of  poles.  In  one  case  (solid  line)  all  poles 

belonging  to  eigenmodes  0-10  for  which  Im{s  .}  <  lOc/b  were  included, 

ni 

while  in  the  second  case  (dotted  line)  only  one,  "dominant"  pole  (the 
pole  closest  to  the  jw  axis)  from  each  eigenmode  collection  was  taken. 

The  third  curve  (dashed  line)  was  obtained  by  using  class  2  coupling 
coefficients  and  dominant  poles  only.  The  time  origin  was  chosen  in 
all  cases  at  0-0°,  the  point  of  the  first  contact  with  the  incident  field. 
As  can  be  seen  in  the  figure,  using  class  1  coupling  coefficient  form 
and  incomplete  collection  of  poles  results  in  a  non-causal  behavior 
of  the  response.  The  class  2  form  is  more  tolerant  and  gives  accurate 
responses  for  this  case.  Since  in  most  practical  situations  the  full 
eigenmode  collections  of  poles  are  not  at  one's  disposal,  this  last 
example  would  suggest  that  the  class  1  coupling  coefficient  is  of  no 
value  and  that  one  must  resort  to  the  class  2  form.  It  should  be  noted. 
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however,  that  the  three  curves  in  Figure  2.2  agree  favorably  for  t<2b/c 
after  the  wave  first  arrives  at  the  observation  point.  The  response  in 
the  period  0<t<2b/c  is  not  of  much  interest,  since  it  is  known  a  priori 
to  be  zero.  Why  not  then  move  the  time  origin  to  the  observation 
coordinate  (port  zone)  and  use  the  simple  class  1  coupling  coefficient? 

This  point  is  supported  by  further  examples  presented  below. 

In  Figure  2.3  and  2.4  are  compared  current  responses  of  a  loop 
antenna  computed  at  <f>  =  90 0  and  <£  =0°,  respectively,  by  using  class 
l  (dotted  line)  and  class  2  (solid  line)  coupling  coefficients  and 
dominant  poles  satisfying  Imfs^]  <  lOc/'b.  The  excitation  conditions 
are  identical  to  that  of  the  previous  example.  Again,  one  observes 
good  agreement  of  the  responses  beginning  from  the  time  when  the  wave 
front  reaches  the  observation  point.  It  is  seen,  however,  that  the 
class  1  responses  deteriorate  slightly  for  observation  points  closer 
to  the  point  of  first  contact  of  the  plane  wave  with  the  loop.  This 
can  be  attributed  to  the  fact  that  in  the  process  of  moving  the  observation 
point  toward  the  point  first  illuminated  by  the  incident  field  increases 
the  area  of  the  object  surface  which  in  early  time  lies  ahead  of  the 
wave  front  in  the  integral  (2.44). 

The  conclusions  drawn  from  the  loop  example  also  appear  to  be 
valid  in  a  thin  cylinder  case,  as  demonstrated  in  Figure  2.5-*2.7. 

The  incident  field  is  a  step  function  plane  wave  whose  angle  of  incidence 
i-' =  30° .  In  the  numerical  procedure  the  length  l  of  the  cylinder  was 
divided  into  72  equal-size  zones  and  the  time  origin  was  placed  at  the 
72nd  zone — the  one  first  illuminated  by  the  incident  field.  In  Figures 
2.5,  2.6,  and  2.7  are  shown,  respectively,  the  current  responses 
evaluated  at  zones  18(1/4)1),  36(l/2£),  and  54(3/4H).  Only  the 
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Figure  2.3  Short-circuit  current  excited  in  the  circular  loop  at  <{>  =  90°  by  a  step  function 

plane  wave.  Time  domain  SEM  expansions  using  class  1  and  class  2  coupling  coefficients 
and  dominant  poles  are  compared. 


Figure  2.4  Short-circuit  current  excited  in  the  circular  loop  at  4>=0°by  a  step  function 
plane  wave.  Time  domain  SEM  expansions  using  class  1  and  class  2  coupling  coefficients 
and  dominant  poles  are  compared. 


Figure  2.6  Short-circuit  current  excited  in  the  thin  cylinder  at  z  =  1/28.  by  a  step  function 
plane  wave.  Time  domain  SEM  expansions  using  class  1  and  class  2  coupling  coefficients  and 
dominant  poles  are  compared. 


dominant  poles  for  which  Im[s^]  <10ttc/?,  were  used  in  the  SEM  expansion 
for  the  current.  It  is  seen  from  these  figures  that  class  1  responses 
(dashed  lines)  agree  favorably  with  class  2  responses  (solid  lines) 
beginning  from  the  time  when  the  leading  edge  of  the  plane  wave  reaches 
the  observation  point.  As  in  the  loop  case  the  agreement  deteriorates 
somewhat  when  the  observation  point  is  moved  toward  the  point  of  first 
contact  of  the  incoming  wave  with  the  cylinder.  The  discrepancy 
between  the  class  1  and  class  2  responses  in  the  present  case  is  greater 
than  in  the  loop  case,  which  can  be  attributed  to  the  fact  that  the  poles 
of  the  cylinder  are  computed  less  accurately  than  the  poles  of  the  loop, 
and  to  the  sensitivity  of  the  class  1  coupling  coefficient  form  to 
errors  in  pole  values  [Baum  and  Pearson ,  1981] . 

The  following  observations  can  be  made  from  the  above  examples. 

1.  Using  an  incomplete  pole  collection  in  the  SEM  expansion 
together  with  class  1  coupling  coefficients  gives  non- 
causal  responses  in  the  period  of  time  before  the  leading 
edge  of  the  incident  field  reaches  the  observation  point. 

2.  The  class  2  coupling  coefficient  form  gives  correct  (causal) 

responses  in  a  similar  situation.  It  is  reasonable  to  assume 

that  it  gives  also  more  accurate  responses  than  the  class  1 

form  for  intermediate  times,  before  t«t  . 

c 

3.  Beginning  from  the  time  when  the  incident  wave  first  reaches 
the  observation  point  the  class  1  and  class  2  responses 
agree  very  favorably,  so  that  the  class  1  coupling  coefficient 
can  be  used  if  the  time  origin  is  placed  at  the  observation 
point. 

4.  The  class  1  response  becomes  less  accurate  when  the  observation 
point  (port  zone)  is  moved  closer  to  the  point  first  contact 
of  the  incident  field  with  the  object.  However,  the  deteriora¬ 
tion  of  the  response  is  not  significant  (of.  figures  2.4  and 
2.7). 

5.  The  responses  obtained  by  using  dominant  poles  only  agree 
favorably  with  the  responses  obtained  by  including  more  complete 
collections  of  poles.  (This  fact  is  exploited  in  the  next 
section  where  the  dominant  pole  approximation  is  discussed  in 
more  detail.) 
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From  the  above  discussion  we  conclude  that  the  class  1  coupling 
coefficient  with  the  time  origin  at  the  port  location  is  the  best 
choice  for  the  SEM  equivalent  circuit  synthesis.  In  Baum's  [1976b] 
formulation  this  is  equivalent  to  choosing  the  turn-on  time  at  the 
instant  when  the  wave  front  reaches  the  observation  point. 

2.7  Dominant  Pole  Approximation 

It  is  noted  that  the  singularity  expansion  for  the  surface  current 
density  on  the  object  (2.26)  contains  a  potentailly  double  infinite 
summation  of  terms  associated  with  an  infinite  number  of  eigenmodes 
and  potentially  infinite  collection  of  poles  belonging  to  a  given  eigen¬ 
mode  [Pearson  and  Wilton,  1981] .  For  many  problems  of  practical  interest 
one  can  recover  by  numerical  methods  [Tesahe,  1973;  Crow  et  al 1972; 
Singaraju  et  al.s  1976],  analytical  methods  [VanBlariaum  and  Mittra , 

1975;  Pearson  and  Roberson,  1980;  Pearson  and  Lee,  1982;  C>.o  and 
Cordaro ,  1979]  a  finite  number  of  poles  which  dominate  the  frequency 
and  time  responses  for  excitations  whose  spectra  are  bandlimited. 

These  poles  are  the  ones  located  closest  to  the  jio-axis  and  are  referred 
to  as  the  dominant  poles  of  the  structure  [Pearson  and  Wilton,  1981]. 

In  the  case  of  the  few  structures  for  which  the  distribution  of  poles 
among  the  eigenmodes  can  be  explicitly  determined  ( e.g the  circular 
loop  and  the  sphere)  the  collection  of  dominant  poles  can  be  obtained 
by  taking  only  one  "dominant"  pole  from  the  number  of  poles  belonging 
to  each  eigenmode.  There  are  strong  indications  that  this  is  true 
in  general.  Denoting  by  sr  the  dominant  pole  of  the  nth  eigenmode, 
one  can  rewrite  (2.26)  in  the  dominant  pole  approximation  as 


(2.48) 


J(v,s)  -  \  Bono(s)JQ(r)f(s )+l  -s  (s3_s  }  Bnnn(s)Jn(r)f(s) 


n  n  n 


with 


fln(s)  -  <Jn(r);Eo(r)p(r,s)>  . 


(2.49) 


It  is  understood  in  (2.48)  that  complex  conjegate  poles  are  included 
in  the  summation.  It  is  noted  that  the  term  associated  with  the  pole 
at  the  origin  is  shown  explicitly  in  (2.48).  As  pointed  out  in  Section 
2.3,  this  term  is  present  only  in  the  case  of  multiply-connected  objects, 
such  as  the  loop  antenna.  It  is  absent,  for  example,  in  the  case  of 
the  thin  dipole  antenna,  which  can  be  considered  as  a  member  of  the 
singly-connected  class  of  structures. 

It  is  important  to  note  that  in  the  case  of  a  multiply-connected 
object  immersed  in  an  external  incident  filed,  one  expects  a  cancellation 
of  the  pole  at  the  origin  to  occur.  This  is  so,  because  in  the  limit 
as  s+0  the  coupling  coefficient  ^Q(s)  is  obtained  by  integrating  along 
a  closed  path  a  conservative  (in  the  limit)  field  against  the  instant 
natural  mode.  One  may  rewrite  (2.48)  in  such  cases  as- 


J (r,s)  -  uo8oJo(r)f (s)  +  £ 

n 


s  (s  -  s  ) 
n  n 


0T,n„(s)J  (r)f (s) 

n  n  n 


(2.50) 


where 


U  *  &im 
0  s-O 


fno(s) 


(2.51) 


For  example,  in  the  plane  wave  case,  in  the  limit  as  s-*-0,  one  arrives 
at 
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+  +  -YP  *  (r  -  r  ) 
fjo(s)  -  <J0(r);Eoe  > 

3  <Jo(0;Eo[l-yp  •  (r-ro)]> 

-  -s<JQ(r) ;Eqp  •  (r-rQ)/c>  (2.52) 

so  that 

UQ  -  -<Jo(?);EQp  •  (Z-rQ)/c>  .  (2.53) 

2.8  Conclusions 

The  SEM  representation  for  the  surface  current  density  on  a  scattering 
object  is  derived  formally  by  way  of  the  eigenfunction  expansion,  and 
the  potenital  limitations  of  this  procedure  are  addressed.  Since  the 
expansion  is  performed  on  the  eigenmode-by-eigenmode  basis,  there  may 
be  branch-integral  contributions  present  [Pearson  and  Wilton,  1981] 
which,  in  the  case  of  finite-extent  perfectly-conducting  objects  in 
lossless  media,  must  cancel  in  the  eigenmode  summation,  to  give  a 
meromorphic  current  response  in  accord  with  [Marin  and  Latham,  1972]. 

It  is  demonstrated  by  specific  examples  in  Appendix  A  that  such  phenomenon 
can,  indeed,  occur.  Only  for  highly  symmetric  objects,  having  frequency 
independent  eigenfunctions,  has  it  been  possible  to  show  (Section  2.3) 
that  the  eigenvalues  are  meromorphic  in  s.  The  omission  of  the  branch- 
integral  constituent  in  the  SEM  exapnsion  can  adversely  affect  the 
realizability  of  the  terminal  eigenadmittances  [Pearson  and  Wilton, 

1981]. 

Two  basic  SEM  expansion  forms  arise  in  the  SEM  expansion  of  the 
current,  the  so-called  modified  and  unmodified  forms.  Since  the  modified 
form  imposes  less  stringent  conditions  on  the  large-s  asymptotic 
properties  of  the  inverse  eigenvalues  [Pearson,  1981],  it  is  suggested 
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that  the  modified  poles  modules  are  preferable  in  the  case  of  problems 
dealth  with  on  a  purely  numerical  basis. 

Two  forms  of  coupling  coefficients  have  been  used  to  date,  the 
so-called  class  1  and  class  2  forms  [Baum,  1976b] .  For  the  SEM 
equivalent  circuit  synthesis  one  is  forced  to  use  the  class  1  SEM 
expansion  form.  It  is  demonstrated  that  this  form  gives  good  results 
provided  the  time  origin  is  placed  at  the  gap  zone.  The  convergence 
of  this  representation  is  assured  since  it  is  shown  that  the  class  1 
expansion  form  is  always  convergent  as  Refs^}-*00  if  the  time  origin 
is  located  on  the  surface  of  the  scattering  object. 

It  is  also  demonstrated  that  for  the  thin-wire  loop  antenna  the 
dominant-pole  approximation  [Pearson  and  Wilton,  1980]  gives  very 
accurate  short-circuit  current  responses.  It  is  conjectured  that  this 
is  true  for  all  highly  resonant  structures. 
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CHAPTER  3 


FORMAL  DEVELOPMENT  OF  EQUIVALENT  CIRCUITS 

3.1  Introduction 

In  this  chapter  is  discussed  a  formal  development^  of  SEM  equivalent 

circuits  for  energy  collecting  structures  with  definable  ports.  With 

some  exceptions,  this  development  essentially  follows  that  presented 

by  Baum  [1976a]  and,  in  a  more  recent  paper,  by  Pearson  and  Wilton 

[1980].  It  is  reiterated  here,  for  the  sake  of  completeness,  in 

notation  consistent  with  the  rest  of  this  work. 

The  objective  of  this  work  is  to  derive  an  equivalent  circuit  at 

2) 

a  port  of  an  antenna  or  scatterer  from  its  SEM  description.  Since 

the  problem  is  linear,  one  can  construct  the  Norton  and  Thevenin 

equivalent  networks,  as  shown  in  Figure  3.1.  It  is  noted  that  the  Norton 

circuit  requires  the  knowledge  of  the  driving-point  admittance  Y(s) 

~  sc 

and  the  short-circuit  current  I  (s).  To  define  the  Thevenin  circuit 
one  needs  to  know  the  driving-point  Impedance  Z(s)  and  the  open-circuit 

"''OC 

voltage  V  (s).  These  quantities  are  not  independent  and  are  related 
by 

Y(s)  -  -TT  (3.D 

Z(s) 

and 

V°C (s) 

Z(s)  =  -z — ^  .  (3.2) 

I  (s) 


^The  issues  of  physical  realizability  of  the  equivalent  circuits  are 
addressed  in  Chapter  4. 

2) 

A  port  is  here  understood  in  the  usual  circuit  sense  as  a  pair  of 
terminals  with  the  property  that  the  current  entering  one  terminal 
is  equal  to  the  current  leaving  the  other  terminal. 
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Figure  3.1  (a)  Single-port  electromagnetic  energy  collecting  structure  and  (b)  its  Norton  and 

(c)  Th^venin  equivalent  circuit  representations. 


The  process  of  constructing  the  Norton  (Thevenln)  equivalent  network 
can  thus  be  naturally  divided  into  two  steps:  the  computation  of 
Y(s)(Z(s))  —  the  immittance  problem,  and  the  computation  of 
Isc(s)(Voc(s))  —  the  source  problem. 

In  order  that  the  circuit  quantities  voltage  and  current,  could 
be  uniquely  defined,  the  field  in  the  gap  region  must  be  conservative. 
This  means  that  the  gap  must  be  electrically  small  in  some  sense. 

In  most  cases  this  requirement  restricts  the  class  of  tractable 
bjects  to  slender  structures  with  narrow  gaps  (in  terms  of  wavelengths 
or,  equivalently,  the  rate  of  variation  of  the  transient  waveshapes 
of  interest). ^ 

A  typical  gap  region  of  width  A  and  perimeter  P  is  shown  in 

Figure  3.2.  The  gap  is  defined  such  that  there  are  two  sides  or 

terminals  separated  by  the  gap  region  with  surface  S  .  The  antenna  or 

8 

scatterer  surface  is  designated  by  S  so  that  S=S  US  .  a  is  a  unit 

7  a  a  g  g 

vector  oriented  from  one  gap  side  to  the  other  on  S^.  Following  Beam 
[197va]  we  define  the  gap  current  and  voltage  as  the  averaged  quantities 

I(s)  -  i  <J(r,s) ;ag(r)>s  (3.3) 

and 

V(s)  m  ~  J  <E(r,s) ;ag(r)>s  .  (3.4) 

Since  a  slow  spatial  variation  of  the  field  in  the  gap  is  assumed,  the 
expressions  (3.3)  and  (3.4)  can  be  simplified  to,  respectively, 

^There  are  exceptions  to  this  rule,  e.g .,  the  symmetrically  excited 
spherical  antenna  or  the  thick  cylindrical  antenna.  One  can  even  allow 
non-symmetric  excitations  (e.g.,  plane  wave)  if  the  results  are 
interpreted  properly. 


(3.5) 


I(s)  s  P  J(rg,s)  •  ag(rg) 


V(s)  a  -  A  E(r  ,s)  •  a  (r  ) 

O  O  o 


(3.6) 


with  r  denoting  a  point  in  the  gap  region. 

8 

The  Norton  and  Thevenin  equivalent  circuits  resulting,  respectively, 
from  the  short-circuit  and  the  open-circuit  boundary  value  problems 
are  derived  in  the  next  sections.  These  circuits  are  based  on  the 
dominant  pole  expansion  (2.48)  for  the  current  density  on  the  object. 

In  his  original  work  Baum  [1976a]  also  derived  the  eigenmode  expansion 
networks  obtained  from  (2.13).  These  networks,  however,  are  of  little 
practical  value,  since  in  most  cases  the  quantities  Jn(r,s)  and  \n(s) 
are  not  explicitly  known. 


3.2  Norton  Equivalent  Circuit 

The  Norton  equivalent  circuit  consists  of  the  driving-point 

•v  ^sc 

admittance  Y(s)  and  the  short-circuit  current  I  (s)  (Figure  3.1(a)). 
The  driving-point  admittance  of  the  structure  is  determined  by  solving 
the  short-circuit  boundary  value  problem  Indicated  in  Figure  3.3(a). 
The  incident  field  is  that  due  to  the  voltage  V(s)*Vof(s)  impressed 
across  the  gap,  so  that  the  coupling  coefficients  (2.49)  take  the  form 


n  ■  <J  (r);  -T-  &  (r)>„  -  V  I 

n  n  A  g  Sg  on 


(3.7) 


where  the  notation 


\  5  I 


S  P  Jr,(0  ' 

n  8  g  g 


(3.8) 


t 


has  been  introduced.  Using  the  coupling  coefficients  (3.7)  in  (2.48) 
the  port  current  can  be  determined  from  (3.3)  as 

iu)  -±e„  x»  «(.)  +  l  — ,  sn  ij  »<.)  .  o.9) 


The  driving-point  admittance  Y(s)  *  I (s) /V(s)  follows  directly  from 


(3.9): 


Y(s)  =  -7-  +  l  Y(s) 


sL  L  n 

o  n 


(3.10) 


where  the  Yn(s)  are  defined  by 

a  a  as 

rj  /  v  n  .  n  n 

Y  (s)  *  -  +  -  =  - - r-  . 

n  s-s  s  s(s-s) 

n  n  n  n 

These  are  the  so-called  modified  pole  admittances,  as  distinguished 


(3.11) 


from  the  unmodified  pole  admittances 


Y  (S)  *  - — 

n  s-s 


(3.12) 


which  result  if  the  expansion  (2.25)  is  used.  The  admittance  residues 


a^  are  given  by 


a  -  6 
n  n  n 


and  the  magnetostatic  inductance  Lq  is  defined  as 


L  =  —  . 

o  a 

o 

One  can  combine  conjugate  pole  terms  in  (3.10)  to  get^ 


(3.13) 


(3.14) 


'For  simplicity,  only  complex  conjugate  poles  are  assumed.  A  trivial 
modification  of  (3.15)  is  necessary  to  avoid  the  double  accounting 
of  purely  real  poles. 
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5<s)  "iT  +1  [V9>  +  5-„<s>] 

on1—  — * 


■  ir +  £  Vs)  «-15> 

o  n 

with  the  conjugate  pole-pair  admittances  defined  as 


Y  (s)  =  Y  (s)  +  Y  (s)  .  (3.16) 

n  n  -n 


In  (3.15)  -  (3.16)  the  notation  is  adopted  that  n  and  -n  constitute 

the  indices  for  a  conjugate  pair  of  poles  (s  »s*).^ 

The  short-circuit  current  can  be  found  from  the  boundary  value 

problem  shown  in  Figure  3.3(b).  The  incident  field  is  now  that  due 

to  a  plane  wave  and  the  coupling  coefficients  n  are  given  by 

n 

.  .  -Y_P  •  (r  -  r  ) 

n  -  <J  (r);E  e  8  >«,  . 

n  n  o  S 


(3.17) 


Note  that  the  class  1  coupling  coefficient  form  is  used  with  the  time 
origin  at  the  gap  region,  as  suggested  in  Section  2.5.  The  short- 
circuit  current  can  be  obtained  from  (2.50)  and  (3.3)  as 


(3.21) 


T 

n  I 

n 


have  been  introduced.  Combining  conjugate  pole  terms  in  (3.19)  gives 

ISC(s)  -  gQf(s)  +  l  ^n(s)Tn  +  Y_n(s)T*jf(s)  .  (3.22) 


It  is  fruitful  to  express  Isc(s)  in  terms  of  the  pole-pair  admittances 

/v 

Y  (s)  as  follows: 
n 


Isc(s)  =  g  f(s)  +  l  Yn(s)Tn(s)f(s) 

=  g  f(s)  +  l  Y  (s)V  (s)  (3.23) 

o  “  n  n 

n 

/N 

where  the  pole-pair  voltage  transfer  functions  T  (s)  are  defined  as 

n 


Y  (s)T  +  Y  (s)T 


Y_(s)  -  Y_n(s) 


f„‘s)  •  '  1  (sVt-?— (s-)~°  ■  SetTn>  '"“n1  — ' T~; 

n  -n  Y  (s) 


(3.24) 


and  the  pole-pair  source  voltages  vn(s)  as 


V  (s)  =•  T  (s)f(s) 
n  n 


(3.25) 


The  driving-point  admittance  (3.15)  and  the  short-circuit  current 
(3.23)  can  be  combined  according  to  the  Norton's  theorem  to  give  the 
equivalent  network  shown  in  Figure  3.4(a).  By  pairing  pole-pair 
admittances  with  the  corresponding  pole-pair  source  voltages  and  applying 
the  Thevenin's  transformation  to  each  pole-pair  circuit  this  network 
can  be  recast  into  the  form  shown  in  Figure  3.4(b).  This  very  fruitful 
transformation  was  first  introduced  by  Baum  [1976a]. 
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3.3  Thevenin  Equivalent  Circuit 

The  Thevenin  equivalent  circuit  results  from  the  open-circuit 
boundary  value  problem  shown  in  Figure  3.5.  The  solution  for  the 
open-circuit  SEM  parameters  is  in  practice  a  much  more  difficult 
task  than  the  solution  for  the  short-circuit  parameters  of  the  same 
structure.  Thus,  it  is  probably  more  practical  to  derive  the  Thevenin 
circuit  indirectly,  from  the  short-circuit  parameters.  This  last 
approach,  which  was  pursued  recently  by  Tesche  and  Gird  [1981],  is 
briefly  discussed  at  the  end  of  this  section. 

The  forcing  function  in  the  impedance  problem  (Figure  3.5(a)) 
is  the  impressed  gap  current 


Jg(£,s)  -  Jg(r)f (s)  s  ±  ag(rg)Iof(s) 

-  i  ag(rg)I(s)  .  r  e Sg  .  (3.26) 

The  field  due  to  this  current  can  be  otained  from  the  integral 
equation  as 


E  (r,s)  =  -  <Z(r,r',s);J  (r,s)> 


-Io  p  <Z(r,r's,);ag(r')>s  f(s) 


(3.27) 


and  the  coupling  coefficients  are  given  by 


nn(s)  -  -I0  j  <Jn(r);<Z(r,r\s);a  (r')>  >g  =  -I  a  (s)  (3.28) 

8  g  a 

where,  for  notational  convenience,  a  new  quantity  has  been  defined: 


an(s)  ■  j  <Jn(r);<Z*’(r,r'  ,s);a  (r’)>  > 

8  g  sa 


(3.29) 
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The  current  on  the  object  can  be  now  expressed  with  the  help  of  (2.48) 

1) 


as 


J(r,s)  *  -  7  8  5  (s)J  (r)I(s)  -  £  ■  8  a  (s)J  (r)I(s) .  (3.30) 

s  o  o  o  £  sn(s-sn)  n  n  n 

The  electric  field  radiated  by  this  current  and  the  gap  current  is 
given  by 


ES(r,s)  =  -  <Z(r,r',s);J(r',s)>s 
-  <Z(r,r',s);Jc(r',s)> 


(3.31) 


g 

Using  the  last  expression  in  (3.4)  gives  the  voltage  across  the  gap  as 

V(s)  *  i  <<z(r,r ' ,s) ;J(r ’ ,s)>s  ;a  (r)> 

a  8  g 

+  4  <<?(r,r's);J (r,s)>  ;a  (r’)>  .  (3.32) 

P  8  g  8  g 

Substituting  (3.26)  and  (3.30)  into  (3.32)  and  using  the  symmetry 

properties  of  the  impedance  kernel  yields 


V(s)  »  a  (s)  -  7  8  a^(s)  -  l 
g  s  o  o  u 

n 

where  a  new  quantity  has  been  defined: 


s 

n 


(s  ”  Sn) 


8na^(s)I(s) 


(3.33) 


ag(s)  *  pr  <ag(r)'Z(r,r,,s);Sg(r,)>s  .  (3.34) 

*► 

Since  J  (r)  is  solenoidal,  it  can  be  seen  from  (3.29)  and  (2.3) 
o 

that  3  (s)  ■*  s  as  s-*-0,  so  that  the  second  term  of  (3.33)  goes  to  zero 
o 

at  this  frequency.  For  n  >  0,  however,  3n(s)  introduces  a  pole  at  zero. 
Also,  3g(s)  has  a  pole  at  the  origin.  Expression  (3.33)  can  be  reexpanded 


It  should  be  kept  in  mind  that  the  SEM  parameters  sn,  ?n(r),  and 
Sn  are  computed  with  the  gap  open  ard  are,  in  general,  different  from 
the  short-circuit  quantities  used  in  the  previous  section. 
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to  give  the  driving-point  impedance  Z(s)  ■  V(s)/I(s)  in  the  form 

■’21 

-  I  "i  ~~  i~  ~~  B_af 


Z(s) 


6  a' 

a'  +  l  -V- 

8  n  s  n 


/  \  M  vA 

s  (s  -  s  )  n  n 
n  n  n 


(3.35) 


where 


[s  a  (s)l 

L  (5)  J 


a'  =  £im  _ 

(8)  L  ® 

and  a  =a  (s  ).  It  is  convenient  to  define  new  quantities,  the 
n  n  n 

electrostatic  capacitance 


(3.36) 


“s  +  2  l  Re{-7^-} 


(3.37) 


and  the  impedance  residues 


B  or 

n  n 


(3.38) 


and  rewrite  (3.35)  as 


Z<»)  *  TF"  +  l  8JS> 


sC  *•  n 

o  n 


(3.39) 


with  the  modified  pole  impedances  Z^Cs)  given  by 

-  a '  s 

z  ( s )  -  ~  77  — T  .  (3.40) 

n  s  (8-s  ) 

n  n 

The  impedance  counterpart  of  (3.12),  the  unmodified  pole  impedances, 
are  defined  as 


Z  (s)  -  - — 

n  s  -  s 


(3.41) 


Combining  the  conjugate  pole  terms  in  (3.39)  gives 
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(3.42) 


zu)  -  jf  *  l  [~VS>  +  l. 

o  n  u 


„<s>] 


*  i  vs> 

o  n 


with  the  conjugate  pole-pair  impedances  defined  as 


Z  (s)  -  Z  (s)  +  Z  (s)  . 
n  n  -n 


(3.43) 


Note  that  for  multiply-connected  objects,  such  as  the  loop  antenna, 

(3.42)  does  not  exhibit  the  proper  behavior  as  s-*-0,  unless  C  «*<». 

o 

This  means  that  the  denominator  of  (3.37)  must  sum  to  zero  in  such  case. 
Also,  note  that  all  poles  contribute  to  the  static  capacitance  Co>  as 
given  by  (3.37).  Numerical  evaluation  of  this  expression  shows  that  if 
only  dominant  poles  are  included,  (3.37)  significantly  underestimates 
the  value  of  Cq.  Much  better  results  can  be  obtained  by  computing 
Cq  apart  from  the  SEM  formulation,  by  solving  for  the  static  charge 
distribution  on  the  structure. 

The  open-circuit  voltage  can  be  found  from  the  boundary  value 
problem  shown  in  Figure  3.5(b).  The  incident  field  is  that  due  to  a 
plane  wave  and  the  coupling  coefficients  are  given  by 


+  .  .  ~Y  P  *  (r-r  ) 

H  -  <5  (r);E  e  8  >e 

n  n  o  c 


(3.44) 


The  current  on  the  object  is  given  by  (2.50)  as 


J(r,s)  »  uoBoJo(r)f  (s)  +  \  -  J  n_ j_ ( r ) f  ( s ) 


/  v  MM*-'  ^ 

s(s-s)  nnn 
n  n  n 


(3.45) 


and  the  open-circuit  voltage  is  given  by  (3.4)  as 


V0v"(s)  ■  i  <<Z*(r,r’,s);J(r',s)>s  ;a  (r)>  . 

a  8  g 

Substituting  (3.45)  into  (3.46)  gives 


(3.46) 
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where  the  pole  pair  transfer  function  T^(s)  are  defined  by 


r(s) 


Zn(s>T;  +  2-n<s>T;* 

Zn<s>+Z4s> 

Z  (s)  -Z  (a) 

Re(T'}  +j  I  {T'}  — — - — - 

n  m  Z  (s) 


(3.55) 


and  the  pole-pair  source  currents  by 


I  (s)  =  T  (s)f(s)  . 
n  n 


(3.56) 


The  driving-point  impedance  and  the  open-circuit  voltage  can  be  now 
combined  according  to  the  Thevenin's  theorem  to  give  the  circuit 
shown  in  Figure  3.6(a).  By  pairing  pole-pair  impedances  with  the 
corresponding  pole-pair  source  currents  and  applying  the  Norton's 
transformation  to  each  pole-pair  circuit  this  network  can  be  recast 
into  the  form  shown  in  Figure  3.6(b). 

As  mentioned  at  the  beginning  of  this  section,  the  open-circuit 
SEM  parameters  are  much  more  difficult  to  obtain  than  the  short-circuit 
parameters  of  the  same  problem.  This  is  so  because  the  highly  localized 
charge  interaction  effects  between  the  two  parts  of  the  structure  in 
the  vicinity  of  the  gap  must  be  modeled  precisely  in  the  open-circuit 
boundary  value  problem.  Furthermore,  it  is  enough  to  solve  only  one 
short-circuit  boundary  value  problem  for  all  possible  gap  locations  on 
the  structure,  whereas  the  open-circuit  problem  must  be  resolved  each 
time  the  gap  is  moved.  Thus,  if  for  any  reason  the  Thevenin  equivalent 
circuit  is  preferred  over  the  Norton  equivalent,  it  is  likely  to  be 
expedient  to  derive  it  from  the  short-circuit  parameters  using  the 
relations  (3.1)  and  (3.2).  Using  (3.15)  in  (3.1)  gives  the  driving- 
point  impedance  as 
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Z(s) 


Y(s) 


1  c  • 

- + 1  Vs) 


(3.57) 


sL 


o  n 

This  expression  can  be  easily  recast  in  a  form  of  a  ratio  of  two 
polynomials  in  s,  and  the  open-circuit  poles  can  be  found  by  solving 
for  the  roots  of  the  denominator  polynomial.  Expanding  Z(s)  in  terms 
of  these  poles  gives  the  input  impedance  in  the  desired  form  (3.39). 

It  should  be  noted  that  the  inductive  term  in  (3.57)  is  not  present 
in  the  case  of  a  simply-connected  structure,  such  as  the  dipole  antenna. 
Such  a  structure  exhibits  a  capacitive  behavior  at  low  frequencies, 
with  the  static  capacitance  given  by^ 


■2 1 Hi) 


(3.58) 


where  a  are  the  admittance  residues  and  s  are  the  short-circuit 
n  n 

poles. 

In  like  fashion,  the  open-circuit  voltage  can  be  otained  by 
substituting  (3.15)  and  (3.23)  into  (3.2)  as 


tsc,  . 

V°c(s)  «  K-&- 
Y(s) 


gof(s)  +  l  Yn(s)Vn(s) 
n 


ir +  £  f„(s) 

o  n 


(3.59) 


In  the  case  of  dipole-like  structures  (simply-connected  objects)  the 
first  terms  in  the  numerator  and  the  denominator  of  (3.59)  are  not 
present,  and  this  expression  can  be  easily  recast  into  the  form  of 


(3.54)  with  the  coefficient  k^  given  by 


^Complex  conjugate  poles  are  assumed.  The  extension  to  purely  real 
poles  is  trivial. 
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where  C  is  given  by  (3.58)  and  T  are  the  source  coefficients 
o  n 


defined  by  (3.21). 
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CHAPTER  4 


SYNTHESIS  OF  THE  POLE-PAIR  CIRCUITS 

4.1  Introduction 

The  formal  equivalent  networks  developed  in  Chapter  3  have  the 

form  of  a  connection  of  simple  modules  associated  with  conjugate  pole 

pairs.  This  network  topology  simplifies  considerably  the  synthesis 

process.  Furthermore,  it  has  the  advantage  that  it  can  be  easily 

expanded  by  adding  more  pole-pair  modules  without  repeating  the  whole 

synthesis  process.  Since  the  pole-pair  modules  are  simple,  one  can 

hope  to  give  simple  prescriptions  for  their  physical  realization.  To 

define  such  prescriptions  is  the  main  objective  of  this  chapter. 

As  the  first  step  toward  this  goal,  the  Norton  and  Thevenin 

networks  from  Figures  3.4(b)  and  3.6(b)  are  recast  into  the  form  given 

in  Figures  4.1  and  4.2,  respectively.  These  networks  are  in  a  physically 

realizable  form  provided  the  admittances  Y  (s)  and  transfer  functions 

n 

A  A 

T  (s)  in  the  Norton  network,  and  the  impedances  Z  (s)  and  transfer 
n  n 

A 

functions  T^(s)  in  the  Thevenin  network,  can  be  physically  realized. 

It  should  be  noted  that  these  networks  contain  controlled  sources. 

In  a  hardware  realization  they  can  be  implemented  by  active  devices, 
such  as  operational  amplifiers. 

A 

The  next  sections  are  devoted  to  the  synthesis  of  the  Y  (s)  and 

n 

A  /\ 

Zn(s)  (driving-point  function  synthesis),  and  the  synthesis  of  Tr(s) 

a 

and  T'(s)  (voltage  transfer  function  synthesis).  Both  modified  and 
n 

unmodified  pole  modules  are  considered. 

Some  of  the  material  of  this  chapter  is  not  new  and  is  included 
here,  for  the  sake  of  continuity,  in  the  notation  consistent  with  the 
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Thevenin-based  equivalent  ci 


rest  of  this  work.  This  is  true  in  particular  of  the  realizability 
conditions  of  the  unmodified  pole-pair  modules  and  their  interpretation. 
Analogous  conditions  which  we  give  for  the  modified  pole-pair  immittances 
are  believed  to  be  new,  however.  Similarly,  the  padding  procedure 
addressed  in  Section  4.3  and  the  approximate  Bott-Duffin  modules  were 
used  previously,  but  the  explicit  expressions  for  the  necessary  padding 
and  network  element  values  in  terms  of  the  SEM  parameters  appear  here 
for  the  first  time.  In  deriving  them  we  benefited  from  an  early  work 
by  Foster  [1955].  Most  of  the  material  on  the  pole-pair  transfer 
function  synthesis  in  Section  4.4  is  also  original. 


4.2  Driving-Point  Function  Synthesis 

The  conjugate  pole-pair  admittance  (2.16)  can  be  written  explicitly 
in  a  form  of  a  biquadratic  function  of  s  as  follows: 

*  A  S2  +B  s  +  C 

y  (S)  *  -2 - 2 - 2.  (4.1) 

n  D  s2  +  E  s  +  F 

n  n  n 

with  the  coefficients  defined  differently  depending  on  whether  the  modified 
or  unmodified  pole  admittances  are  used.  For  the  unmodified  modules 
the  numerator  coefficients  of  (4.1)  are  given  by 


A  ■  0  , 

n  * 

2c 


and 


c  -  d  q 
n  n*  n 


and  the  denominator  coefficients  by 


(4.2a) 
(A. 2b) 

(4.2c) 
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-1 


(4.3  a) 


D 

n 

E 

n 


=  Q. 


F  = 
n 


(4.3  b) 
(4.3  c) 


For  the  modified  modules  the  numerator  coefficients  of  (4.1)  are  given  bv 


A 

n 


B 

n 


d  q  -  c 
n  n  n 

Q  Is  |2 

n 1  n 


c  (2Q2  -  1)  +  d  q 
n  _ n^n 

Q2|s  | 

n 1  n  1 


(4.4a) 


(4.4b) 


and 


(4.4c) 


with  the  denominator  coefficients  still  given  by  (4.3).  In  (4.2)-(4.4) 
the  following  notation  is  introduced: 


s  = 

n 

-°n  +  jWn  * 

(4.5) 

a  * 
n 

Cn+Jdn  ’ 

(4.6) 

|s| 

«n  " 

n 

(4.7) 

2a  ’ 

n 


and 

U)  _ 

■  f  ■  ■X  - 1  •  <4-« 

n 

where  the  residues  a  are  defined  by  (3.13)  and  where  Q  are  the 

n  J  xn 

quality  factors  of  the  poles  s^. 

The  conjugate  pole-pair  impedances  Z  (s),  defined  in  (3.43),  can 

n 

also  be  written  in  the  form  (4.1)  with  the  only  difference  that  the 
residues  a^  are  given  by  (3.38)  and  the  poles  are,  of  course,  the 
open-circuit  poles.  In  what  follows  we  deal  specifically  with  the 
driving-point  admittance  function  (4.1).  However,  all  conclusions 
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hold  also  for  the  impedance  function  Z  (s).  The  circuits  which  realize 

n 

A  A 

Z  (s)  are  just  the  duals  of  the  circuits  for  Y  (s). 
n  J  n 

The  necessary  (but  not  sufficient)  condition  for  the  physical 

A 

realizability  of  Y^Cs)  is  that  the  denominator  coefficients  of  (4.1) 
be  positive  and  that  the  numerator  coefficients  be  non-negative. 

This  requirement  leads  to  the  following  restriction  on  the  residue 
location  if  the  unmodified  pole  modules  are  used: 


and 


c  >  0 
n  — 


c  -  d  q  >  0 
n  n^n  — 


(4.9a) 


(4.9b) 


If  the  modified  pole  modules  are  employed,  the  conditions  are 


and 


d  q  -  c  >  0 
n^n  n  — 


c n(2Q„-D  +  d  q  >0 
n  n  n  n  — 


(4.10a) 


(4.10b) 


The  following  discussion  can  be  simplified  considerably  by  assuming 


Q  >  —  s  0.707  . 

n  fi 


(4.11) 


This  assumption  does  not  seem  to  be  restrictive  since  the  dominant  poles 
of  all  structures  considered  so  far  satisfy  (4.11),  Even  the  dominant 
poles  of  a  sphere,  which  is  an  extremely  low-Q  structure,  satisfy  Qn^l* 
Assuming  that  (4.11)  holds,  (4.10)  can  be  rewritten  as 

9. 


n  _  n  „ 

-  c  —  <  q 

2  ,  -  d  —  nn 


and 


2Q  -  1  n 
n 


d  >  0 
n  — 


(4.12a) 


(4.12b) 
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The  driving-point  admittance  function  Y  (s)  is  physically  realizable 

n 

if  it  is  positive  real  (PR)  [Bmcne„  1931],  i.e.a  if  it  satisfies  the 
following  conditions: 

1.  Re{Yn(s)}>0  for  Re{s}  >  0;  and 

2.  Y  (s)  is  real  when  s  is  real. 

n 

Since  it  is  extremely  difficult  to  apply  these  conditions  directly, 
the  following  equivalent  requirements,  which  are  easier  to  check,  can 
be  employed: 

A.  Yn(s)  cannot  have  any  poles  or  zeros  in  the  right  half  of 
the  complex  frequency  plane; 

B.  any  ja)-axis  poles  of  Yn(s)  must  be  simple  and  have  positive 
real  residues;  and 

C.  Re{Yn(jai) }  0  for  all  <a. 

In  the  case  of  the  biquadratic  function  (4.1)  the  conditions  A  and 
B  above  are  automatically  satisfied  if  (4.9)  and  (4.12)  hold,  respectively, 
for  the  case  of  unmodified  and  modified  pole  modules.  It  only  remains 
to  check  the  condition  C. 

A 

The  behavior  of  Re{Yn(ju)}  can  be  investigated  with  the  help  of 
the  function 


P„(s)  -■![?„(,) +?„(-s)]  (4.13) 

which  on  the  jto  axis  is  equal  to  the  real  part  of  Y^(s)  there 
t Guillemin ,  1977].  For  the  unmodified  pole  modules  case  Fn(s)  has  the 
form 


Vs> 


(CD  -  B  £  )sJ 
n  n  n  n 


+  C  F 
n  n 


D2sl*  +  (2D  F  -  E2)s2  +  F2 
n  n  n  n  n 


(4.14) 


In  addition  to  a  double  zero  at  infinity  this  function  has  two  zeros 


given  by 
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/  c"f 

c  a  +  /  - 

o  BE  -CD 

n  n  n  n 

If  the  zeros  are  real,  i.e.,  if 

B  E  -C  D  >0 
n  n  n  n  — 


(4.15) 


(4.16) 


the  pole-zero  pattern  of  F  (s)  is  that  shown  in  Figure  4.3(a)  and 

n 

the  corresponding  behavior  of  Re{Y  (ju)}  is  that  shown  in  Figure 

n 

(4. 3(b). ^  If  the  condition  (4.16)  is  violated,  the  zeros  (4.15) 
become  purely  imaginary,  as  shown  in  Figure  4.3(c),  and  the  character 

A 

of  the  real  part  of  Y^tju))  changes  to  that  shown  in  Figure  4.3(d). 

It  is  evident  from  this  figure  that  in  this  case  there  is  a  negative 

A 

excursion  in  the  real  part  of  Y  (s)  on  the  ju-axis,  and  the  condition 

n 

C  is  violated.  Hence,  (4.16)  is  a  necessary  condition  for  the  PR-ness 
of  (4,1)  for  the  case  of  unmodified  pole  modules.  Expressed  in  terms 
of  the  SEM  parameters  this  condition  takes  the  form 


c  +  d  q  >  0  .  (4.17) 

n  n  n  — 

Combining  (4.9)  and  (4.17)  yields 

c 

>  q  (4.18) 

—  n 

which  is,  therefore,  the  necessary  and  sufficient  condition  for  the 
physical  realizability  of  the  unmodified  pole  pair  admittances.  This 
condition  was  stated  previously  by  Guillemin  [1977]. 

For  the  modified  pole  modules,  Fn(s)  takes  the  form 


K! 


If  (4.16)  is  satisfied  with  the  equality  sign, ..the  zeros  (4.12) 
move  to  infinity.  However,  the  character  of  Re{Yn(jco)}  is  still 
that  shown  in  the  Figure  4.3(b). 
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(4.19) 


A  D  sz  +  (A  F  -BE) 
2  n  n  nn  nn 

F  (s)  ■  s  - 

n  -2„4 


Dzs“  +  (2D  F  -  E2)s2  +F2 
n  n  n  n  n 

This  function  has  a  double  zero  at  the  origin  and  a  pair  of  zeros 
given  by 


/BE  -A  F 
+  inn  n  n 

~  /  AD 
n  n 


(4.20) 


If  the  condition 


B  E  -  A  F  >  0  (4.21) 

n  n  n  n  — 

is  satisfied,  the  pole-zero  pattern  of  F^(s)  and  the  character  of 

A 

the  real  part  of  Yn(ju>)  are  that  shown  in  Figure  4.4(a)  and  (4.4(b), 
respectively.  If  the  condition  (4.21)  is  violated,  the  pole-zero 
pattern  of  FR(S)  and  the  character  of  the  real  part  become  that  shown 
In  Figure  4.4(c)  and  (4.4(d),  respectively.  Again  in  this  case,  there 
is  a  negative  excursion  of  the  real  part  on  the  ju>-axis.  Hence,  (4.21) 
is  a  necessary  condition  for  the  realizability  of  (4.1)  in  the  modified 
pole-pair  admittances  case.  This  condition  may  be  written  in  terms 
of  the  poles  and  residues  as 


c  (3Q2  -1)  -  d  (Q2  -  1)  q  >  0 
n  n  n  n  n  — 


(4.22) 


or,  through  the  use  of  (4,11)  and  (4.12a)  as 


c  Q  - 1 

— 2 -  q 

dn  3Q2  -  1  n 
n 


(4.23) 


Combining  (4.12)  and  (4.23)  gives  the  necessary  and  sufficient  condition 
for  the  physical  realizability  of  the  modified  pole-pair  admittances  as 


Q2  -1  c 

_2 _  q  <  <  q 

3Q2  -  1  n  -  dn  "  n 


(4.24a) 
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function  of  the  modified  pole  module  and  the 


with 


(4.24b) 


d  >  0  . 

n  — 

If  Qn  » 1,  which  is  satisfied  for  most  dominant  poles  of  highly- 
resonant  structures,  the  first  condition  above  is  given  approximately  by 

?  c 

tQ  </<2Q  (4.25) 

j  n  —  d  —  n 
n 

The  conditions  for  realizability  of  the  unmodified  and  modified 

pole-pair  modules  can  be  interpreted  geometrically  in  the  a  -plane  as 

n 

shown  in  Figure  4.5.  If  the  residue  falls  into  region  I,  the 
unmodified  pole-pair  admittance  is  realizable;  if  it  falls  in  region 
II,  the  modified  pole-pair  admittance  is  realizable;  if  the  residue 
lies  outside  the  regions  I  and  II,  neither  form  is  realizable. 

The  realizability  conditions  for  the  modified  pole-pair  admittances 
were  given  previously  by  Streable  and  Pearson  [1981]  in  a  form  not 
amenable  to  any  geometrical  interpretation.  Similar  conditions  derived 
by  Beam  and  Singaraju  [1980]  are  necessary  but  not  sufficient. 

If  the  PR-ness  conditions  are  satisfied,  the  unmodified  and  modified 
pole-pair  admittances  can  be  synthesized  by  a  continued-fraction 
expansion  [Matthaei ,  1954;  Miehalski  and  Pearson ,  1979]  in  a  form 
shown  in  Figures  4.6  and  4.7,  respectively.  (The  index  n  on  the 
coefficients  of  (4.1)  is  dropped  in  the  circuits  given  in  Figures  4.6 
and  4.7,  and  other  figures  in  this  chapter.)  It  should  be  noted  that 
two  equivalent  realizations  are  possible  in  each  case. 

4.3  Synthesis  of  Approximate  Driving-Point  Circuits 

It  has  been  found  that  for  most  poles  of  the  thin-wire  dipole  and 

loop  antennas,  as  well  as  the  spherical  antenna,  the  residue  a  falls 

n  , 

t.1.^  .  ...  T  ..J  TT  _  C  ri  AC  J  —  ^1 _ _ 1  _ _ J  _ 


1*1*1 


I-UNMODIFIED  ADMITTANCE  REALIZABLE 
IT-MODIFIED  ADMITTANCE  REALIZABLE 


Figure  4.5  Geometrical  interpretation  of  the  realizability  conditions 
for  the  modified  and  unmodified  pole-pair  admittances. 
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*.6  Practical  realizations  of  the  pole-pair  admittance  function  when  the  residue  lies  in 
[  of  Figure  4.5.  (Q = BE-CD) 


admittance  is,  in  most  cases,  not  physically  realizable.  This  manifests, 
itself  as  a  negative  excursion  in  the  real  part  of  the  driving-point 
function  on  the  ju-axis,  as  discussed  in  Section  4.2.  It  has  been 
observed,  however,  that  in  the  case  of  modified  pole  modules  this 
negative  excursion  is  in  most  cases  negligible,  or  at  least  small,  in 
comparison  with  the  value  of  the  real  part  near  the  resonance.  Thus, 
if  a  small  conductance  G  is  added  Y  (s),  the  resulting  function  cun  be 
rendered  PR  without  introducing  appreciable  error.  This  procedure  is 
called  a  real-part  padding  by  Guillemin  [1977]  and  was  first  applied 
in  the  present  context  by  Streable  and  Pearson  [1981]  in  synthesizing 
driving-point  equivalent  circuits  for  the  thin-wire  dipole  and  loop 
antennas.  It  should  be  noted  that  neither  the  poles  nor  the  residues 
are  affected  in  the  padding  process.  Only  the  zeros  of  the  pole-pair 
module  are  changed. 

The  amount  of  padding  which  must  be  added  to  a  pole-pair 
admittance,  so  that  the  condition  C  of  Section  4.2  is  satisfied, 
can  be  found  explicitly  in  terms  of  the  poles  and  residues  to  be 

Q  (2Q  |  a  )  +  d  -cq) 

_  <n  ^n1  n'  n  n'n  ,,  ,n 

Gn - ;  ;  (4-26) 

qJv 

for  the  unmodified  pole-pair  module,  and 

2Q3 1 a  |  -  d  (3Q2  -  1)  -  c  (Q2  -  1 ) q 
Gn  -  —  n-  q - - - - - 2—2 - -  (4.27) 

Qn|s>n 


(4.27) 


for  the  modified  module.  The  numerator  coefficients  of  the  padded 
pole-pair  admittance  (4.1)  can  be  expressed  in  terms  of  the  poles  and 
residues  as  follows: 


Q  ( 2Q  |  a  |  +  d  -  c  q  ) 
n  n  1  i)  _ n _ n  n 

’  „  I  ..  I  z 

n  i  !■> 

n  n 


<4.28.-i) 


i 


2Q  la  +  d  +  c  q 
n  n1  n  n^n 


q  s 
^n '  n  1 


(4.28b) 


and 


2Q3  a  -c  (Q2  -  l)q  -d  (3Q2  -  1) 
n'  n1  n  n  n  n 

q  0 

n 


(4.28c) 


The  denominator  coefficients  are  not  affected  and  are  given  by  (4.3). 
It  should  be  noted  that  after  padding  all  coefficients  of  the 
biquadratic  function  (4.1)  are  nonzero,  which  makes  the  synthesis  of 
it  more  difficult. 

It  can  be  shown  [Foster,  1955]  that  the  necessary  and  sufficient 
condition  for  the  physical  realizability  of  (4.1)  is 


/a"  F  -  Jc  D 


n  n 


n  n 


<  B  E 
—  n  n 


(4.29) 


It  is  easy  to  see  that  this  condition  is  consistent  with  (4.16)  in 
the  case  when  A^  *  0,  and  with  (4.21)  when  0^  =  0.  The  padded  pole-pair 
admittance  is  a  minimum-real-part  function,  i.e.,  the  real  part  is  zero  at 
some  finite  frequency  [ Balabanian ,  1958],  For  such  functions  (4.29) 
is  satisfied  with  the  equality  sign,  i.e.. 


A  F  -  /C  D  -  ±  /b  E  .  (4.30) 

n  n  n  n  n  n 

It  appears  [Miohalski  and.  Pearson,  1979]  that  the  only  general  synthesis 
methods  which  can  be  applied  in  this  case  are  the  Brune  method  [ Bruno , 
1931]  and  the  Bott-Duffin  method  [Bott  and  Duffin,  1949].  The  synthesis 
of  (4.1)  by  these  methods  is  described  in  considerable  detail  in 
[Miohalski  and  Pearson,  1979].  Here,  we  quote  only  the  principal 
results. 
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The  Brune  network  realizing  (4.1)  is  shown  in  Figure  4.8  with 
explicit  expressions  for  the  element  values  in  terms  of  the  coefficients 
of  (4.1).  The  Brune  network  is  simple  and  can  be  successfully  used 
for  analysis  purposes.  It  is  less  attractive  for  a  physical  realization, 
because  of  the  presence  of  a  unity-coupled  transformer.  Even  a  slight 
deviation  of  the  coupling  coefficient  from  unity  gives  rise  to  an 

A 

additional  zero  of  Yn(s)  at  infinity.  This  does  not  seem  to  have  a 
very  drastic  effect  on  the  behavior  of  the  function,  however. 

An  application  of  the  Bott-Duffin  procedure  to  (4.1)  results 
in  two  different  network  topologies  depending  on  the  sign  on  the  right 
hand  side  of  (4.30).  For  the  "+"  sign  the  network  shown  in  Figure 
4.9(a)  results,  whereas  for  the  sign  the  network  in  Figure  4.9(b) 
is  applicable.  These  two  situations  will  be  referred  to  as  cases  A 
and  B,  respectively.  An  explicit  network  for  the  case  B  (or,  rather, 
the  dual  of  it)  was  given  previously  by  Foster  [1955].  Note  that  the 
Bott-Duffin  networks  are  in  a  form  of  a  balanced  bridge,  i.e.,  the 
points  "A"  and  in  Figures  4.9(a)  and  4.9(b)  are  at  the  same  potential, 
and  any  impedance  (even  a  short-circuit)  can  be  inserted  between  them 
without  affecting  the  terminal  admittance.  This  fact  can  be  used  to 
reduce  by  one  the  number  of  reactive  elements  by  inserting  a  capacitor 
or  an  inductor  of  a  proper  value  into  the  bridge  arm,  and  performing 
a  delta-wye  or  a  wye-delta  transformation.  The  resulting  networks  for 
the  two  cases  above  are  shown  in  Figure  4,10  and  4.11.  Again,  the 
explicit  reduced  networks  for  the  case  B  were  given  previously  by 
Foster  [1955]. 

Even  in  the  simplified  form  the  Bott-Duffin  networks  are  rather 
complicated.  They  are  generally  considered  unattractive  from  the  point 
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Figure  4.9  Bott-Duffin  realizations  of  the  padded  pole-pair 
admittance  function  for  (a)  case  A  and  (b)  case  B. 
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Q=A/CE+B/BD 


(b) 


Figure  4.10  Modified  Bott-Duffin  networks  for  the  padded  pole- 
pair  admittance  function  —  case  A. 
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Ra  BVDF  +  E  /AC 

(b) 


Figure  4.11  Modified  Bott-Duffin  networks  for  the  padded  pole- 
pair  admittance  function  —  case  B. 


of  view  of  sensitivity,  as  element  variations  from  their  nominal 
values  give  rise  to  three  additional  poles  and  zeros  in  the  driving- 
point  function. 

In  order  to  determine  which  of  the  two  cases  A  or  B  is  applicable 
for  a  given  set  of  SEM  parameters,  it  is  helpful  to  consider  the 
expression 


d  q  -  c 

A  F  -  C  D  »  —  - - — 

n  11  n  n  |s  |Q 
n  n 

It  is  easy  to  see  from  (4.30)  that  the  case  A  arises  when  (4.31)  is 
positive;  if  it  is  negative,  the  case  B  is  applicable.  This  can  be 
summarized  as  follows: 


(4.31) 


d  q  -  c 
n  n  n 


>  0 
<  0 


case  A 

case.  B 


(4.32) 


The  geometrical  interpretation  of  this  condition  in  the  a^-plane  is 
given  in  Figure  4.5  with  the  regions  of  applicability  of  the  two  cases 
explicitly  indicated. 

It  can  be  shown  that  the  resistances  in  Figure  4.9(a)  and 
R^  in  Figure  4.9(b)  are  reciprocals  of  Gn  -  the  amount  of  padding  used. 
Thus,  if  the  pole-pair  admittance  is  almost  PR,  so  that  the  padding 
is  negligible,  these  resistances  are  large  and  the  entire  circuit  legs 
in  series  with  them  can  be  neglected.  The  Bott-Duffin  networks  simplified 
in  this  way  are  shown  in  Figure  4.12.  The  approximate  pole-pair  module 
for  the  case  A  (Figure  4.12(a))  was  first  derived  by  Stredble  and 
Peax’son  [1981].  They  did  not  give  explicit  expressions  for  the  element 
values,  however. 
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\/£F  -v5d  =  +i/BE 

(a) 


i/5F  -  Vcd  =  -VIE 

(b) 


Figure  4.12  Simplified  Bctt-Duffin  circuits  approximating  the  padded 
pole-pair  admittance  function  in  the  case  of  negligible  padding. 

(a)  circuit  for  case  A  and  (b)  circuit  for  case  B. 


In  all  cases  considered  to  date  (the  thin-wire  dipole  and  loop 
antennas,  the  spherical  antenna)  only  the  case  A  has  arisen.  It  seems, 
however,  that  this  cannot  be  generalized,  and  situations  where  the 
case  B  is  appropriate  cannot  be  excluded  a  priori. 


4.4  Voltage  Transfer  Function  Synthesis 

The  voltage  transfer  function  of  the  Norton  network  (3.24)  and 
its  Thevenin  network  counterpart  (3.55)  can  both  be  written  in  the 
form 


with  the  coefficients 
of  (4.1))  given  by 


Tn(s) 


A  s  +  B 
n  n 

C  s  +D 


(4.33) 


n  n 

(which  should  not  be  confused  with  the  coefficients 


'  Eela„Vj}  ■ 


and 


Is  I  (c  -  d  q  ) 
1  n1  n  n  n 

2Q_ 


if  the  unmodified  pole-pair  modules  are  used,  and  by 


(4.34a) 

(4.34b) 

(4.34c) 


(4.34d) 


and 
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d  q  +C  (2Q*-1) 
n  n  n  n 


(4.35c) 


(4.35d) 
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if  the  modified  pole-pair  modules  are  employed.  In  (4.34)-(4.35)  a^ 
are  the  immittance  residues  given  either  by  (3.13)  for  the  Norton 
formulation,  or  by  (3.38)  for  the  Thevenin  formulation.  The  source 
coefficients  are  defined  in  these  two  cases  by  (3.21)  and  (3.50), 
respectively. 

Before  proceeding  with  the  synthesis  of  (4.33),  several  comments 
are  in  order.  First,  note  that  the  equivalent  circuits  transformations 
due  to  Baum  [1976a]  indicated  in  Figures  3.4  and  3.6  lead  to  very  simple 
pole-pair  voltage  transfer  functions  which  are  bilinear  in  s.  Without 
this  transformation  one  would  have  to  deal  with  biquadratic  transfer 
functions,  which  would  be  a  more  difficult  task.  Second,  note  that 
each  transfer  network  in  Figures  4.1  and  4.2  is  followed  by  a  controlled 
source.  The  main  purpose  of  these  controlled  sources  is  to  provide 
the  impedance  buffering  which  is  required  because  synthesis  techniques 
do  not  exist  which  realize  a  given  transfer  function  for  a  general, 
unspecified  loading  conditions.  In  addition,  these  controlled  sources 
serve  as  amplifiers.  The  amplifier  function  is  necessary  because, 
in  general,  most  synthesis  techniques  realize  a  given  transfer  function 
only  to  within  a  constant  multiplier,  with  the  expectation  that  the 
gain  can  be  restored  at  a  later  stage. 

The  theory  of  transfer  function  synthesis  is  a  vast  and  mature 
field  [e.g.,  Balabanian,  1958].  Fortunately,  the  function  (4.33) 
is  so  simple  that  only  a  small  part  of  the  accumulated  knowledge  needs 
to  be  applied.  Before  attempting  to  synthesize  (4.33)  it  is  appropriate 
to  quote  some  of  the  conditions  which  a  general  voltage  transfer 
function,  say  T(s),  must  satisfy  in  order  to  be  physically  realizable 
[Balabanian,  1958].  First  of  all,  no  poles  of  T(s)  can  lie  in  the 
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right  half  plane.  In  addition,  no  poles  of  T(s)  can  lie  at  zero  or 
infinity.  The  zeros,  however,  can  lie  anywhere  in  the  s  plane.  If 
all  the  zeros  of  ?(s)  are  confined  to  the  left  half  plane,  t(s)  is 
called  a  minimum-phase  function;  otherwise  it  is  called  a  nonminimum- 
phase  function.  The  conditions  imposed  on  T(s)  depend  on  the  network 
topology  desired.  Here  we  will  seek  a  transformerless  realization. 

Also,  we  will  prefer  an  unbalanced  (common  ground)  topology  over  a 
balanced  (no  common  input  and  output  terminals)  realization,  whenever 
possible.  Following  Balabanian  [1958]  we  summarize  these  conditions 
as  follows. 

For  unbalanced  two-ports  with  no  mutual  inductance: 

1.  Transmission  zeros  may  lie  anywhere  except  on  the  positive 
real  axis. 

2.  For  real  positive  values  of  s,  the  value  of  the  transfer 
function  lies  between  0  and  1.  The  maximum  value  of  unity 
can  be  attained  only  at  zero  or  infinity,  or  both. 

3.  The  numerator  coefficients  of  the  transfer  function  are  positive 
(or  zero)  and  no  greater  than  the  corresponding  denominator 
coefficients,  assuming  no  cancellation  of  factors. 

4.  For  a  ladder  network  realization,  which  is  a  special  form  of 
an  unbalanced  network,  the  transmission  zeros  cannot  lie  in 
the  right  half  plane. 

For  balanced  two-ports  with  no  mutual  inductance: 

1.  Transmission  zeros  may  lie  anywhere,  including  the  positive 
real  axis. 

2.  For  real  positive  values  of  s  the  value  of  the  transfer  function 
lies  between  -1  and  +1.  The  extremities  of  the  range  can 

be  attained  only  at  zero  or  infinity,  or  both. 

3.  The  numerator  coefficients  of  the  transfer  function  may  be 
negative,  but  they  are  no  greater  in  magnitude  than  the 
corresponding  denominator  coefficients,  assuming  no  cancellation 
of  factors. 
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As  discussed  above,  for  physical  realizability  one  requires  that 
the  pole  of  (4.33)  lie  in  the  left  half  of  the  complex  frequency 
plane.  This  is  satisfied  only  if  the  denominator  coefficients  of 
(4.33)  are  of  equal  sign.  In  the  unmodified  pole  module  case  this 
imposes  the  following  restrictions  on  the  location  of  the  residue  a  : 


-s.  <  -L 


(4.36) 


The  geometrical  interpretation  of  this  condition  is  given  in  Figure 

4.13(a).  The  transfer  function  is  realizable  if  the  a  lies  in  the 

n 

shaded  region  of  the  figure.  If  the  modified  pole  modules  are  employed, 
and  (4.11)  is  satisfied,  the  residue  location  is  restricted  as  follows: 


- - —  <  -p  <  q 

2Q2  -  1  dn  ~  n 
n 


(4.37) 


The  realizability  region  associated  with  this  condition  is  the  shaded 
region  in  Figure  4.13(b). 

The  comparison  of  Figures  4.5  and  4.13  reveals  that  the  voltage 
transfer  function  is  always  realizable  if  the  associated  pole-pair 
immittance  function  is  realizable.  If  the  pole-pair  immittance  is 
initially  non-PR,  i.e. ,  the  residue  falls  outside  the  regions  I  and 
II  of  Figure  4.5  so  that  the  padding  procedure  described  in  Section  4.3 
is  employed,  the  associated  voltage  transfer  function  is  still  realizable 
in  at  least  one  of  the  two  forms:  modified  or  unmodified.  There  is  a 
region  in  the  a^-plane  where  the  realizability  regions  of  Figures 


4.13(a)  and  4.13(b)  overlap,  so  that  both  forms  could  be  realized. 

We  have  found  that  for  all  dominant  poles  of  the  thin-wire  dipole  and 
loop  antennas  the  residues  a^  fall  into  that  part  of  the  realizability 
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region  lying  in  the  first  quadrant  of  Figure  4.13(b),  so  that  the 
modified  form  of  the  voltage  transfer  function  is  realizable. 


For  notational  convenience  the  index  n  on  the  coefficients  of 
(4.33)  is  dropped  throughout  the  rest  of  this  section.  Also  it  is 
convenient  to  introduce  a  scaling  factor  k  and  rewrite  (4.33)  as 


1  kAs  +  kB  _  *12  t  ^12  (4.40) 

k  Cs+D  Sn  y22 


where  z  and  y  denote,  respectively,  the  open-circuit  impedance 
parameters  and  the  short-circuit  parameters  of  the  two-port  [Balabanian, 
1958],  If  (4.40)  is  realizable,  the  coefficients  C  and  D  are  of  the 


same  sign,  as  discussed  above.  It  will  be  assumed  here  that  they  are 


positive,  which  can  always  be  accomplished  by  multiplying  the  numerator 
and  the  denominator  of  (4.40)  by  -1.  The  following  discussion  can  also 


be  simplified  if  the  sign  of  k  is  chosen  such  that  kA  is  positive. 

Then  it  is  only  necessary  to  consider  separately  two  situations: 

kB > 0  (T  (s)  is  a  minimum-phase  function),  and  kB <  0  (T  (s)  is  a 
—  n  n 

nonminimum-phase  function).  In  each  case  two  circuit  topologies  are 


obtained  depending  on  the  choice  of  the  impedance  or  admittance 
parameters  representation  in  (4.40). 


Case  1 .  kB  >,  0 

In  this  case  one  can  make  the  associations 

z^2  ■  kAs  +  kB  (4.41a) 

and 

z'u  -  Cs+D  -  (kAs  +  kB)  +  (C-kA)s  +  (D  -  kB)  .  (4.41b) 
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To  this  set  of  parameters  corresponds  the  RL  ladder  circuit  shown  in 
Figure  4.14(a).  This  circuit  is  physically  realizable  if  the  value 
of  the  unspecified  constant  k  is  restricted  to 

fc.  D 

k  <_  min  I—  ,  -j  .  (4.42) 

It  is  easy  to  see  that  by  choosing  the  maximum  allowed  value  of  k 
one  can  eliminate  either  L,  or  R  from  the  circuit  in  Figure  4.14(a). 

A  A 

If  the  admittance  representation  is  used  in  (4.40)  one  can 

write 

-y^2  =  kAs + kB  (4.43a) 

and 

y22  =  Cs  +  D  =  (kAs  +  kB)  +  (C-kA)s  +  (D  -  kB)  .  (4.43b) 

This  set  of  parameters  can  be  realized  by  the  RC  ladder  network  shown 
in  Figure  4.14(b).  This  network  is  physically  realizable  xf  k  satisfies 
(4.42).  Again,  by  a  judicious  choice  of  the  value  of  k  one  can  eliminate 
either  Cgor  R^. 

Case  2.  kB  <  0 

Since  in  this  case  the  voltage  transfer  function  has  a  positive 

real  zero,  it  cannot  be  realized  in  an  unbalanced  form  [ Balabanian , 

1958].  Instead,  one  must  resort  to  the  symmetrical  lattice  topology 

shown  in  Figure  4.15(a).  If  this  network  is  unwrapped  to  the  equivalent 

form  of  Figure  4.15(b),  it  is  recognized  that  the  symmetrical  lattice 

is  identical  with  a  bridge  network. 

The  lattice  arm  impedances  Z  and  Z,  can  be  expressed  in  terms 

a  b 

of  the  open-circuit  impedance  parameters  as 
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4 

4 


(b) 


Figure  4.14  Ladder  network  realizations  of  (a)  the  RL  and  (b)  RC 

type  for  the  pole-pair  voltage  transfer  function  —  negative  transmission 

zero  case. 
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Figure  4.15  (a)  Standard  and  (b)  bridge  f 

symmetrical  lattice. 
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(4.44a) 


*a  *  *11“  *12 

and 

Zb  *  z ^  + z12  .  (4.44b) 

Alternatively,  the  lattice  arm  admittances  » 1/Za  and  Y^  ■  1/Z^  can 
be  expressed  in  terms  of  the  short-circuit  parameters  as 

and 

Using  (4.41)  in  (4.44) 


Ya  =  y22  ”  y12 


Yb“y22  +  y12  * 


gives 


(4.45a) 


(4.45b) 


Z  *  (C-kA)s  +  (D-kB)  (4.46a) 

a 

and 

1  =  (C  +  kA)s  +  (D  +  kB)  .  (4.46b) 

□ 

This  set  of  parameters  can  be  realized  by  the  network  shown  in  Figure 
4.16(a).  For  physical  realizability  it  is  required  that 


k  <  min 


C 

A’ 


(4.47) 


By  giving  k  the  maximum  value  allowed  by  (4.47)  one  can  make  either 
La  or  disappear. 

In  like  fashion,  using  (4.43)  in  (4.45)  gives 


Y  -  (C  +  kA)s  +  (D  +  kB)  (4.48a) 

a 

and 

Y.  -  (C-kA)s  +  (D-kB)  .  (4.48b) 

D 

To  this  set  of  parameters  corresponds  the  network  shown  in  Figure 
4.16(b).  The  value  of  k  is  still  restricted  by  (4.47).  Again,  either 
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LA=C‘kA  RA=D'kA 


LB=C+kA/ 


CbmXrasM§ 


p  .  ■ 

"B  D-kB 


(b) 


Figure  4.16  Symmetrical  lattice  realizations  of  (a)  the  RL  and  (b) 
RC  type  for  the  pole-pair  voltage  transfer  function  —  positive 
transmission  zero  case. 
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R^  or  Cg  can  be  removed  by  a  judicious  choice  of  the  value  of  k. 

The  circuits  given  above  do  not  realize  the  multiplicative  factor 
1/k  in  (4.40).  The  gain  lost  can  be  restored  in  the  controlled 
source  stage  following  each  voltage  transfer  network.  The  amplification 
factor  of  this  controlled  source  must  have  the  value  p  •  1/k. 

Still  different  realizations  in  the  case  when  kB < 0  can  be 
obtained  by  writing  (4.40)  as  a  product  of  an  all-pass  function^ 
and  a  minimum-phase  function  as  follows 

t  (e)  m  —  M§±M  kAs  -  kB  4nx 

VS}  k  kA  -kB  C  +D  • 

Now  the  all-pass  function  can  be  realized  as  a  simple  constant- 
2) 

resistance  LC  lattice,  and  the  minimum-phase  function  can  be  synthesized 
as  a  constant-resistance  ladder  [ Balabanian ,  1958].  A  cascade  connection 
of  these  two-ports  realizes  (4.49).  This  realization,  however,  is 
more  complicated  than  the  networks  shown  in  Figure  4.16. 

4,5  Conclusions 

In  this  chapter  simple  procedures  are  developed  for  the  synthesis 
of  physically  realizable  active  equivalent  circuits  for  energy  collecting 
structures  from  their  SEM  description.  The  general  network  topology 
is  that  shown  in  Figure  4.1  for  the  short-circuit  boundary  value  problem 
(the  Norton  equivalent),  or  that  shown  in  Figure  4,2  for  the  open- 
circuit  boundary  value  problem  (the  Thevenin  equivalent).  These  networks 

^An  all-pass  function  is  defined  as  a  transfer  function  whose  zeros 
are  all  in  the  right  half  plane  and  whose  poles  are  the  negatives 
of  its  zeros. 

2) 

A  constant-resistance  two-port  is  defined  as  one  whose  driving- 
point  impedance  is  equal  to  a  constant  R  when  the  two-port  is 
terminated  in  a  resistance  R. 
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have  the  form  of  a  connection  of  simple  modules  associated  with  conjugate 
pole  pairs  of  the  structure.  These  pole-pair  circuits  are  synthesized 
one  at  a  time,  so  if  the  need  arises  the  network  can  be  easily  expanded 
by  adding  more  pole-pair  modules  without  repreating  the  whole  synthesis 
process.  By  using  the  network  topologies  of  Figures  4.1  and  4.2  the 
synthesis  process  if  reduced  to  the  realization  of  a  driving-point 
function  biquadratic  in  s  and  a  voltage  transfer  function  bilinear  in  s. 

The  realizability  conditions  and  explicit  circuit  realizations  are 
given  for  these  functions.  If  the  pole-pair  driving-point  function  is 
initially  nonrealizable,  a  simple  padding  procedure  is  described  to 
synthesize  an  approximating  network.  It  is  shown  that  the  pole-pair 
voltage  transfer  function  can  be  always  realized  either  in  the  form 
of  a  ladder  circuit  or  of  a  symmetrical  lattice  circuit.  Unfortunately, 
the  applicability  of  each  of  the  two  topologies  is  dependent  on  the 
direction  of  arrival  of  the  incident  plane  wave,  so  that  the  circuit 
topology  may  change  when  excitation  conditions  are  changed.  The 
equivalent  network  does  not  depend  on  the  time  history  carried  by  the 
incident  field,  however.  Also,  the  autonomic  voltage  source  is 
connected  to  only  one  port,  which  is  a  desirable  feature.  As  the 
Figures  4.1  and  4.2  indicate,  the  controlled  sources  could  not  be 
avoided.  They  can  be  implemented  in  practice  with  active  devices  which, 
however,  usually  can  not  sustain  very  high  voltages. 

In  the  next  chapters  the  procedures  developed  here  are  used  to 
synthesize  practical  equivalent  circuits  for  the  loop,  dipole,  and 
spherical  antenna. 
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CHAPTER  5 


EQUIVALENT  CIRCUIT  SYNTHESIS  FOR  A  CIRCULAR  LOOP 

5.1  Introduction 

In  this  chapter  the  procedures  developed  in  Chapters  2-4  are 
used  to  synthesize  an  SEM  Nor ton- type,  dominant  pole-pair  equivalent 
circuit  for  a  thin-wire  loop  antenna.  The  thin-wire  circular  loop 
is  one  of  the  few  structures  for  which  an  explicit  solution  for  the 
current  can  be  found  (subject  to  the  thin-wire  approximation).  The 
development  presented  here  differs  in  details  from  the  general  theory 
developed  in  Chapter  3  because  of  the  occurrence  of  the  eigenvalue 
degeneracy  in  the  present  case. 

The  performance  of  the  developed  SEM  equivalent  circuits  is 
compared  with  responses  obtained  by  the  classical  frequency  domain — 
inverse  Fourier  transform  approach.  The  application  of  these  circuits 
to  the  analysis  of  antennas  with  nonlinear  loads  is  demonstrated  in 
Appendix  B. 

5.2  Preliminary  Theory 

The  geometry  of  the  loop  is  defined  in  Figure  5,1.  It  is  assumed 

that  the  radius  of  the  wire  a  is  much  smaller  than  the  radius  of  the 

loop  b(a«b)  and  that  oj  a/c  «  1,  where  co  *  go  is  the  largest 
r  max  max  ° 

significant  spectral  component  of  the  excitation.  As  a  consequence, 
it  can  be  assumed  that  there  is  only  a  ({(-component  of  the  surface 

•w 

current  density,  t.e.,  J(r,s)  * J^(r,s)0,  where  $  denotes  the  unit 
vector  in  the  $  direction,  and  that  the  current  density  is  uniform  around 

the  wire,  i.e.,  J ,  (r,s)  ■  J .  (0,s)  .  The  total  current  can  be  defined  as 

W  <P 
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I(<f>,s)  ■  2iraJ ,(4»*s)  . 


(5.1) 


'V 

The  EFIE  for  the  total  current  becomes  one-dimensional  with  the 
kernel  defined  by 


Z(«M'.s) 


Z  r— 
O 

T  Ls 


S  COS  ($-<$') 


s  3$' 


3  1 

-sR 

e 

3£j 

4trbR 

(5.2) 


s  =  yb  is  the  normalized  complex  frequency,  y  being  the  complex 
propagation  constant,  and  R  is  the  normalized  distance  between  the 
source  point  and  the  observation  point  on  the  surface  of  the  wire 
(see  Figure  5.1)  given  by 


with 


S./4siny*^l)  +  (|): 


A  ■  2a  sin(^)  . 


(5.3) 


(5.4) 


The  suitable  symmetric  product  is  defined  as 

TT 


<[  ];[  ]>  5 


[  ]  •  [  ]b  d<p  . 


(5.5) 


j 

-IT 


The  eigenvalues  of  the  EFIE  can  be  found  as  [Harrington,  1968; 
King,  1969] 


where 


a  (s)  *  -Is 
n  J 


An(s)  *  IT  an(s) 


K  .  (s)  +  K  ,  (s)  , 

— - 2^—  +  (J)  ?„(.) 


(5.6) 


(5.7) 


and 
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Kn(.) 


n  t 


Jn(<M>')  -sR 


d(4H>’)  • 


(5.8) 


-7T  -IT 


As  shown  by  Wu  [1962],  eq.  (5.8)  may  be  approximately  written  in 
terms  of  integrals  of  Anger-Weber  functions  for  unrestricted  n,  as 


Vs)  =  7  -i  |  [no(z)+j  Jo(z)]dz 


(5.9a) 


K  (s)  =  -jl  (M)k  f^)  +  cl  -  ~ 

n  7T|_ol'b;o'-b;  nj  2 


[Sn(z)  +j  J2n(z)_  dz 


(5.9b) 


C  =  £n(4F)  -2  l  (2m+l)' 
m=0 


(5.10) 


where  C,  =  -j2s  and  I  and  K  are  the  modified  Bessel  functions  of  the 
o  o 

first  and  second  kind,  respectively,  and  F= 1.781...  is  the  Euler's 
constant.  This  representation  was  originally  obtained  by  Wu  for  the  real 
frequency  case  (s*jm);  it  was  generalized  to  the  complex  frequency 
case  by  Umashankar  and  Wilton  [1974]. 

The  eigenfunctions  associated  with  the  eigenvalues  (5.6)  are 


I  e(<t>> 

n.S 


cos (n$) 


sin(n<f>) 


(5.11) 


so  that  there  is  an  eigenvalue  degeneracy  for  n  > 1. 

The  current  I(s)  =  I(4>  ,s)  at  the  location  0  *  cf>  can  be  expanded 

8  8 

in  terms  of  the  eigenfunctions  as 


I(s)  »  ^  -4rr  - * - 


Zo  ao(s)  <1  o(cf>)  ,1  (<*>)> 

o,e  o,e 


.  <1  _(*),E*(*,s)> 

+  I  I  — — — - ^ i -  I  ) 

n-l5-n<«)  <1  G(*).I  „(♦>>  n’°  8 

n,G  n,a 


(5.12) 
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where  the  subscript  a  indicates  e  (for  "even")  or  o  (for  "odd"). 
Noting  that 


and 


<1  (<f>),I  ($)>  -  27Tb 

o,e  o,e 


<I_  e(<f»).I„  e(<P)>  =  Kb 


n»0  n»o 

one  can  rewrite  (5.12)  as 


(5.13a) 


(5.13b) 


I(s) 


ttZ 


1  <1  a(4>),Ei(4>,s)> 


aQ(s)  n,e 


.  .  ±  1 

n“i  a  n 


(5.14) 


5.3  Derivation  of  the  Driving- Point  Admittance  of  the  Loop  Antenna 

For  the  driving-point  admittance  computation  a  delta  function 

generator  located  at  $  *  <J>  is  assumed,  with  the  field  defined  by 

§ 


Ej(<fr,s)  -  5(<P-$g)  . 


(5.15) 

Using  (5.15)  in  (5.14)  the  admittance  Y(s)  ■  I(s) /V(s)  can  be  written 


as 


Y(s) 


ttZ 


+ ^  1 4* 


(5.16) 


Umashankar  and  Wilton  [1974]  have  shown  that  the  term  l/a^s) 


can  be  expanded  in  a  residue  series 


.  R  00 

_-J—  „  _s°  +  y 

in(s)  3  i-1 


1  .  R* 

ni  ni 


s  -  s 


ni 


s  -  s 


ni 


(5.17) 


where  the  term  associated  with  the  pole  at  s  *  0  is  present  only  for 
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n-0,  i.e.,  R  *0  for  n  >  1,  and  the  residues  R  .  are  given  by 

nO  —  nl  - 


s  -  s 


R  .  ”  2-im  — . 7^7 

ni  .  a_(s) 

s  sni  n 


dan(sni) 


ni 

-1 


ds 


(5.18) 


In  actual  computations  this  series  must  be  truncated  after  some 
finite  number  of  terms  (say  M) .  Umashankar  and  Wilton  [1974]  have 
found  that  the  truncated  series  (5.17)  converged  very  slowly  at  low 
frequencies  and  they  were  forced,  in  effect,  to  use  the  modified  pole 
expansion  (see  Section  2.3) 


an(s) 


(5.19) 


They  also  observed  that  the  poles  belonging  to  a  given  eigenmode  could 
be  divided  into  three  categories: 

1.  Type  I  -  A  single  pole  of  each  eigenset  nearest  the  jio-axis 

at  approximately  oj*n.  This  pole  is  the  principal  contributor 
to  the  response  of  the  loop  -  it  is  the  dominant  pole  (see 
Section  2.7). 

2.  Type  II  -  A  set  of  (n+1)  poles  (including  conjugate  pairs) 
which  lie  roughly  along  an  elliptic  arc  centered  at  s  *  0 
and  with  a  semi-major  axis  somewhat  larger  than  n. 

3.  Type  III  -  A  group  of  an  infinite  number  of  poles  lying 
almost  parallel  to  the  jw-axis.  The  spacing  of  these  poles 
is  asymptotic  to  Auj-irc/b. 

This  classification  is  illustrated  in  Figure  5.2  for  a  loop  with  the 

shape  factor  fl  ■  28,n(2Trb/a)  .  In  this  figure  all  fourth  eigenmode  poles 

with  Im{s  . } < 20b/c  are  shown  together  with  the  dominant  (Type  I) 
ni  — 

poles  of  the  eigenmodes  0-19.  More  recently  Blackburn  [1976]  (see 
also  [ Blackburn  and  Wilton ,  1978])  has  provided  an  extensive  tabulation 
of  the  poles  and  associated  residues  for  the  loop  for  several 
different  shape  parameters  ft. 
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Denoting  by  the  dominant  pole  of  the  n  eigenmode  and  by 

R'  the  associated  residue,  one  can  write 
n 


a0(s)  “  s 


(5.20a) 


R'  R'* 

-  rJ  =  s  _ 2 _ +  _x— n _ 

a  (s)  s  (s-s  )  s*(s  -  s*) 

n  n  n  n  n 


,  n  >  1 . 


(5.20b) 


To  simplify  notation,  it  is  convenient  to  define  a  normalized  residue 


n  ttZ 


(5.21) 


Using  (5.19)  and  (5.20)  in  (5.16)  yields 


?<a)  +  s  ?„<•> 

o  n=l 


(5.22) 


where  Lq = 1/Rq,  N  is  the  number  of  dominant  poles  used,  and 


Yn(s)  =  Yn  +  Yn  o(s)’ 

n  n , e  n,  o 


(5.23) 


Y  (s)  *  Y  (s)  +  Y  (s), 

n  jG  n  jO  *"n  ^  o 


a  (4>  ) 

Y  (s)  =  8  -a  f-  , 

n,o  Vs  "V 


(5.25) 


a  =  2R  ,  n  >  1 
n  n  — 


(5.26) 


As  discussed  in  Section  3.2,  the  index  -n  indicates  the  term  associated 
with  the  conjugate  pole.  Note  that  since 


e(<^  +  X„  o(4>)  S  1 
n ,  e  n,o 


(5.27) 
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for  all  values  of  <$,  it  follows  that 


and 


Y  (s)  -  Y  (s)  +  Y  (s)  -  s  — - r 

n  n,e  n,o  s  (s  -  s  ) 


Y 

n 


(s) 


Y  (s)  +  Y  (s) 
n  -n 


(5.28) 


(5.29) 


is  independent  on  the  port  location.  This  could  be  anticipated,  of 
course,  because  of  the  symmetry  of  the  loop. 

The  quality  of  the  dominant-pole  approximation  to  the  driving 
point  admittance  of  the  loop  antenna  with  fi*15.0  is  illustrated  in 
Figure  5.3  which  shows  the  comparison  of  the  behavior  of  Y(s)  on  the 
jo)-axis  computed  both  from  the  approximate  formula  (5.22),  with  N *  10, 
and  from  the  exact  expression  (5.16).  It  is  seen  that  the  agreement 
is  satisfactory,  in  particular  in  the  real  parts.  The  agreement  in 
the  imaginary  parts  deteriorates  for  w>3c/b. 


5.4  Derivation  of  the  Short-Circuit  Current  of  the  Loop  Antenna 

For  the  short-circuit  current  computation  at  <J >  *  $  a  plane  wave 

g 

incident  field  is  assumed  with  the  geometry  defined  in  Figure  5.4. 

The  component  of  the  incident  field  on  the  wire  is  given  by 

eJ<4>,s)  -  E*Q(4>)p(4>,s)f(s)  (5.30a) 

with 


and 


Q(<t>)  ■  sini^  cos8  s  in  (<)>-<}>  )  -  cost!  cos((j>-4>i) 
s  sin0[cos(<J>-4>  )  -  cos(4>  -  d>  )  ] 

P (4>,s)  -  e  1  g 


(5.30b) 


(5.30c) 


In  (5.30a)  E^  denotes  the  magnitude  of  the  field  strength  and  f(s) 
is  the  Laplace  transform  of  the  time  history  carried  by  the  plane  wave. 
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Figure  5.3  (a)  The  real  part  and  (b)  the  imaginary  part  of  the  spectrum 

of  the  driving-point  admittance  of  the  loop  antenna.  The  dominant  pole 
approximation  (5.22)  is  compared  with  the  exact  expression  (5.16). 
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(b) 


Figure  5.4  Circular  loop  in  a  plane-wave  incident  field,  (a)  The 
general  view;  (b)  the  plane  containing  the  loop. 


The  time  origin  is  placed  at  the  port  location  4>“4>g»  as  discussed 
in  Section  2.5.  Using  (5.30)  in  (5.14)  and  normalizing  the  current 


to  2iTbE  yields 
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with  the  coupling  coefficients  given  by 
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(5.32a) 
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(5.32b) 


In  the  dominant-pole  approximation  the  short-circuit  current  becomes 


Is<=(s)  =  Rq  Uj 
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Combining  N  conjugate  pole-pairs  in  (5.33)  results  in 
ISC(s)  *  g„f(s)  + 


n  r  ~  .  “|~ 

y  y  Y  (s)T  +  Y  (s)T*  f(s) 

4,  4  ln,cr  '  n,a  -n,o  n,a| 
n=l  a  *- 


where 


e  »  R  £im 

rt  _  ..A 


o  o  s-K)  (  8 


n(s) 

O 


cosi^  sin© 

R  o 

o  2 


(5.34) 


(5.35) 


and 


(5.36) 
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Eq.  (5.34)  can  be  rewritten  with  the  help  of  (5.32)  as 


ISC(s)  = 


g0f(s) 


N  r~ 

£  [Vn 

n*l  *- 


(s)T  +  Y  (s)T; 
n  -n 


:] 


f(s) 


(5.37) 


where  the  source  coefficients  T  are  defined  by 
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T  =  I2  (4)  )T  +  I2  (4>  )T 
n  n,e  g  n,e  n,o  g  n,o 
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>e 


■s  sin0  cosO  -  (bj ) 
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(5.38) 


As  discussed  in  Section  3.2,  eq.  (5.37)  can  be  further  transformed  to 

N 

ibC(s)  3  gof(s)  +  l  Yn(s)Tn(s)f(s)  (5.39) 

n*l 

with  the  voltage  transfer  functions  Tn(s)  defined  by  (3.24). 

The  quality  of  the  dominant  pole  approximation  to  the  short-circuit 
current  of  a  loop  antenna  with  0*15.0  is  illustrated  in  Figure  5.5 
which  shows  the  comparison  of  the  behavior  of  I  (s)  on  the  ju-axis 
computed  both  from  the  approximate  formula  (5.39),  with  N *  10,  and 
from  the  exact  expression  (5.31).  It  is  seen  that  the  agreement  is 
satisfactory. 

The  expressions  (5.22)  for  the  driving-point  admittance,  and 
(5.39)  for  the  short-circuit  current,  are  now  in  a  form  suitable  for 
the  application  of  the  synthesis  procedures  developed  in  Chapter  4. 


5.5  Equivalent  Circuits  for  the  Loop  Antenna  With  SI  *  15 

In  this  section  Norton-type  equivalent  circuits  are  presented  for 
a  loop  antenna  with  the  shape  factor  (1*15.  The  general  layout  of  the 
equivalent  circuit  is  shown  in  Figure  4.1.  It  was  found  that  only  10 
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Figure  5.5  (a)  The  real  part  and  (b)  the  imaginary  part  of  the  spectrum 

of  the  short-circuit  current  of  the  loop  antenna.  The  dominant  pole 
approximation  (5.31)  is  compared  to  the  exact  expression  (5.39). 
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or  less  dominant  pole-pair  modules  model  the  structure  adequately  for 
the  excitations  used.  The  SEM  data  employed  (the  poles  and  the 
admittance  residues)  were  taken  from  [Blackburn,  1976]  and  are  listed 
in  Table  5.1.  Also  indicated  in  this  table  for  each  pole  are  the  quality 
factor  Q^,  the  residue  location — as  referred  to  Figure  4.5,  and  the  padding 
necessary  to  restore  the  PR-ness  of  the  associated  pole-pair  admittance 
Y^Cs).  The  column  denoted  "%"  gives  the  ratio  (in  %)  of  the  padding 
conductance  Gn  to  the  maximum  value  of  the  real  part  of  Yn(ju)). 

A 

It  is  evident  from  Table  5.1  that  only  Y^(s)  is  initially  PR  and  can  be 
realized  by  one  of  the  ladder  circuits  shown  in  Figure  4.7.  For  all 
other  pole-pair  admittances  the  residues  lie  outside  the  regions  I  and 
II,  in  the  area  denoted  "case  A"  (see  Figure  4.5),  so  that  padding  is 
required.  However,  the  amount  of  padding  is  insignificant  as  compared 

/A 

to  the  value  of  the  real  part  of  Y^(juj)  near  the  resonance.  The  padded 
pole-pair  admittances  can  be  realized  by  the  simplified  Bott-Duffin 
modules,  as  discussed  in  Section  4.3.  The  element  values  of  the  first 
10  dominant  pole-pair  modules  are  listed  in  Table  5.2.  Similar  driving- 
point  admittance  network  for  the  loop  was  synthesized  previously  by 
Streable  and  Pearson  [1981]. 

The  current  response  of  the  equivalent  network  defined  in  Table 
5.2  due  to  a  Gaussian  voltage  pulse,  computed  by  the  SCEPTRE  [Jensen 
and  MaNamee 3  1976]  circuit  analysis  program,  is  shown  in  Figure  5.6. 

This  response  is  compared  to  a  waveform  obtained  from  the  TWTD 
[VanBlaricum  and  Miller ,  1972]  computer  code,  which  solves  an  integral 
equation  for  the  current  in  the  time  domain.  It  requires,  however, 
that  the  structure  be  modeled  by  a  connection  of  straight  wire  segments. 
Only  about  50  segments  were  employed  to  model  the  loop  for  the  TWTD  run 
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The  Poles  s  and  residues  a  are  normalized  to  c/b 


Table  5.2  Element  values  of  the  dominant  pole-pair  admittance 
networks  for  a  loop  antenna  with  Q*15.*^ 


LADDER  CIRCUIT  (Figure  4.7(a)) 


n 

C1/b[pF/m] 

i^tn] 

L^/btyH/m] 

R2(kft] 

n 

3.3034 

_ 

47.374 

3.0865 

_ j 

10.164 

BOTT-DUFFIN  CIRCUIT-CASE  A  (Figure  4.12(a)) 

n 

Co/b[pF/m] 

L^/bfpH/m] 

C1/b[pF/m] 

Rjtkfl] 

H 

0.9333 

2.8764 

44.3967 

16.144 

0.5095 

2.7291 

2.8702 

20.727 

H 

0.3355 

2.6209 

1.0626 

24.225 

5 

0.2443 

2.5355 

0.5527 

26.998 

6 

0.1893 

2.4648 

0.3384 

29.254 

7 

0.1531 

2.4045 

0.2284 

31.125 

8 

0.1276 

2.352 

0.1645 

32.701 

0.1089 

2.3055 

0.1242 

35.042 

n 

0.0946 

2.2637 

0.0971 

35.197 

^The  static  inductance  L  /b « 7.215  [yH/m]. 

o 
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with  the  gap  width  equal  to  the  length  of  one  segment.  This  crude 
modeling  is  a  possible  explanation  of  the  quite  significant  dis¬ 
agreement  (about  15%  in  the  worse  case)  of  the  SCEPTRE  and  TWTD 


curves . 

Because  of  the  symmetry  of  the  loop,  the  driving-point  part  of  the 

equivalent  network  does  not  depend  on  the  gap  location  ■$  ;  neither 

does  it  depend  on  the  direction  of  propagation  of  the  plane  wave  incident 

field,  of  course.  The  transfer  function  part  of  the  equivalent  network 

does  depend  on  the  angles  of  incidence  and  the  polarization  of  the 

incident  plane  wave,  as  discussed  in  Chapter  4.  Not  only  the  element 

values  of  the  active  part  of  the  network  change  when  the  direction  of 

incidence  is  changed,  but  also  the  network  topology  can  be  affected. 

However,  it  is  shown  in  Chapter  4  that  at  most  two  network  topologies 

can  result:  a  ladder  network  or  a  symmetrical  lattice  network.  For 

both  cases  explicit  expressions  for  the  element  values  in  terms  of  the 

poles  s,  the  admittance  residues  a  and  the  source  coefficients  T 

n  n 

are  given  in  Chapter  4.  The  expression  for  T  in  the  loop  case  is 

n 

derived  in  Section  5.4,  eq.  (5.38).  The  topologies  and  element  values 
of  the  corresponding  transfer  networks  are  listed  in  Tables  5. 3-5. 6. 

Note  that  we  have  chosen  to  use  the  RC  networks  rather  than  the  RL 
networks.  This  choice  was  made  because  the  element  values  of  the  RC 
networks  turned  out  to  be  more  realistic  than  the  element  values  of  the 
equivalent  RL  networks. 

The  short-circuit  current  responses  of  the  complete  equivalent 
networks  computed  by  SCEPTRE  are  compared  in  Figures  5.7-5.12  with 
responses  obtained  by  other  means  for  different  excitation  conditions. 
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Table  5.4  Element  values  of  the  dominant  pole-pair  voltage  transfer  networks  for  a  loop  antenna 
(0  =  15)  with  a  port  at  .f)g  =  90° ,  excited  by  a  plane  wave  with  the  angles  of  incidence  0  =  90  , 
iu  =0°.  and  v-  180°.*) 


The  transfer  networks  associated  with  poles  s  for  n= 1,3,5,...  are  not  present. 
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:  5.7  Short-circuit  current  excited  in  the  circular  loop  at  0=0°  by  a  step  function 
wave.  SCEPTRE  and  TWTD  results  are  compared. 


:  5.9  Short-circuit  current  excited  in  the  circular  loop  at  <p  = 180°  by  a  step  function 
wave.  SCEPTRE  and  TWTD  results  are  compared. 


orc-circuit  current  excited  in  the  circular  loop  at  >{>  =  0°  by  a  double-exponential 
wave.  SCEPTRE  and  FFT  results  are  compared. 


!  5.11  Short-circuit  current  excited  in  the  circular  loop  at  4>  =  90°  by  a  double-exponential 
wave.  SCEPTRE  and  FFT  results  are  compared. 


*  5.12  Short-circuit  current  excited  in  the  circular  loop  at  <|>  =  180°  by  a  double-exponential 
wave.  SCEPTRE  and  FFT  results  are  compared. 


SCEPTRE  responses  are  compared  with  TWTD  results  in  Figures 
5. 7-5.9.  The  time  history  carried  by  the  plane  wave  is  a  step  function. 
As  in  the  case  of  the  admittance  network,  there  is  a  significant 
discrepancy  between  these  two  curves.  Again,  one  can  attribute  this 
discrepancy  to  the  relatively  crude  modeling  of  the  loop  by  the 
straight  pieces  of  wire.  It  is  of  interest  to  note  that  the  current 
waveform  approaches  a  constant  non-zero  value  for  late  times.  This 
phenomenon  is  possible  because  there  is  a  d-c  (direct-current) 
component  in  the  incident  waveform  spectrum  and  the  loop  is  perfectly 
conducting. 

The  circuit  responses  are  compared  with  the  results  of  the 
classical  frequency  domain — Fourier  transform  method  in  Figures  5.10- 
5.12.  The  time  history  of  the  incident  wave  is  a  double  exponential 
function.  In  computing  the  frequency  domain  response  the  series 
(5,31)  for  the  short-circuit  current  was  evaluated  for  s = ju)  with  30 
terms  included.  The  subroutine  developed  by  Blackburn  et  cl.  [1978] 
was  used  to  evaluate  the  a^jw)  factors.  The  time  domain  response  was 
obtained  by  applying  the  Fast  Fourier  Transform  (FFT)  algorithm 
[ Brigham ,  1974].  It  is  seen  fromFigures  5.10-5.12  that  the  SCEPTRE  and 
FFT  responses  agree  favorably  (the  maximum  discrepancy  is  less  than  8%). 

5.6  Conclusions 

The  responses  of  the  dominant  pole-pair  equivalent  circuits  for 
a  thin-wire  loop  antenna  agree  favorably  with  the  exact  responses 
obtained  from  the  frequency  domain — inverse  FFT  approach.  The  amount 
of  the  real-part  padding  necessary  for  physical  realizability  of  the 
pole  pair  admittances  is  negligible  in  comparison  to  the  value  of  the 
real  part  at  the  resonance. 


CHAPTER  6 


EQUIVALENT  CIRCUIT  SYNTHESIS  FOR  A  CYLINDRICAL  DIPOLE 

6.1  Introduction 

In  this  chapter  the  procedures  developed  in  Chapters  2-4  are  applied 
to  the  synthesis  of  the  dominant  pole-pair,  Norton-based,  equivalent 
circuits  for  a  thin  cylindrical  dipole  antenna.  The  cylindrical  dipole 
is  a  member  of  the  class  of  structures  for  vhich  an  exact  analytical 
solution  is  not  available.  Instead,  the  SEM  data  are  derived  numerically 
by  the  MoM  technique.  Simplified  equivalent  circuits  based  on  the 
sinusoidal  mode  approximation  are  also  derived.  For  a  symmetrical  dipole 
antenna  or  a  cylindrical  post  over  a  ground  plane,  particularly  simple 
and  yet  reasonably  accurate  equivalent  circuits  result.  The  responses 
of  the  SEM  equivalent  circuits  are  compared  with  the  results  of  the 
classical  frequency  domain — inverse  FFT  approach.  An  application  of  the 
derived  equivalent  circuits  to  the  analysis  of  antennas  with  nonlinear 
loads  is  illustrated  in  Appendix  B. 

6.2  Preliminary  Theory 

The  geometry  of  the  dipole  is  defined  in  Figure  6.1.  It  is  assumed 

that  the  cylinder  Is  slender  (a  «  l)  and  electrically  thin,  i.e., 

a^^a/c  «  1,  where  u  ■  is  the  largest  significant  spectral  component 

of  the  excitation.  As  a  consequence,  it  can  be  assumed  that  there  is 

only  a  z-component  of  the  surface  current  density,  i.e.,  J(r,s) * Jz(r,s)z, 

where  z  denotes  the  unit  vector  in  the  z  direction.  It  can  be  further 

assumed  that  the  current  density  is  uniform  around  the  cylinder,  i.e., 

J  (r,s)=J  (z,s),  so  that  the  total  current  is  given  by 
2  2 
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(6.1) 


I(z,s)  *  2na  J  (z,s)  . 

z 

The  current  on  the  dipole  can  be  found  from  the  so-called  "Pockling ton's 
equation" 


where 


l 

I(z,s)K(z,z ' ,s)dz ’ 

0 


E^(z,s) 


(6.2) 


with 


K(z,z’ ,s) 


2tt 


7 r 

. 

-7 T 


d<t>' 


(6.3) 


R  “  /(z-z')2  +  4a2  sin2(4>’/2) . 


(6.4) 


E*(z,s)  denotes  the  z  component  of  the  incident  electric  field  and 

Y  *  s/c  is  the  complex  propagation  constant. 

In  order  to  find  the  short-circuit  SEM  parameters  of  the  dipole 

the  gap  is  closed  (Figure  6.1)  and  the  integro-differential  operator 

in  (6.2)  is  matricized  by  using  the  method  of  moments  (MoM)  technique 

[Harrington,  1968],  In  a  simple  variation  of  this  method,  as  applied 

to  the  present  problem,  the  cylinder  is  divided  into  a  number,  say  M, 

of  equal  size  zones  of  width  A  »  £/M  and  the  current  distribution  within 

a  zone  is  approximated  by  a  constant.  By  replacing  the  differential 

operator  in  (6.2)  by  finite  differences  f Wilton  and  Butler ,  1976]  and 

enforcing  the  integral  equation  at  the  match  points  z^,  k*l,2,...,N 

(collocation  testing),  a  linear  matrix  equation  is  obtained  for  the 

expansion  coefficients  of  the  current.  The  natural  frequencies  (the 

poles)  are  found  as  the  zeros  of  the  determinant  of  the  MoM  matrix. 

At  these  freqeuncies  the  homogeneous  eq.  (6.2)  has  nontrivial  solutions — 

the  natural  modes.  The  natural  mode  associated  with  the  pole  s  can 

n 

be  expressed  as 
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(6.5) 


V2>  *  X  !n.k  V2> 

k*l 


where  ^  are  the  current  expansion  coefficients  obtained  from  the 
MoM  solution  and  P^(z)  are  the  pulse  expansion  functions  defined  by 


pk(z)  ■ 


fl  ,  -  A/2  £  z  +  A/2 


|0  ,  otherwise 


(6.6) 


From  the  MoM  matrix  and  the  natural  modes  one  can  compute  the  normali¬ 
zation  constants  8  according  to  eq.  (2.20). 
n 

6.3  Derivation  of  the  Driving-Point  Admittance  of  the  Dipole  Antenna 

The  driving-point  admittance  of  the  dipole  is  determined  by  computing 
the  current  I(s)  in  the  gap  region  due  to  a  voltage  V(s)  * V^f (s)  impressed 
across  the  gap,  as  discussed  in  Section  3.2.  The  incident  field  due 
to  this  voltage  can  be  approximated  by 


E^(z,s)  s  -2-  Pg(z)f(s) 


(6.7) 


where  k  =*  g  is  the  index  of  the  gap  zone.  The  coupling  coefficients 
can  be  found  from  (see  Section  3.2) 


I  (z)  P  (z)dz  s  V  1 
n  A  g  o  n,g 


(6.8) 


and  the  admittance  residues  from 


a  s  3  i' 

n  n  n,g 


(6.9) 


The  driving-point  admittance  Y(s)  » I(s)/V(s)  has  the  form 


N  „ 

Y(s)  3  l  Yn(s) 
n“l 


(6.10) 


146 


where  Yn(s)  are  the  pole-pair  admittances  defined  in  (3.16)  and  N  is 
the  number  of  dominant  poles  employed.  Note  that  the  inductive  term 
of  eq.  (3.15)  is  not  present  in  the  dipole  case. 

The  poles  and  natural  modes  of  a  thin  cylinder  with  the  shape 
factor  £2  *  2£n(£/a)  *  10.6  were  first  found  by  Tesahe  [1973]  by  the 
procedure  described  in  the  previous  section.  In  this  work,  however, 
the  natural  modes  are  given  in  a  graphical  form  only.  Since  the  numerical 
values  of  the  coefficients  1^  ^  were  needed  for  the  circuit  synthesis 
purposes,  Tesche's  data  were  reproduced  using  72  pulses  in  the  current 
expansion.  The  gap  width  was  taken  to  be  one  zone.  The  results  are 
summarized  in  Figure  6.2  and  Table  6.1.  It  can  be  seen  from  Figure 
6.2  that  the  poles  lie  in  layers  roughly  parallel  to  the  jw-axis.  The 
first  layer  poles  are  the  dominant  ones.  Also  shown  in  this  figure 
is  a  grouping  of  poles  according  to  eigensets,  conjectured  by  D.  R.  Wilton 
[Streable  and  Pearson ,  1981],  Note  that,  unlike  in  the  loop  case,  there 
are  only  a  finite  number  of  poles  belonging  to  each  eigenmode  in  this 
representation. 

The  quality  of  the  dominant-pole  approximation  to  the  driving- 
point  admittance  of  a  cylindrical  dipole  antenna  with  £2*10.6  is 
illustrated  in  Figure  6.3  which  shows  the  comparison  of  the  behavior 
of  Y(s)  on  the  jw-axis  computed  from  the  approximate  formula  (6.10) 
with  N»10,  and  by  solving  the  integral  equation  (6.2)  numerically. 

It  is  seen  that  the  agreement  is  quite  good  in  the  real  parts  case 
while  it  is  not  as  good  in  the  imaginary  parts,  especially  for 
co  >  2rr  c /l. 
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Figure  6.3  (a)  The  real  part  and  (b)  the  imaginary  part  of  the  spectrum 

of  the  driving-point  admittance  of  the  dipole  antenna.  The  dominant  pole 
approximation  (6.10)  is  compared  with  the  "exact"  curves  obtained  by 
solving  the  integral  eq.  (6.2)  numerically. 
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6.4  Derivation  of  the  Short-Circuit  Current  of  the  Dipole  Antenna 
For  the  short-circuit  current  computation  the  incident  field  is 
that  due  to  a  plane  wave  carrying  a  time  history  f(t).  The  z-component 
of  the  electric  field  evaluated  on  the  z-i-jcis  is  given  by 


_ i  *  Y(z  -  z2)cos0_ 

Ez(z,s)  =  sin6  e  f(s) 


(6.11) 


with  the  time  origin  chosen  in  the  gap  zone,  as  discussed  in  Section  2.6. 
With  the  current  normalized  to  £E*,  the  coupling  coefficients  can  be 
computed  as 


sinfi 


y  (z-z  )cos9 

I  (z)e  n  ®  dz 
n 


.  .  r,  M  Y  (z,  -  z  )cos0 

Asin0  r  *  'n  k  g 

l  *  n,ke 
k»l 


(6.12) 


with  the  source  coefficients  given  by 


n»g 


(6.13) 


The  short-circuit  current  can  be  expressed  as 


N 

iSC(s)  S  l  Y  (s)T  (s)f(s)  (6.14) 

n  n 

n=l 

A 

with  the  voltage  transfer  functions  Tn(s)  defined  in  (3.24).  Note  that 
the  first  term  of  eq.  (3.23)  is  not  present  in  the  dipole  case. 

The  quality  of  the  dominant-pole  approximation  to  the  short-circuit 
current  of  a  cylindrical  dipole  antenna  with  ft -10.6  is  illustrated  in 
Figure  6.4,  in  which  the  behavior  of  I  (jtu)  computed  from  the  approximate 
formula  (6.14)  with  N -  10  is  compared  with  the  "exact"  waveform  obtained 
from  the  solution  of  the  integral  equation  (6.2)  numerically  by  the  MoM 
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I 

Figure  6.4  (a)  The  real  part  and  (b)  the  imaginary  part  of  the  spectrum 

of  the  short-circuit  current  of  the  dipole  antenna.  The  dominant  pole 
approximation  (6.14)  is  compared  with  the  "exact"  curves  obtained  by 
solving  the  integral  eq.  (6.2)  numerically. 
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technique.  The  time  history  is  the  delta  function  and  the  angle  of 
incidence  9  “60°.  It  is  seen  that  the  agreement  is  excellent,  both  in 
the  real  and  imaginary  parts. 


6.5  Equivalent  Circuits  for  the  Cylindrical  Dipole  with  ft- 10.6 

In  Table  6.1  are  listed  for  each  pole  Sn  the  quality  factor  Q^, 
the  admittance  residue  a^  and  its  location — as  referred  to  Figure  4.5, 
and  the  padding  necessary  to  restore  the  PR-ness  of  the  associated 

A 

pole-pair  admittance  Y  (s).  The  column  denoted  "%"  gives  the  ratio 

n 

(in  %)  of  the  padding  conductance  to  the  maximum  value  of  the  real 

A 

part  of  Y^Cjoi).  Data  are  provided  for  two  gap  location:  z^  =  l/2£ 
and  z^  * l/4£.  As  can  be  seen  from  this  table,  for  most  poles  the 
residues  fall  outside  the  regions  I  and  II,  in  the  area  denoted  "case 
A"  (see  Figure  4.5),  so  that  after  padding  the  Bott-Duffin  module  can 

A 

be  employed  to  realize  Y^(s).  For  two  dominant  poles  the  residues  lie 
in  region  II,  so  that  one  of  the  ladder  networks  from  Figure  4.7  can  be 
directly  used.  The  element  values  of  the  driving-point  equivalent 
networks  for  the  two  gap  locations  are  listed  in  Tables  6.2  and  6.3, 
respectively.  These  networks  were  first  obtained  by  Streable  and 
Pearson  [1981].  We  include  them  here  for  the  sake  of  completeness, 
since  they  are  part  of  the  active  equivalent  networks.  The  responses  of 
the  equivalent  networks  for  a  center-  and  a  quarter-driven  dipole  due 
to  a  Gaussian  voltage  pulse,  computed  by  the  SCEPTRE  [Jensen  and 
MoNamee ,  1976]  circuit  analysis  program,  are  compared  to  the  TWTD 
[VanBlaricum  and  Miller t  1972]  time  domain  integral  equation  solutions 
in  Figures  6.5  and  6.6,  respectively.  Except  for  the  early  time,  the 
agreement  is  seen  to  be  quite  good.  The  early  time  response  of  the  SOI 
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Table  6.2  Element  values  of  the  dominant  pole-pair  admittance 
networks  for  a  center-driven  dipole  antenna  (£2=10.6). 


BOTT-DUFFIN 

CIRCUIT-CASE  A 

(Figure  4.12(a)) 

n*> 

Co/£[pF/m] 

I^/MpH/m] 

C1/2[pF/m] 

R^kD) 

1 

2.9728 

0.4584 

131.99 

2.1290 

3 

0.5063 

0.3875 

1.134 

4.8650 

5 

0.2412 

0.3595 

0.2966 

6.3673 

m 

0.1519 

0.3417 

i 

0.1329 

7.3458 

H 

0.1089 

0.3285 

7 . 529*10~2  | 

8.0824 

^Resonant  frequencies  sn,  n*2,4,...  do  not  contribute  to  the 
driving-point  network. 
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Table  6.3  Element  values  of  the  dominant  pole-pair  admittance 
networks  for  a  quarter-driven  dipole  antenna  with  SI  = 10.6. 


LADDER  CIRCUIT  (Figure  4.7(a)) 

n*> 

C^/£[pF/m] 

Rx[fi] 

L1/«,[yH/m] 

R2tkT2] 

1 

1.6149 

7.307 

0.8243 

4.244 

9 

2.998*10-2 

67.612 

0.4897 

103.08 

BOTT-DUFFIN 

CIRCUIT-CASE  A  (Figure  4.12(a)) 

n 

CQ/£[pF/m] 

L1/4[yH/a] 

Cj/ilpF/m] 

a, [kfi] 

2 

1.0326 

0.4009 

3.1783 

3.2045 

3 

0.3246 

0.8869 

0.2850 

5.0706 

5 

8.738*10~2 

0.5696 

2.4309 

30.882 

6 

0.2171 

0.3375 

0.1780 

5.3000 

m 

0.1209 

0.8301 

3. 797*10-2 

4.5332 

H 

0.1166 

0.3095 

5.693*10-2 

5.7955 

*) 

4th  and  8th  resonant  frequencies  do  not  contribute  to  the 
driving-point  network. 
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Figure  6.5  Input  current  of  a  symmetrical  dipole  antenna  excited  by  a  Gaussian  voltage 
pulse.  SCEPTRE  and  TWTD  results  are  compared. 


tc /i 

.6  Input  current  of  a  quarter-driven  dipole  antenna  excited  by  a  Gaussian  voltage 
SCEPTRE  and  TWIT)  results  are  compared. 


equivalent  circuit  can  be  significantly  improved  by  augmenting  it  by  a 
shunt  capacitor  of  a  proper  value  as  discussed  in  Section  3.3.  The 
value  of  the  capacitor  must  be  computed  apart  from  the  SEM  representation, 
e.g.j  by  solving  for  the  static  charge  distribution  on  the  object.  The 
effectiveness  of  this  procedure  is  illustrated  by  Figure  6.5  for  the 
center-driven  dipole  case. 


The  active  part  of  the  equivalent  network  depends  not  only  on  the 
gap  location  z^,  but  also  on  the  angle  of  incidence  0  of  the  plane  wave. 
As  discussed  in  Chapter  4,  the  voltage  transfer  network  realizing  Tn(s) 
can  be  always  synthesized  in  one  of  the  two  forms:  as  a  ladder  circuit 


or  as  a  symmetrical  lattice  circuit.  For  both  cases  explicit  expressions 


for  the  element  values  in  terms  of  the  poles  sn»  the  admittance  residues 

a^,  the  the  source  coefficients  T^  are  given  in  Chapter  4.  The  expressions 

for  a  and  T  in  the  dipole  case  are  given  in  Sections  6.3  and  6.4, 
n  n 

respectively. 


The  topologies  and  element  values  of  the  first  ten  dominant  pole- 


pair  voltage  transfer  networks  of  a  dipole  (ft =10.6)  for  three  gap 

locations  (z  ■  l/4£,  1/2)1,  and  3/4£)  and  three  angles  of  incidence 
8 

(0-30°,  60°,  and  90°)  are  listed  in  Tables  6.4-6.11.  As  in  the  loop 


case,  we  have  chosen  to  use  the  RC  networks  rather  than  the  equivalent 
RL  realizations  (see  Section  4.4).  It  was  found  that  less  than  10 


dominant  pole-pair  modules  were  enough  to  model  properly  the  dipole 

behavior  for  the  excitations  used.  To  give  the  reader  the  idea  of 

the  complexity  of  a  typical  equivalent  circuit  for  a  cylindrical  dipole, 

we  show  explicitly  in  Figure  6.7  the  topology  of  the  equivalent  circuit 

for  the  case  z  =3/4)1  and  0*60°. 
g 
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Figure  6.7  Topology  of  the  equivalent  circuit  for  a  plane-wave  excited  cylindrical  dipole  antenna 
with  the  port  located  at  z = 3/4H  and  the  angle  of  incidence  0  s1  60°. 


* 


The  short-circuit  current  responses  of  the  equivalent  circuits  for 
a  plane  wave  excited  dipole  (H"10.6)  for  various  angles  of  incidence 
and  various  gap  locations  are  compared  in  Figures  6.8-6.15  to  results 
obtained  by  the  TWTD  code.  The  time  history  of  the  plane  wave  is  a 
Gaussian  pulse.  As  can  be  seen  from  these  figures ,  the  circuit  responses 
computed  by  SCEPTRE  agree  favorably  with  the  responses  obtained  from  the 
time  domain  integral  equation  solution. 

Finally,  in  Figure  6.16  is  shown  the  port  current  of  a  symmetrical 
dipole  antenna  (fi-10.6)  loaded  by  a  parallel  RLC  circuit  with  the 
quality  factor  Q“40,  the  characteristic  resistance  Ro*25ft,  and  the 
resonant  frequency  u)q  ■  it,  excited  by  a  plane  wave  with  the  angle  of 
incidence  0  » 60°  and  a  double  exponential  time  history .  The  solid  line 
represents  the  response  obtained  from  the  numerical  solution  of  the 
integral  equation  (6.2)  (modified  by  the  loading  impedance  in  the  port 
zone)  followed  by  an  FFT  inversion.  The  dashed  line  is  the  equivalent 
circuit  response  obtained  from  the  SCEPTRE  analysis.  It  is  seen  that 
the  agreement  between  these  two  curves  is  satisfactory,  but  not  as  good 
as  in  the  short-circuit  conditions. 

6.6  Simplified  Equivalent  Circuits  for  a  Cylindrical  Dipole  Based 

on  the  Sinusoidal  Mode  Approximation 

It  has  long  been  known  that  the  natural  current  modes  of  a  thin 
cylinder  are  approximately  sinusoidal  [King,  1967).  That  this  is 
indeed  the  case  was  recently  confirmed  by  the  numerical  work  by 
Teeahe  [1973]  where  he  plotted  the  "exact"  natural  modes  of  a  cylinder 
with  Q“10.6.  Hallin  [1930]  used  the  sinusoidal  approximation 
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Figure  6.8  Short-circuit  current  excited  in  the  symmetrical  dipole  antenna  by  a  Gaussian 
function  plane  wave  with  the  angle  of  incidence  9=30°.  SCEPTRE  and  TWTD  results  are 
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Figure  6.9  Short-circuit  current  excited  in  the  symmetrical  dipole  antenna  by  a  Gaussian 
function  plane  wave  with  the  angle  of  incidence  0=60°.  SCEPTRE  and  TWTD  results  are 
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Figure  6.12  Short-circuit  current  excited  in  a  dipole  antenna  with  a  port  located  at 
z = 3/4£  by  a  Gaussian  function  plane  wave  with  the  angle  of  incidence  8=60°.  SCEPTRE 
and  TWTD  results  are  compared. 


Figure  6.15  Short-circuit  current  excited  in  a  dipole  antenna  with  a  port  located  at  z  =  1/42 
by  a  Gaussian  function  plane  wave  with  the  angle  of  incidence  0=60°.  SCEPTRE  and  TWTD  results 
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Figure  6.16  Port  current  of  a  symmetrical  dipole  antenna  loaded  by  a  parallel  RLC  resonant 
circuit  and  excited  by  a  double-exponential  plane  wave  with  the  angle  of  incidence  6=60°. 
The  SCEPTRE  response  is  compared  with  the  response  obtained  by  solving  numerically  the 
integral  eq.  (6.2)  (modified  for  loading)  in  the  frequency  domain  and  FFT  inversion. 


t 


Vz) 


+  oaf1) 


(6.15) 


to  derive  by  an  Iterative  procedure  an  analytical  expression  for  the 
dominant  poles  of  a  thin  cylinder.  In  a  first  order  approximation  he 
obtained 

fn(2nrrD  -Ci(2nn) 

8„ - 5 - a - +  Jn" 


1  - 


Si(2mr) 

rntft 


+  0(Q-2) 


(6.16) 


where  Si  and  Ci  denote  the  sine  and  cosine  integrals  [Abramowitz  and 

Stegun,  1965],  respectively,  F»  1.781...  is  the  Euler's  constant,  and  the 

poles  s^  are  normalized  to  c/£.  Hallen  [1930]  (and  more  recently  Marin 

[1974])  also  carried  the  iteration  one  step  further  and  obtained 
~  -2 

expressions  for  I  (z)  which  are  correct  up  to  ft  and  expressions  for 
n 

_3 

correct  up  to  ft  .  These  expressions  are  rather  complicated, 
however.  It  was  found  by  Langenberg  [1978]  that  the  first  order  formula 
(6.16)  can  be  used  for  the  first  few  resonances  for  cylinders  with 
ft  >12.  The  second  order  formula  is  good  for  ft  >8. 

The  zeroth  order  approximation  to  the  natural  modes  (6.15)  and 
the  first  order  approximation  to  the  poles  (6.16)  were  used  successfully 
by  Marin  and  Liu  [1976]  and  Langenberg  [1978]  for  the  solution  of 
transient  thin-wire  problems  and  a  similar  approach  was  used  by 
Singaraju  and  Gardner  [1976]  for  the  computation  of  a  transient  response 
of  a  helical  antenna. 

By  using  the  so-called  "reduced  kernel"  [King,  1956]  in  (6.2) 
and  the  sinusoidal  current  modes  (6.15)  it  can  be  shown  [Marin,  1974] 
that  the  normalization  constants  (normalized  to  c/i )  can  be  approximated 
by 


179 


(6.17) 


5  ~  ^TT  ~  1 

n  Z  ft  30ft 
o 


with  Zq  =  120tt .  The  admittance  residues  are  thus  given  by 


sin2 


a  £  S  I2(z  )  a 
n  n  n  g 


nir 

t  ZgJ 


30ft 


(6.18) 


For  a  symmetrical  dipole  (z  *  1/2JI)  (6.18)  reduces  to 


an  S  30ft  »  n“1>3*5*-*- 


(6.19) 


Note  that  these  residues  are  real  and  equal  for  all  poles.  For 
*  10.6,  a^*  (3.145  +j0.0)  mS,  which  differs  rather  significantly 
from  the  "exact"  residues  listed  in  Table  6.1. 

Since  the  residues  (6.19)  are  purely  real,  they  fall  into  region 
I  of  Figure  4.5.  Thus,  the  pole-pair  admittances  of  a  symmetrical 
dipole  can  be  realized  in  the  unmodified  form 


Yn(s) 


1 

15ft 


s  +  o 

_ n _ 

s2  +  2o  s  +  |  s  |  2 
n  1  n  1 


(6.20) 


A  continued-fraction  expansion  [ Matthaei ,  1954]  of  (6.20)  leads  to  the 
circuit  in  Figure  6.17(a)  or  to  its  equivalent  shown  in  Figure  6.17(b). 
The  first  of  these  circuits  was  derived  previously  by  Beam  and  Singaraju 
[1980]. 

By  a  slightly  different  procedure  still  another  approximation  to 
the  pole-pair  admittance  can  be  derived: 
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Figure  6.17  Simplified  pole-pair  admittance  circuits  derived  by  the  sinusoidal  mode  approximation 
(frequency  normalization  to  c/2,  and  impedance  level  normalization  to  15f2  is  understood). 


which  can  be  realized  by  the  circuit  shown  in  Figure  6.17(c).  This 
circuit  resembles  the  circuit  used  by  Barnes  [1973],  but  the  expressions 
for  the  element  values  are  different. 

The  input  current  response  of  the  simplified  equivalent  circuit 
for  a  symmetrical  dipole  (£2=10.6)  due  to  a  Gaussian  voltage  pulse  is 
compared  in  Figure  6.18  to  the  result  of  the  TWTD  code  analysis.  Also 
shown  in  this  figure  is  the  response  of  the  "exact"  equivalent 
circuit  using  the  reduced  Bott-Duffin  modules.  As  can  be  seen  from  this 
comparison,  the  response  of  the  simplified  circuit  is  significantly 
worse  than  the  response  of  the  Bott-Duffin  network.  Since  the  complex¬ 
ities  of  the  two  circuit  topologies  are  comparable,  one  can  conclude 
that  it  is  more  expedient  to  employ  the  padding  procedure  and  the 
Bott-Duffin  realization  rather  than  the  sinusoidal  mode  approximation  in 
the  development  of  the  driving-point  equivalent  circuits. 

The  sinusoidal  natural-mode  approximation  can  be  also  used  to 
simplify  the  source  synthesis  problem  of  a  plane  wave  illuminated 
dipole.  Namely,  the  substitution  of  (6.15)  into  (6.12)  gives 


(6.22a) 


for  odd  values  of  n,  and 


o  rr  Q  fS  COBS'* 

n^2  -VV^lcoaS  .i„h  -Ar— 

(-1)  2  f  sinei - ‘ - L_2— L 

sin[nir(5  -%)]  fS  cosS-v 2  ; 

1  '^g  /J  n  ,  mr 

2  J  +  [~J 


(6.22b) 


sin[mr(Cg-%)]  /•s^osS', 
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Figure  6.18  Input  current  of  a  symmetrical  dipole  antenna  excited  by  a  Gaussian  voltage 
pulse.  The  TWTD  response  if  compared  with  the  responses  of  the  equivalent  circuits  using 
the  Bott-Duffin  modules  and  the  simplified  modules  resulting  from  the  sinsoidal  mode 
approximation. 


These  are  still  complex  numbers,  but  a  further  approximation 

S  =  1n7T 

n  J 

leads  to  purely  real  source  coefficients 


(6.24) 


(6.25) 


Since  the  T  ’s  are  real  numbers,  an  examination  of  (3.24)  reveals  that 
n 

the  voltage  transfer  functions  T^Cs)  degenerate  to  real  constants  which 
can  be  realized  by  VCVS  with  gain  factors  given  by  (6.25). 

The  gain  constants  of  the  simplified  dominant  pole-pair  transfer 
circuits  for  a  symmetrical  dipole  (fi-10.6)  are  listed  for  three 
different  angles  of  incidence  6  in  Table  6.12.  The  topology  o5  a 
simplified  equivalent  circuit  using  the  reduced  Bott-Duffin  modules  is 
shown  in  Figure  6.19.  The  short-circuit  current  response  of  this 
network  for  a  plane  wave  with  6  ■  30°  and  a  Gaussian  time  history  is 
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-T  1  <- 


Table  6.12  Gain  constants  Tn'  of  the  voltage  controlled 
voltage  sources  associated  with  the  dominant  pole-pairs  of 
a  symmetrical  dipole  antenna  (ft=10.6)  excited  by  a  plane 
wave  with  the  angle  of  incidence  6*30°,  60°,  and  90°. 


30° 

60° 

90° 

0.26598 

0.51980 

0.63662 

0.25050 

0.17327 

-0.21221 

0.22117 

-0.10396 

0.12732 

0.18103 

-0.07426 

-0.09095 

0.13415 

0.05776 

0.07074 

^Normalized  to  (£E*). 

o 

Ik*) 

Resonant  frequencies  s  ,  n«2,4,...,  do  not  contribute  to 
the  equivalent  network? 


compared  in  Figure  6.20  to  the  response  obtained  from  the  TWTD  code. 

The  agreement  is  seen  to  be  very  satisfactory.  Short-circuit  current 
responses  of  the  network  from  Figure  6.19  due  to  a  double-exponential 
wave  with  three  angles  of  incidence:  0*90°,  60°,  and  30°  are  plotted 
in  Figure  6.21.  These  curves  agree  favorably  with  the  responses 
computed  by  the  WT-MBA/LLL1B  code  [Landt  et  aim,  1974],  included  in 
[Cho  and  Cordaro ,  1979]. 

Finally,  in  Figure  6.22  we  compare  the  response  of  the  circuit 
from  Figure  6.19  with  0  ■  90°  and  a  Gaussian  time  history  to  the  response 
of  the  same  circuit  with  the  Bott-Duffin  modules  replaced  by  the 
simplified  circuits  from  Figure  6.17(a).  It  is  seen  that  there  is  a 
significant  deterioration  of  the  circuit  performance  when  the  sinusoidal 
mode  approximation  is  used  in  the  admittance  modules.  It  appears, 
however,  that  this  approximation  can  be  used  successfully  to  simplify 
the  source  modules  of  an  equivalent  circuit  for  a  symmetrical  dipole 
or  a  cylindrical  post  over  a  ground  plane. 

6.7  Conclusions 

The  responses  of  the  dominant  pole-pair  equivalent  circuits  for 
the  thin-wire  dipole  antenna  agree  favorably  with  the  "exact"  responses 
obtained  from  the  classical  frequency  domain — inverse  FFT  approach. 

The  amount  of  padding  necessary  for  physical  realizability  of  the  pole- 
pair  admittances  is  negligible  in  comparison  with  the  value  of  the  real 
part  of  the  admittances  at  the  resonance. 

It  is  demonstrated  that  the  sinusoidal  mode  approximation  gives 
poor  results  when  used  in  the  synthesis  of  the  dominant  pole-pair 
admittances.  However,  it  can  be  employed  successfully  in  the  synthesis 
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Figure  6.20  Short-circuit  current  excited  In  a  symmetrical  dipole  antenna  by  a  Gaussian 
function  plane  wave  with  the  angle  of  incidence  6=  30°.  The  response  of  the  equivalent 
circuit  with  simplified  transfer  networks  Is  compared  with  the  result  of  the  TWTD  analysis. 


Figure  6.21  Short  circuit  current  responses  excited  In  a  symmetrical  dipole  antenna  by 
double-exponential  plane  wave  with  various  angles  of  incidence. 


Figure  6,22  Short-circuit  current  excited  in  a  symmetrical  dipole  antenna  by  a  Gaussian 
function  plane  wave  and  the  angle  of  incidence  0=90°.  The  response  of  the  equivalent 
circuit  with  the  Bott-Duffin  modules  is  compared  witli  the  response  of  the  circuit  using 
simplified  modules  resulting  from  the  sinusoidal  mode  approximation. 
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of  the  pole-pair  voltage  transfer  networks.  For  a  symmetrical  cylin¬ 
drical  dipole  or  a  cylindrical  post  over  a  ground  plane,  the  pole- 
pair  transfer  networks  reduce  to  just  voltage-controlled  voltage 
sources  with  gain  constants  dependent  on  the  angle  of  incidence  of  the 
impinging  plane  wave. 
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CHAPTER  7 


EQUIVALENT  CIRCUIT  SYNTHESIS  FOR  A  SPHERICAL  ANTENNA 

7.1  Introduction 

Although  the  spherical  antenna  configuration  is  not  one  which 
has  been  employed  frequently  in  practice  (its  main  application  has 
been  in  sensors  for  the  electromagnetic  pulse  (EMP)  measurements 
[Baum,  1969;  Baum  et  at.,  1978]),  it  is  nevertheless  worthwhile  studying 
this  structure  since  it  is  one  of  the  few  finite  shapes  which  can  be 
analyzed  exactly.  Also,  with  its  resonances  so  highly  damped,  it  provides 
a  stringent  test  for  the  dominant  pole  approximation  introduced  in 
Chapter  2. 

In  the  next  sections  SEM  equivalent  circuits  are  derived  for  a 
spherical  dipole  with  a  slot  along  a  longitudinal  line.  The  develop¬ 
ment  presented  here  differs  in  details  from  the  general  theory 
developed  in  Chapter  3  due  to  the  explicit  appearance  of  the  internal 
resonances  of  the  structure  in  the  formulation.  This  phenomenon, 
which  has  not  been  taken  into  consideration  in  the  derivations  of 
Chapter  3,  posed  no  problem  in  the  thin-wire  loop  and  dipole  cases, 
since  the  cavity  resonances  of  these  structures  lie  well  outside  the 
frequency  range  of  interest.  Another  difficulty  not  encountered 
before  is  the  fact  that  the  extent  of  the  port  region  of  the  spherical 
antenna  is  not  electrically  small.  Therefore  the  time  origin  cannot 
be  located  in  the  gap  zone  Independently  of  the  direction  of  arrival 
of  the  incident  field,  as  suggested  in  Chapter  2.  This  difficulty 
is  overcome  by  moving  the  time  origin  to  the  point  of  the  first  contact 
of  the  gap  region  with  the  leading  edge  of  the  impinging  plane  wave. 
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The  static  capacitance  of  the  structure  can  be  considered,  along 
with  the  poles  and  residues,  as  an  important  SEM  parameter.  If  known 


it  can  often  be  used  to  improve  the  early  time  response  of  the  SEM 
equivalent  circuit  by  augmenting  it  by  a  corrective  capacitor  of  a 
proper  value.  Two  methods  of  computation  of  the  static  capacitance 
of  the  spherical  antenna  are  compared  in  Appendix  C. 

The  analysis  presented  in  this  chapter  is  perhaps  more  detailed 
than  the  analyses  of  the  loop  and  cylindrical  dipole  antennas  included 
in  Chapters  5  and  6.  Some  of  the  reasons  for  this  are  that  the 
spherical  geometry  is  amenable  to  an  exact  analysis,  and  that  it 
provides  the  "worst  case"  estimate  of  the  performance  of  the  dominant 
pole-pair-at-a-time  equivalent  circuits. 


7.2  Preliminary  Theory 

The  surface  current  density  J(9,iJ>,s)  on  a  perfectly  conducting 
sphere  of  radius  a  centered  at  the  origin  of  the  spherical  coordinate 
system  (r,0,<(>)  is  given  by  an  integral  equation  of  the  form 

<T(0,<J),e,,(t),,s);J(8',({!’,s)>  *  (T-nn)  •  E^a.Q.^.s)  (7.1) 


with  the  symmetric  product  defined  as 

IT  TT 


<[  1 


iM>5|  I 


[  ]•[  Jazsin9d9d<J>  . 


-tt  0 


(7.2) 


As  in  Chapter  2,  1  denotes  the  identity  dyadic,  n“n(0,4>)  is  a  unit 
vector  normal  at  (0,4>)  to  the  surface  of  the  sphere,  and  E^fa.P.^.s) 
is  the  factorable  incident  field  evaluated  on  the  sphere  surface.  As 
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before,  the  time  history  of  the  incident  field  is  f(t)  with  the  Laplace 
transform  f(s).  The  dyadic  kernel  T  can  be  expanded  in  terms  of  tne 
spherical  wave  functions,  as  shown  by  Tai  [1971].  This  representation, 
when  particularized  to  r*a,  as  above,  constitutes  an  eigenfunction 
expansion  for  F  in  the  eigenfunctions  of  the  integral  operator  in 
(7.1).  From  this  expansion  the  resolvent  kernel  can  be  con¬ 

structed  by  means  of  eigenvalue  reciprocation,  which  leads  to  the 
solution  of  (7.1)  in  the  form 


J(9«<J',s)  —  £ 
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lf(,) 


jj 


(7.3) 

where  the  index  a  takes  on  the  values  "e"  (for  "even")  and  "o" 

(for  "odd")  and  the  superscripts  TM  and  TE  stand  for  "transverse 

magnetic  to  r"  and  "transverse  electric  to  r",  respectively.  The 

meaning  of  the  other  symbols  is  explaned  below. 

~tm  -TE 

The  X^  (s)  and  X^  (s)  are  the  eigenvalues  of  the  integral  operator 
given  by 

-TM 

X  (s)  -  -Z  [s  i  (s)]’[s  k  (s)]’  (7.4) 

n  on  n 

and 

~TF 

X*E,(s)  -  Zo[s  in(s)J[s  kn(s)]  (7.3) 
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1 


where  Zq  is  the  intrinsic  impedance  of  the  medium  and  the  prime  denotes 

differentiation  with  respect  of  the  normalized  complex  frequency 

s**Ya*  Y  being  the  complex  propagation  constant.  The  functions 

i  (s)  and  k  (s)  have  been  introduced  by  Baum  [1971 J  and  are  related 
n  n 

to  the  spherical  Bessel  and  Hankel  functions  as 
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I 

4 


yo  -  r  y-jo 
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(7.6) 

(7.7) 


The  eigenfunctions  Q 


n,m,§ 


(© , 4>)  and  R  e(6,<J>)  are  given  by 
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(7.8) 


<t>  (7.9) 


where  P^  denotes  the  associated  Legendre  functions  of  degree  n  and 
order  m  [Abramowitz  and  Stegun ,  1965],  Note  that  the  eigenfunctions 
are  real  and  frequency  independent.  Also,  there  is  a  manifold  eigen¬ 
value  degeneracy  due  to  the  sphere  symmetry. 

The  normalization  constants  3  are  defined  as 

n,m 
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n,m 


2  2n(n+l)  (iH-m) ! 
(1  +  6  )+a  7-.  P — r4 — 

o  9  m  (2n+l)(n-m)! 


-1 


(7.10) 


where  5  is  the  Kronecker  delta  function  and  the  coupling  coefficients 
are  given  by 


i™}  ■ « 

n,m,o 


5„,.,  <,<«•♦> 


;Eo(a,6,4»)p(a,6,<|),s)> 


(7.11) 
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where  Eq  is  the  spatial  distribution  factor  and  p  is  the  propagation 
factor  of  the  incident  field  ( of .  eq.  (2.6)  of  Chapter  2). 

For  the  purposes  of  equivalent  circuit  synthesis  the  port  current 
I(s)  can  be  defined  as  the  net  surface  current  passing  the  specified 
longitudinal  line  0*9o  which  coincides  with  the  edge  of  a  narrow 
circumferential  slot  in  the  spherical  shell  (Figure  7.1).  This  is  the 
component  of  current  which  would  excite  a  load  attached  to  the  spherical 
antenna  via  a  radial  or  a  biconical  transmission  line.  Thus,  the  port 
current  is  given  by 

TT 

I(s)  *  -  f  §*J(0  ,0,s)sin8  ad$  .  (7.12) 

Jo  o 

-IT 


An  examination  of  (7.3)  reveals  that  only  the  m=0,  cr  -  e  TM  terms 
survive  the  integration  in  (7.12),  so  that 


I(s) 


TT  r  2n+l 

Z  ^  n(n+l) 
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(7.13) 


with  the  coupling  coefficient 


nn(s)  -  <P^(cos0)0;Eo(a, 0,(f)p(a,e,(j>,s)> 


(7.14) 


7.3  Derivation  of  the  Driving-Point  Admittance  of  the  Spherical 
Antenna 

We  consider  a  spherical  shell  of  radius  a  with  a  circumferential 
slot  of  width  d(d/a»2aQ)  with  the  center  located  at  0  *  0^  (Figure  7.1). 
Let  this  structure  be  excited  by  a  ^-independent  electric  field  at  the 
surface  of  the  sphere 

V 

£0(e,s) - f  g(9)f(s)  (7.15) 


where 


and  where 
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|g(0)d6  -  1  (7.16) 
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V(s)  -  VQf (s) 
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(7.17) 


is  the  applied  voltage.  From  (7.13)  and  (7.17)  the  driving-point 
admittance  Y(s)  * I(s)/V(s)  takes  the  form 
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(7.18) 


where  0  *  9  -  a  .  By  using  the  Wronskian  relation  for  the  functions 
o  g  o 

i  (s)  and  k  (s)  [Beam,  1971] 
n  n 
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it  can  shown  that 
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Sin(3) 


skn(s) 


[sin(s)],[skn(s)]'  [sin(s)]'  [skn(s)]' 


(7.20) 


Thus,  the  driving-point  admittance  (7.18)  can  be  divided  into  its 
internal  and  external  parts  [Baum,  1969] 
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and 


198 


-  f  l 


o  n-1 


f  nk  (s) 

H - 2 - v  (0  ,ct  ) 

I  [sk  <3)1 'J  n  8  0 


(7.22) 


where 
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(7.23) 


can  be  considered  as  the  excitation  voltage  of  the  n  eigenmode 
[Franoesohetti ,  1976]. 

The  series  (7.21)  and  (7.22)  assume  the  knowledge  of  the  gap 
field  distribution  g(0).  This  function  depends  on  the  details  of  the 
feeder  goemetry  and  requires  a  solution  of  a  generally  complicated 
boundary  value  problem.  Therefore,  it  is  customary  to  assume  a  priori 
the  field  distribution  across  the  gap.  However,  this  can  be  justified 
only  for  very  narrow  slots  since,  as  demonstrated  by  Franoesohetti 
[1976],  the  admittance  of  the  spherical  antenna  is  rather  sensitive 
to  the  field  distribution  across  the  gap.  Some  of  the  distribution 
functions  used  in  the  literature  are  listed  below. 

1.  Delta  function  [Stratton  and  Chu,  1941;  Karr ,  1951;  Van  Bladelj 
1964;  Ramo  et  at .,  1965] 


g(0)  -  5(0-6  )  ,  0<6<n 

o 

2.  Step  function  [Infold,  1947;  Weeks,  1964] 


(7.24) 


■7  ,  9  -d  <e<0  +d 

d  g  o—  —  g  o 


0  ,  otherwise 


(7.25) 


3.  Distribution  associated  with  a  biconical  transmission  line 
feed  [Schelkunoff,  1952] 
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where  Kq  is  the  characteristic  impedance  of  the  transmission 
line  given  by 
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4.  Distribution  exhibiting  the  edge  singularity  [Baum,  1969; 
Franoesahetti ,  1976] 
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(7.27) 


5.  Distribution  approaching  the  delta  function  as  the  parameter 
m-v  x 


g(8)  -  (sinO)2”*1  ,  O<0<it.  (7.28) 

2^“H‘1(m!)2 

This  distribution  was  used  by  InfeZd  [1947]  to  investigate 
the  convergence  properties  of  (7.22)  as  a  function  of  the 
gap  model.  A  peculiar  feature  of  (7.28)  is  that  it  reduces 
the  number  of  terms  in  the  series  (7.21)  and  (7.22)  to  m. 

The  series  (7.21)  and  (7.22)  converge  very  slowly  (for  the  delta 
gap  model  they  diverge).  In  some  cases  the  nonlinear  transformation 
technique  [Shanks,  1955]  can  be  used  to  accelerate  the  series 
convergence.  This  method  was  successfully  applied  in  a  related  problem 
by  Tesahe  et  at.  ]1976],  It  seems  to  be  more  expedient,  however,  to 
follow  the  technique  employed  by  Sahelkunoff  [1952]  and,  more  recently, 
by  Fvcmcesahetti  [1976]  (see  also  Buaai  and  Franceschetti ,  1974]). 
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They  observed  Chat  the  series  (7.21-7.22)  represents  a  pure  capacitance 
in  the  limit  as  s-*-0.  If  this  series  is  subtracted  term-by-term  from 
the  admittance  series,  its  convergence  properties  improve  dramatically. 


For  the  symmetrical  antenna  case  (0^  ■  90°)  they  were  able  to  sum  the 
capacitance  series  explicitly.  An  application  of  this  technique  to 
(7.21-7.22)  results  in 
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and 
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with  the  unnormalized  interior  and  exterior  capacitances  given  by 


s2k  ^ (s) 
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(7.31) 


and 
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It  is  shown  in  Appendix  C  that  for  a  narrow  equatorial  slot  (aQ  «  1) 
with  a  constant  field  distribution  the  capacitances  C*nt  and  Cext  are 
given  by 


and 


(7.33) 


(7.34) 
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Perhaps  the  only  practical  way  to  drive  or  load  a  spherical  antenna 
is  through  a  biconical  or  a  radial  transmission  line.  In  such 
configuration  only  the  external  surface  of  the  sphere  is  excited. 

•vgXt 

Therefore,  the  following  discussion  will  concentrate  on  Y  (s) 

(the  superscript  "ext"  is  dropped  in  the  rest  of  this  chapter). 

By  defining 


b_  =  V  (6  ,a  ) 

n  g  o 


Z  n 
o 


(7.35) 


(7.22)  and  (7.30)  can  be  rewritten,  respectively,  as 
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(7.38) 


(7.39) 


It  is  of  interest  to  point  out  that  the  normalized  eigenadmittances 
yn(s)  and  y^(s)  have  their  exact  circuit  equivalents.  Namely,  using 
the  recurrence  relations  for  the  spherical  Bessel  function  [Abramcwi tz 
and  Stegun ,  1965]  it  can  be  shown  that  [Cku,  1948;  That,  1978] 
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This  continued- fraction  expansion  can  be  realized  by  the  circuit 
shown  in  Figure  7.2(a).  In  like  fashion,  it  can  be  shown  that 


y;(s) 


-s 


,  2n-l  A 
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2n-3 
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which  corresponds  to  the  circuit  shown  in  Figure  7.2(b).  This  represen¬ 
tation,  which  employs  a  gyrator  with  the  gyration  impedance  j/s, 
was  suggested  by  Fra.ncesch.etti  [1976]. 

The  normalized  eigenadmittance  y  (s)  can  be  explicitly  expressed 

n 

as  a  rational  function  of  s.  This  can  be  done  by  noting  that  [Beam, 
1971] 
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where  wn(l/s)  denotes  a  Bessel  polynomial  in  1/s  [ Balabanian ,  1958]. 
For  notational  simplicity  it  is  helpful  to  define 
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(7.44) 
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m,n 


( 2n-m) ! 

2n  mm! (n-m) ! 


(7.45) 


Note  that  the  numerator  and  the  denominator  polynomials  in  (7.45) 
are  of  the  same  order  and  that  y^(0)  *  0.  It  follows  from  this  that 
the  modified  pole  expansion  (2.24)  is  applicable  in  the  spherical 
antenna  case. 

From  (7.45)  the  zeros  and  poles  of  y^(s)  can  be  easily  found  as 
the  roots  of  the  numerator  and  the  denominator  polynomial,  respectively. 
The  poles  and  zeros  of  the  first  six  eigenadmittances  occuring  in 
the  upper  left  quadrant  of  the  complex  frequency  plane  are  shown  in 
Figure  7.3.  It  should  be  noted  that  the  first  layer  poles  of  the  sphere 
are  much  more  damped  than  the  dominant  poles  of  the  thin-wire  loop 
and  dipole  antennas.  Incidentally,  the  zeros  of  the  TM  admittances, 
except  for  the  zero  at  the  origin,  are  also  the  poles  of  the  terms 
associated  with  the  TE  eigenmodes  in  eq.  (7.3).  Extensive  listings 
of  the  sphere  poles  are  included  in  [Martinez  et  at.,  1972;  Granzow , 
1966]. 

The  residue  r^^  at  the  i^  pole  of  the  nC^  eigenadmittance  s^^ 
can  be  found  as 
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Figure  7.3  Pole  (x)  and  zero  (o)  distributions  of  the  first  six 
TM  eigenmodes  of  the  sphere  (the  zero  at  the  origin  belongs  to  all 
eigenmodes) . 
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This  simple  expression  follows  from  the  Ricatti-Bessel  equation 
[Baum,  1971] 

„2 


sTTn(n+r)fs  kn<s)r'  ~  8  Vs)  *  0  * 


(7.47) 


Denoting  the  dominant  pole  of  the  nC^  eigenadmittance  by  sn  and 

the  corresponding  residue  by  r  ,  the  nC^  normalized  eigenadmittance 

n 


can  be  approximated  by 

yn(s)  s  s 


r  r* 

n  +  n 

s  (s  -  s  )  s*(s  -  s*) 
n  n'  n  n' 


(7.48) 


Using  this  in  (7.36)  and  taking  N  dominant  pole  pairs  gives 

N  N  . 

Y(s)  a  l  [Yn(s)  +  Y_n(s) ]  -  l  Yn(s) 


(7.49) 


n*l 


n»l 


where  the  modified  pole  admittance  Yn(s)  and  the  modified  pole-pair 

A 

admittance  Yn(s)  are  defined  by  (3.11)  and  (3.16),  respectively,  with 
the  admittance  residue  given  by 


b  r 
n  n 


(7.50) 


or,  more  explicitly,  by 


ira  2n+l 


6  +a 

9  c 


■  -  ~r  .  A  -  r  TTY  1  p  (cose  )sin6 

n  Z  d  n(n+l)  [s2+n(n+l)]  n  o  c 


Pn(cos9)sin6d0 


(7.51) 
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In  Figure  7.4  is  shown  the  behavior  of  Y(jui)  of  a  spherical  antenna 

(0  *90°,  d/a  "0.05)  computed  from  the  quickly  convergent  formula  (7.38) 

S 

with  the  static  capacitance  given  by  (7.34)  as  C*  104.74  a  [pF],  It 
is  seen  from  this  figure  that  the  driving-point  admittance  of  the 
spherical  antenna  evaluated  on  the  jui-axis  does  not  exhibit  resonances 
and  remains  capacitive  at  all  frequencies.  On  the  same  figure  is  plotted 
the  admittance  obtained  from  the  dominant -pole  approximation  (7.49). 

Shown  is  also  the  character  of  the  imaginary  part  when  the  dominant 
pole-pair  admittance  is  augmented  by  a  corrective  shunt  capacitor 
adjusting  the  static  capacitance  of  the  circuit  to  the  proper  value. 

It  is  evident  from  the  comparison  of  thase  plots  that  the  dominant- 
pole  formula,  even  with  the  capacitive  correction,  gives  a  poor 
approximation  to  the  admittance  of  the  spherical  antenna. 

7.4  Derivation  of  the  Short-Circuit  Current  of  the  Spherical  Antenna 
For  the  short-circuit  current  computation  the  spherical  antenna, 
with  its  slot  closed,  is  excited  by  a  plane  wave  incident  field  with 
the  polarization  as  indicated  in  Figure  7.1.  No  generality  is 
sacrificed  by  limiting  the  discussion  to  this  polarization  since  the 
field  component  parallel  to  the  x,y  plane  does  not  contribute  to  the 
net  current  across  the  gap.  Also,  the  discussion  can  be  simplified 
by  making  the  propagation  vector  p  parallel  to  the  x,z  plane.  Since 
the  antenna  geometry  is  independent  of  $,  this  represents  no  loss  in 
generality,  either. 

The  plane  wave  incident  field  is  given  by  eq.  (2.7)  of  Chapter  2 
with  Eq * u  Eq,  where  G  is  the  unit  polarization  vector,  as  shown  in 
Figure  7.1.  Since  the  extent  of  the  gap  region  of  the  spherical  antenna 
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Figure  7,4  The  spectrum  of  the  driving-point  admittance  of  a  spherical  antenna  with  an  equatorial 
slot.  The  dominant  pole  approximation  (7.49)  is  compared  with  the  exact  expression  (7.38). 


cannot  be  considered  electrically  small,  one  cannot  choose  the  time 
origin  location  tq  independently  of  the  direction  of  arrival  of  the 
incident  field.  This  difficulty  is  overcome  by  moving  the  time  origin 
to  the  point  of  the  first  contact  of  the  incoming  wave  with  the  gap 
region.  Thus,  with 


r  *  a  p  cos(9  +  6, ) 
o  o  1 


(7.52) 


and 


P  '  * 


—  cos(0  +  6. ) 
c  o  1 


(7.53) 


the  incident  field  becomes 


-  -  -Y  P*  (r~ro) 

Eo(r’3)  ‘  Eoe 


-ST 


G  e 
o 


u  e 


f(s) 
-Y  P  *  r 


(7.54) 


where  s  *  ya  is  the  normalized  complex  frequency  and  tq  ■  tQc/a.  The  last 
term  in  (7.54)  can  be  expanded  in  terms  of  the  spherical  vector  wave 
functions  as  follows  [Beam,  1971;  Morse  and  Feshbaoh,  1953] 

-*■  a>  n 


u  e 
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-Y  P 


w  n  r~  — I 

•  f  r*  n  '►  -►  *+• 

■  7  T  a  M  (yr)  +  b  N  (yr) 
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n-i.  m*u  *—  — 1 


(7.55) 


Mn  n,  rr(Yr)  "  i(Yr)R  _  _(©»4>) 

n,m,o  n  n>nijC7 


(7.56) 
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[(yr)i  (yr)]’ 

<9’*>+ - (W) - Qn,m,o(W)  (7-57> 


with  $  (8 ,<J>)  and  X  (9,4>)  defined  in  (7.8)  and  (7.9)  and 

ii}iQy{j  n,m,  o 


p 

n,m, 


e(0,0) 


e 


Pm(cos9) 

n 


fcos(m<j>) 

\sin(m4>) 


(7.58) 
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The  coefficients  in  (7.55)  are  given  by  [Baum,  1969] 


and 


n_K  nn+l  2n+l  (n-m)  l  Pn(c0s9l) 

n,m  o,nr  n(n+l)  (n+m)!  m  sin0^ 
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n,m  ^  o,mM  ;  n(n+l)  (n+m)! 


90, 


(7.59) 


(7.60) 


With  the  help  of  (7.54)  and  (7.55)  the  coupling  coefficient 
(normalized  to  EQa)  can  be  found  from  (7.14)  as 


[s  i  (s)]'  -ST 

fi  (s)  ■  b  -  e 

n  n,o  s 


TTn: 


9  P  (cos0) 
n 


IF 


sin0d8 


[s  i  (s)  ] '  -sr  . 

2(_i)n  - 2 -  e  0  pl(cos0i)  . 


(7.61) 


When  (7.61)  is  substituted  for  the  coupling  coefficient  in  (7.13),  the 

term  [s  i^s)]',  which  is  associated  with  the  purely  imaginary  cavity 

resonances,  is  cancelled  out.  The  short-circuit  current,  normalized 

to  E  a,  takes  the  form 
o 


-ST 


ISC(s)  «  l  p-  (-l)n  -~TT  P1  (cos0  )sin0  P1(cos0  )  -  e-  rr 
n-1  Zo  n(n+l)  n  o  on  1  s[s  kn(s)J 


7  f(s). 


(7.62) 

To  expand  (7.62)  in  terms  of  its  singularities  it  is  necessary  to  examine 
the  term 

-ST  -ST 

o  o 

e _  _ _ e _ _ 

s[s  kn(s)]'  “  “  n  s  kn(s)  +  s^kn_1(s) 


npn(s)  +  s2pn_i(s) 


(7.63) 
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where  the  polynomial  p  (s)  is  defined  in  (7.43).  The  comparison  of 


(7.63)  with  (7.45)  shows  that  ISC(s)  has  the  same  poles  s  .  as  Y(s). 

ni 

The  residues  of  (7.63)  at  these  poles  can  be  found  as 
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(s  -  S  .) 


ni  s[s  ^n(s)  I ' 


-s  ,T 
ni  o 


Snifs  kn(sni)]” 


"S  4T 

ni  o 


n+l  sni(1 ~  V 
Sni  6 
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(7.64) 


In  the  dominant  pole  approximation  the  short-circuit  current 
becomes 
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(7.65) 


which  can  be  rewritten  as 


I  <s)  ~  l  Y  (s)T  f(s) 
n=l  n 


(7.66) 


with  the  source  coefficients  T  given  by 

SnTo  P^(cos0.) 
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Combining  N  conjugate  pole-pairs  in  (7.66)  results  in 

N 

ISC(s)  ^  l  [Yn(s)Tn  +  Y_n(s)T*]f(s) 
n=l 

N  ^ 

=  l  Y  (s)T  (s)?(s)  (7.68) 

n  n 

n=l 

A 

with  the  pole-pair  voltage  transfer  functions  T^(s)  defined  in  (3.24). 

In  Figure  7.5  the  spectrum  of  the  short-circuit  current  of  the 

spherical  antenna  (0  =90°,  d/a -0.05)  computed  from  (7.62)  with  the 

8 

delta  function  time  history  (f(s)  -  1)  is  compared  with  the  spectrum 
obtained  from  the  dominant  pole  approximation  (7.68).  The  agreement 
between  these  two  spectra,  although  quite  poor,  is  nevertheless  much  more 
satisfactory  than  the  agreement  in  analogous  admittance  results 
(Figure  7.4). 

7.5  Equivalent  Circuits  for  the  Spherical  Antenna  With  an  Equatorial 
Slot  of  Width  d/a  =  0.05 

In  Table  7.1  are  listed  the  admittance  residues  a  associated  with 

n 

the  first  20  dominant  poles  of  a  spherical  antenna  with  an  equatorial 

slot  (0  -90°)  of  width  d/a  =  0.05.  Included  also  for  each  pole  s  are 
8  n 

the  quality  factor  and  the  amount  of  the  real  part  padding  Gn 

required  for  the  realizability  of  the  associated  pole-pair  admittance. 

The  column  denoted  gives  the  ratio  (in  %)  of  the  padding  conductance 

G  to  the  resonance  value  of  the  real  part  of  Y  (jo>).  It  is  seen  from 
n  n 

this  table  that,  except  for  the  first  few  poles,  the  amount  of  padding 

G  is  a  significant  fraction  of  the  value  of  the  real  part  at  the 
n 

resonance.  Only  the  first  pole  pair  admittance  does  not  require 
padding  and  is  realizable  in  the  modified  form.  This  is  further 
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Figure  7.5  The  spectrum  of  the  short-circuit  current  excited  in  a  spherical  antenna  with  an 
equatorial  slot  by  a  broad-side-incident  delta  function  plane  wave.  The  dominant  pole  approximation 
(7.68)  is  compared  with  the  exact  expression  (7.62). 
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illustrated  in  Figures  7.6-7.11  where  the  behavior  of  the  real  part  of 
some  of  the  dominant  pole  pair  admittances  is  plotted  as  a  function  of 
the  real  frequency  to.  It  is  seen  that,  except  for  the  first  pole-pair 
admittance,  there  is  a  significant  negative  excursion  of  the  real  part 
on  the  jto-axis. 

The  element  values  and  the  circuit  topology  of  the  dominant 

pole-pair  admittance  networks  for  the  spherical  antenna  with  0^  «  90° 

and  d/a *0.05  are  included  in  Table  7.2.  The  transient  response  of  this 

equivalent  network  due  to  a  Gaussian  voltage  pulse  is  compared  in 

Figure  7 . 12  with  the  "exact"  response  obtained  from  (7.38)  and  an 

FFT  inversion.  Shown  also  is  the  response  of  the  equivalent  circuit 

augmented  by  a  corrective  capacitor.  In  both  cases  the  discrepancy 

between  the  circuit  response  and  the  exact  response  is  significant. 

In  Table  7.3  are  defined  the  dominant  pole-pair  voltage  transfer 

networks  of  the  spherical  antenna  (9  =90°,  d/a  =  0.05)  for  the  case 

8 

of  a  broad-side  plane  wave  incident  (0^=90°). 

The  short-circuit  current  response  of  the  complete  equivalent 
network  for  the  spherical  antenna  due  to  a  plane  wave  incident  field 
(0^  =  90°)  carrying  a  step  function  time  history  is  compared  in  Figure 
7.13withthe  "exact"  response  obtained  from  (7.67)  by  an  FFT  inversion. 
Again,  the  agreement  between  these  curves  is  not  satisfactory. 

7.6  Conclusions 

The  spherical  antenna,  with  its  resonances  so  highly  damped, 
provides  an  extremely  stringent  test  for  the  dominnnt-pole-at-at-time 
equivalent  circuits.  The  dominant-pole  approximation,  which  works 
so  well  for  the  thin  wire  loop  and  dipole  antennas  (Chapters  5  and 
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Figure  7.11 
of  the  spherical  antenna  with 


Table  7.2  Element  values  of  the  dominant  pole-pair  admittance 
networks  for  a  spherical  antenna  with  an  equatorial  slot  (0  *  90°) 
of  width  d/a  =  0.05.  ® 


LADDER  CIRCUIT  (Figure  4.7(a)) 


C^/afpF/m] 

R^n] 

L^/a[yH/m] 

r2[«] 

41.66 

0.0 

266.71 

80.0 

BOTT-DUFFIN  CIRCUIT-CASE 

(Figure  4.12(a)) 

CQ/a[pF/m] 

L^/a[yH/m] 

C^/a[pF/m] 

R1tfi] 

11.923 

132.45 

33.456 

115.84 

7.269 

93.48 

11.884 

123.54 

5.348 

74.257 

6.416 

125.38 

4.282 

62.613 

4.136 

125.34 

3.597 

54.737 

2.941 

124.61 

3.115 

49.027 

2.225 

123.66 

2.755 

44.69 

1.756 

122.67 

2.474 

41.287 

1.430 

121.75 

2.246 

38.547 

1.191 

120.96 

in 

cvi 


o 

c\i 


in 


in 

6 


n 


i 


Figure  7.12  Input  current  of  the  spherical  antenna  with  an  equatorial  slot  excited  by  a  Gaussian 
voltage  pulse.  The  responses  of  the  dominant  pole  equivalent  circuit  with  and  without  the 
corrective  capacitor  are  compared  with  the  exact  response. 
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Figure  7.13  Short-circuit  current  excited  in  a  spherical  antenna  with  an  equatorial  slot  by 
broad-side-incident  step  function  plane  wave.  Tl  e  response  of  the  dominant  pole  equivalent 
circuit  is  compared  with  the  exact  response. 


6),  is  a  rather  crude  one  in  the  spherical  antenna  case.  Also,  the 
dominant  pole-pair  admittances  of  the  spherical  antenna  require  a 
much  larger  amount  of  real-part  padding  (relative  to  the  value  of  the 
real  part  of  the  admittance  at  the  resonance)  to  render  them  positive 

real  than  do  the  loop  and  dipole  antennas.  As  a  consequence  of  this, 
the  performance  of  the  dominant  pole-pair  equivalent  circuits  for 

the  spherical  antenna  is  considerably  poorer  than  the  performance  of 
similar  circuits  for  the  thin-wire  loop  and  dipole  structures,  as 
demonstrated  in  Chapters  5  and  6.  It  seems  that  most  of  the  inaccuracy 
in  the  circuit  response  is  introduced  by  the  driving-point  admittance 
part  of  the  equivalent  circuit.  The  performance  of  the  source  part 
is  significantly  better  than  that  of  the  admittance  part  of  the  equivalent 
network,  so  that  the  short-circuit  current  waveform  does  not  differ 
significantly  from  the  true  response. 

Due  to  its  extreme  low-Q  character,  the  spherical  antenna  provides 
a  worst  case  estimate  of  the  performance  of  the  dominant-pole  approximation 
in  the  sense  that  for  poles  residing  so  far  from  the  ju)-axis  the  dominance 
of  the  poles  nearest  the  axis  is,  at  best,  weak.  Although  the  dominant 
pole  results  are  relatively  poor  in  Figures  7.12  and  7.13,  they  agree 
within  35  percent  with  more  accurately  computed  results.  We  therefore 
conclude  that  in  the  low-Q  limit,  the  approximation  fails  gracefully 
and  that  for  structures  whose  dominant  pole  quality  factors  are  inter¬ 
mediate  between  those  of  wire  structures  and  the  sphere,  the  dominant 
pole-approximation  may  be  sensibly  applied. 
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CHAPTER  8 


SUMMARY  AND  CONCLUSIONS 

Practical  and  systematic  procedures  have  been  developed  for  the 
synthesis  from  the  SEM  (Singularity  Expansion  Method)  description  of 
acitve,  physically  realizable  equivalent  circuits  for  electromagnetic 
energy  collecting  structures  with  identifiable  ports.  The  necessary 
SEM  information  includes  the  poles  {s^}  of  the  structure  and  the 
associated  natural  modes  and  normalization  coefficients.  From  the  poles 
and  the  associated  modal  distributions  one  can  compute  for  a  given 
incident  field  the  coupling  coefficients  which  weight  a  given  singularity's 
contribution  to  the  total  response.  Two  different  coupling  coefficient 
forms  have  been  used  in  the  SEM  work  to  date — the  so-called  class  1 
and  class  2  forms  [Baum,  1976b].  In  the  SEM  equivalent  circuit  synthesis 
only  the  class  1  form  is  applicable,  which  is  known  to  exhibit  numerical 
instability  in  early  times  [Marin,  1972;  Pearson,  1979].  However,  it 
is  demonstrated  in  Chapter  2  that  the  simple  class  1  coupling  coefficient 
can  be  used  with  success  provided  the  time  origin  of  the  problem  is 
placed  in  the  port  zone  (gap  region) ;  it  is  also  argued  that  the  class 
1  expansion  is  always  convergent  as  Refs^]  -*--<»  for  this  time  origin 
location. 

The  SEM  expansion  for  the  current  density  J(r,s)  on  a  scattering 
object  is  derived  formally  in  Chapter  2  from  the  eigenfunction  expansion 
of  J(r,s)  in  terms  of  the  eigenvalues  and  eigenvectors  of  the  electric 
field  integral  equation  (EFIE)  operator.  It  is  assumed  in  this 
development  that  a  complete  set  of  eigenvectors  exists  which  is  not 
necessarily  true  since  the  EFIE  operator  is  non-selfadjoint.  In  fact, 
it  is  demonstrated  in  Appendix  A  that  the  eigenvalues  of  the  EFIE 
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operator  may  have  branch  points  in  the  complex  frequency  plane  at 
which  two  or  more  eigenvalues  become  degenerate  so  that  the  inclusion 
of  root  vectors  in  the  expansion  may  be  necessary  [Ream,  1980],  If 
an  eigenvalue  has  a  branch  point,  a  branch  cut  must  be  chosen  in  the 
s-plane,  resulting  in  a  branch-integral  term  in  its  SEM  expansion. 

This  can  have  an  important  bearing  on  the  SEM  equivalent  circuit 
synthesis  because  in  most  practical  cases  these  branch-integral 
constituents  are  not  explicity  identifiable.  Their  omission,  however, 
can  result  in  non-real izable  terminal  eigenadmittances  [Pearson  and 
Wilton,  1981]. 

In  the  process  of  expanding  the  inverse  eigenvalue  in  terms  of  its 
singularities  one  refers  to  the  Mittag-Lef f ler  theorem  from  the  theory 
of  functions  of  a  complex  variable.  Two  basic  SEM  expansion  forms 
result  depending  on  the  large-s  asymptotic  behavior  of  the  eigenvalue — 
the  so-called  modified  and  unmodified  forms  [Beam,  1976a],  It  is 
observed  in  Chapter  2  that  the  modified  form  imposes  less  stringent 
conditions  on  the  asymptotic  behavior  of  the  eigenvalues  of  the  integral 
operator.  Therefore  this  form  appears  to  be  preferable  in  the  case  of 
problems  solved  numerically.  This  suggestion  is  supported  by  the  fact 
that,  as  shown  in  Chapter  7,  the  modified  SEM  expansion  is  the  proper 
one  in  the  case  of  the  sphere,  which  Is  one  of  the  few  structures 
for  which  an  exact  analytical  solution  is  available. 

The  SOI  equivalent  network  can  be  derived  either  from  the  Norton 
or  from  the  Thevenin  equivalence,  leading  to  two  different  but  equivalent 
networks.  Of  course,  the  SEM  parameters  which  must  be  supplied  are 
different  in  those  two  cases.  The  short-circuit  SEM  data  are  usually 
found  by  solving  the  short-circuit  boundary  value  problem,  whereas  the 
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open-circuit  parameters  result  from  the  solution  of  the  open-circuit 
boundary  value  problem.  The  open-circuit  SEM  parameters  are  much  more 
difficult  to  obtain  than  the  short-circuit  parameters  of  the  same 
object  because  of  the  highly  localized  charge  interaction  effects 
between  the  two  parts  of  the  structure  in  the  vicinity  of  the  gap 
which  must  be  modeled  precisely  in  the  open-circuit  problem.  Further¬ 
more,  it  is  enough  to  solve  only  one  short-circuit  boundary  value 
problem  for  all  possible  gap  locations  on  the  structure,  whereas  the 
open-circuit  problem  must  be  solved  for  each  location.  Thus,  if  for 
some  reason  the  Thevenin-based  equivalent  circuit  is  preferred  over 
the  Norton-based  equivalent,  it  is  perhaps  more  expedient  to  derive 
it  from  the  short-circuit  parameters,  as  described  in  Chapter  3. 

The  SEM  equivalent  circuit  synthesis  procedures  developed  in  this 
work  are  based  on  the  first-layer,  dominant  pole-pairs  because  they 
dominate  the  response  of  the  structure  and  they  are  the  ones  easiest 
to  extract  by  numerical  or  experimental  means.  An  alternative 
approach  of  grouping  the  poles  according  to  eignvalues  is  not 

practical,  however,  because  in  most  cases  this  collection 
of  poles  is  not  explicitly  identifiable.  Also,  the  number  of  poles 

grows  rapidly  for  higher  order  modes,  making  the  corresponding  circuits 
extremely  complicated  [Streable  and  Pearson ,  1981]. 

The  SEM  equivalent  circuits  developed  in  Chapter  4  have  the  form 
of  a  connection  of  simple  modules  associated  with  the  conjugate  pole 
pairs  of  the  structure.  These  pole-pair  circuits  are  synthesized  one 
at  a  time,  so  if  the  need  arises  the  network  can  be  easily  expanded 
by  adding  more  pole-pair  modules  without  repeating  the  whole  synthesis 
process.  By  using  the  network  topologies  developed  in  Chapter  4  the 
synthesis  process  is  reduced  to  the  realization  of  a  driving-point 
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function  biquadratic  in  s  and  a  voltage  transfer  function  bilinear 
in  s.  The  realizability  conditions  and  explicit  circuit  realizations 
are  given  for  these  functions. 

If  the  driving-point  function  is  realizable,  it  can  be  synthesized 
by  simple  continued  fraction  expansion.  If  a  pole-pair  driving- 
point  function  is  initially  nonrealizable,  a  simple  padding  procedure 
is  described  to  synthesize  an  approximating  circuit.  The  padded  pole- 
pair  driving-point  function  is  a  minimum-real-part  function.  It  appears 
that  the  only  transformerless  realization  applicable  in  this  case  is 
the  Bott-Duffin  network.  If  the  amount  of  padding  is  negligible,  this 
network  can  be  simplified  by  deleting  one  leg  without  introducing 
appreciable  error. 

It  is  shown  that  the  pole-pair  voltage  transfer  function  can  be 
always  realized  either  in  the  form  of  a  ladder  circuit  or  in  the  form 
of  a  symmetrical  lattice  circuit.  Unfortunately,  the  applicability  of 
each  of  the  two  topologies  is  dependent  on  the  direction  of  arrival 
of  the  incident  plane  wave,  so  that  the  circuit  topology  may  change  when 
excitation  conditions  are  changed.  In  each  of  the  two  cases,  two 
equivalent  realizations  are  possible:  an  RC  realization  and  an  RL 
realization.  In  the  cases  tested  it  was  found  that  the  element  values 
of  the  RC  netv,  rk  were  more  realistic  than  the  element  values  of  the 
equivalent  RL  realization. 

Each  of  the  transfer  networks  must  be  followed  by  a  voltage- 
controlled  voltage  source  which  serves  as  an  impedance  buffering 
stage  and  as  an  amplifier  to  rest., -re  the  gain  lost  in  the  associated 
transfer  network.  The  controlled  sources  can  be  implemented  in  practice 
with  active  devices  which,  however,  usually  cannot  sustain  high  voltages. 
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The  synthesis  procedures  developed  are  tested  on  thin-wire  loop 
and  dipole  antennas  (Chapters  5  and  6),  as  well  as  on  a  spherical 
antenna  (Chapter  7).  The  synthesized  equivalent  circuits  are  analyzed 
by  a  general  circuit  analysis  program  (SCEPTRE)  and  the  responses 
compared  with  classical  frequency  domain  and  numerical  inverse  Fourier 
transform  results  or  with  results  obtained  from  the  Thin  Wire  Time 
Domain  (TWTD)  code.  In  the  case  of  the  thin-wire  structures,  the  results 
are  very  encouraging.  The  circuit  responses  agree  favorably  with  the 
responses  obtained  by  other  means.  Of  course,  the  quality  of  the 
circuit  responses  is  commensurate  with  the  quality  of  the  dominant- 
pole  approximation,  which  is  very  accurate  in  the  case  of  highly  resonant 
(slender)  structures.  The  only  additional  approximation  in  the  SEM 
circuit  is  the  possible  padding  introduced  in  the  synthesis  of  the 
pole-pair  driving-point  function.  This  padding  is  negligible  in  the 
thin-wire  loop  and  dipole  cases;  it  is  rather  significant  for  the 
spherical  antenna.  In  the  case  of  the  spherical  antenna  the  dominant 
pole  approximation  breaks  down,  too.  For  these  reasons  the  SEM  equivalent 
circuit  of  the  spherical  antenna  cannot  be  expected  to  perform  well, 
unless  modules  associated  with  the  second  and,  perhaps,  third  layer 

of  poles  are  included  in  the  circuit. 

The  driving-point  responses  of  the  equivalent  circuits  are  consider¬ 
ably  poorer  than  the  short-circuit  responses  excited  by  a  plane  wave 
incident  field.  This  can  be  attributed  mainly  to  the  capacit've 
deficiency  of  the  SEM  equivalent  circuit.  Namely,  the  capacitance  of 
the  finite  equivalent  circuit  in  the  limit  as  s -*■  0  underestimates  the 
true  static  capacitance  of  the  structure.  Therefore,  the  early- 
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time  response  of  an  SEM  equivalent  circuit  can  often  be  improved  by 
augmenting  it  by  a  shunt  capacitor  of  a  proper  value.  This  phenomenon 
does  not  influence  the  short-circuit  response  because  the  input  terminals 
of  the  structure  are  shorted  for  that  case.  Most  of  the  comparisons 
included  in  this  work  involve  the  short-circuit  conditions  primarily 
because  only  for  those  cases  were  benchmark  results  readily  available. 

In  view  of  the  above,  these  tests  may  appear  to  be,  perhaps',  too  for¬ 
giving.  However,  the  loaded  dipole  example  included  in  Chapter  6 
indicates  that  our  equivalent  circuits  also  perform  well  under  finite 
loading  conditions.  The  SEM  equivalent  circuits  with  nonlinear  loads 
are  analyzed  in  Appendix  B.  No  data  were  available  for  comparison 
in  those  cases,  but  the  results  lend  themselves  to  a  straightforward 
physical  interpretation. 

In  the  thin-wire  dipole  case  the  sinusoidal  natural-mode  approxi¬ 
mation  can  be  employed  to  simplify  the  SEM  equivalent  circuits.  It 
is  demonstrated  in  Chapter  6  that  this  approximation  is  poor  in  the  case 
of  the  driving-point  circuits  but  it  can  be  applied  successfully  to 
the  voltage  transfer  networks.  For  a  symmetrical  dipole  antenna  or 
a  cylindrical  post  over  a  ground  plane,  particularly  simple  and  yet 
reasonably  accurate  circuits  result. 

In  conclusion,  the  main  objective  of  this  work  has  been  fulfilled. 
The  SEM  equivalent  circuits  developed  are  physically  realizable,  simple, 
and  reasonably  accurate.  These  circuits  possess  most,  but  not  all  of 
the  desirable  features  listed  in  the  Introduction.  In  particular,  the 
network  topology  does  change  with  the  change  of  the  direction  of 
arrival  of  the  incident  field.  However,  the  equivalent  circuit  does 
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not  depend  on  the  time  history  of  the  incident  field.  Also,  the  autonomic 
voltage  source  is  connected  to  only  one  port,  and  the  equivalent  circuit 
is  transformerless,  but  it  does  include  controlled  sources. 
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APPENDIX  A 


ON  THE  EXISTENCE  OF  BRANCH-POINTS  IN  THE  EIGENVALUES 
OF  THE  ELECTRIC  FIELD  INTEGRAL  EQUATION  (EFIE)  OPERATOR 
IN  THE  COMPLEX  FREQUENCY  PLANE 


A.  1  Introduction 

The  Singularity  Expansion  Method  (SEM)  [Baum,  1971]  representation 
for  the  current  density  on  a  scattering  object  can  be  derived  formally 
from  the  eigenfunction  expansion  involving  eigenvalues  and  eigenfunctions 
of  the  associated  integral  operator  [Baum,  1975].  It  has  generally 
been  assumed  that  the  inverse  eigenvalues  are  analytic  in  the  complex 
frequency  plane  (s-plane)  except  for  pole  singularities  [Pearson  and 
Wilton,  1981].  This  is  suggested  by  Marin  and  Latham's  result 
[Marin  and  Latham,  1972]  that  in  the  case  of  finite  extent  perfectly 
conducting  bodies  in  lossless  media,  poles  are  the  only  singularities 
in  the  SEM  expansion.  ^  However,  as  pointed  out  by  Pearson  and  Wilton 
[1981],  this  does  not  preclude  the  possibility  of  the  occurrence  of 
other  singularities  which  cancel  in  the  complete  sum  representing 
the  singularity  expansion  of  the  current.  As  it  turns  out,  in  fact, 
only  for  bodies  such  as  the  sphere  or  the  circular  loop  antenna,  in 
which  geometrical  symmetry  completely  determines  the  eigenfunctions 
(and  hence  they  are  frequency  independent),  has  it  been  possible  to 
show  (Section  2.3)  that  such  singularities  cannot  occur. 


^In  contrast  to  finite  extent  bodies,  infinite  cylindrical  objects 
always  manifest  solutions  with  a  branch  point  singularity  of  well- 
understood  origin  at  zero  frequency.  Schafer  and  Kouyoumjian  [1975] 
discuss  this  branch  point  for  the  case  of  a  conducting  circular  cylinder. 
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The  presence  or  absence  of  branch  points  has  an  important  bearing  on 
the  validity  of  the  Eigenmode  Expansion  Method  (EEM)  [Baum,  1975] 
representation  of  the  resolvent  kernel,  because  at  a  branch-point  in 
the  s-plane  two  or  more  eigenvalues  become  degenerate  and  the  inclusion 
of  root  vectors  [Bamm,  1980]  in  the  expansion  may  be  necessary  at  that 
point.  The  expansion  of  individual  reciprocal  eigenvalues  or  any 
associated  "eigen-quantities"  to  obtain  their  singularity  expansion 
can  also  require  integrals  along  the  concomitant  branch  cuts  in  order  for 
the  expansion  to  be  complete.  The  previously-mentioned  meromorphicity 
result  of  Marin  and  Latham  for  finite  extent  objects  ensures  that  when 
the  eigenvalue  expansions  are  summed,  the  branch-integral  constituents 
cancel  and  the  traditionally-used  SEM  forms  result.  However,  in  some 
applications,  singularity  expansions  are  applied  on  an  eigenmode-at- 
a-time  basis  and  the  branch  integrals  must  be  retained.  For  example, 
the  eigenvalues  and  their  reciprocals  manifest  a  "generalized  positive¬ 
realness"  [Pearson  and  Wilton,  1981]  which  is  important  in  SEM  circuit 
synthesis.  To  discard  the  branch-integral  constituent  of  their 
expansion  may  compromise  this  property. 

Although  the  existence  of  branch-integral  constituents  in  the  SEM 
representation  has  been  speculated  upon  [Baum,  1978;  Pearson  and 
Wilton,  1981],  no  specific  examples  of  their  occurrence  have  been  cited, 
to  date  (except  for  the  well-known  infinite  cylinder  case,  mentioned 
above).  The  purpose  of  this  paper  is  to  fill  this  gap. 

In  Section  A.  2  we  demonstrate  simple  circuit  and  transmission  line 
examples  in  which  the  eigenvalues  of  the  impedance  matrix  have  branch¬ 
points.  Although  only  finite-dimensional  impedance  operators  are 
involved  in  these  examples,  they  provide  valuable  insight  into  the 
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of  the  branch-point  occurrence.  In  particular,  they  clearly  demonstrate 
the  role  of  the  object  symmetry  in  the  occurrence  of  branch-points  in 
the  eigenfunction  expansion  of  the  inverse  operator.  In  Section  A. 3 
we  consider  the  problem  of  scalar  scattering  by  a  prolate  spheroid 
and  demonstrate  the  presence  of  branch  points  in  the  eigenvalues  of 
the  associated  integral  operator  which  follows  from  the  analysis  of 
the  spheroidal  wave  equation  in  the  complex  frequency  plane  [Oguchi, 

1970].  In  Section  A. 4,  similar  analysis  of  the  Mathieu's  differential 
equation  [Blanch  and  Clemm ,  1969]  leads  us  to  the  conclusion  that 
branch  points  besides  that  associated  with  the  infinite  extent  of  the 
object  exist  in  the  frequency  plane  behavior  of  the  eigenvalues  of 
the  EFIE  for  an  infinite  perfectly-conducting  elliptic  cylinder. 

A.  2  Branch  Points  in  Circuit  and  Transmission  Line  Problems 
A.  2.1  General  Remarks 

The  terminal  quantities  of  an  N-port  linear,  passive,  and  bilateral 
electrical  network  are  related  via  the  equation 

(Z][I)  «  [V]  ,  (A . 1) 

where  fl]  and  [V]  are  column  vectors  containing  N-port  Laplace-transformed 

currents  and  voltages,  respectively,  and  [Z]  is  the  N*N  symmetric 

open-circuit  impedance  matrix.  Eq.  (A.l)  is  analogous  to  the  electric 

field  integral  equation  (2.1).  The  primary  difference  is  that  the  electro¬ 
magnetic  problem  consists  of  an  uncountably  infinite  number  of  ports. 

Even  this  distinction  vanishes,  however,  when  a  numerical  solution  of 

(2.1)  is  undertaken  because  then  a  matrix  approximation  replaces  T  and 

column  vectors  analogous  to  voltage  and  current  replace  and  J, 

respectively  [Hai'X'ington,  1968].  In  the  following,  some  conjectures 
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concerning  (2.1)  are  made  based  on  observations  of  properties  of  (A.l). 

In  general,  an  eigenvector  expansion  of  the  inverse  of  [ZJ  can 
be  found  as  follows: 


r  -ip  -i  t 

M  I  (s)  I  (s) 

p  ,  -i  N  n  n 

[z‘  (s)]  ■  Jx  f  ~w  1 

I  (a)  I  (s) 
n  In 


(A.  2) 


where  the  superscript  T  denotes  transpose  and  X^  and  [I  ]  are 
eigenvalues  and  eigenvectors  of  [Z].  In  the  regions  in  the  s-plane 
where  (A. 2)  holds,  solution  of  (A.l)  can  be  found  in  a  form  analogous  to 
(2.13).  However,  it  is  demonstrated  in  simple  examples  below  that  there 
may  be  points  in  the  s-plane  where  the  eigenvalues  become  degenerate  and 
the  simple  expansion  (A. 2)  is  no  longer  valid. 

Before  presenting  specific  examples  it  is  helpful  to  investigate 
what  general  conclusions  can  be  drawn  concerning  the  eigenvalues  of 
a  two-port.  The  characteristic  equation  of  a  lumped  reciprocal  two- 
port  network  described  in  the  form  of  (A.l)  is 


det 


Z11~X  Z12 


12 


z22-X 


X‘  -  (i-11  +  ::22)X  +  ZllZ22_Z12  =  ° 


(A. 3) 


where  an<*  z22  are  raC:'-ona^  functions  in  s.  From  (A. 2. 3) 

eigenvalues  become 

X ^  (s)  -  +  z22  t  ^ ^ s »  (A. 4 ) 

where  the  discriminant  A(s)  is  given  by 
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(A. 5) 


The  radical  in  (A. 4)  suggests  that,  in  general,  the  eigenvalues  have 
branch  points  in  the  s-plane.  Branch  points  do  not  occur,  however,  if 
one  of  the  following  conditions  is  satisfied: 

1.  (uncoupled  ports),  which  is  a  trivial  case; 

2.  A(s)  =0,  which  can  happen  only  if  (zji  -  z22-*  =  —  j 2z^ 2 »  h°wever, 
this  situation  is  not  possible  since  zjj  are  rational  functions 
with  real  coefficients; 

3.  A(s)  is  frequency  independent  or  of  the  form  [P(s)/Q(s)]2 
where  P(s)  and  Q(s)  are  polynomials  in  s; 

4.  zn  ”  z22’  f°r  a  symmetric  network. 

The  role  of  the  network  symmetry  in  eliminating  branch  points  in 
the  impedance  matrix  can  also  be  illustrated  in  the  example  of  a  three- 
port  wich  rotational  symmetry.  The  impedance  matrix  of  such  a  network 


is  given  by 


(A.  6) 


The  characteristic  equation 

det  ([Z]  -  Ml])  ■  +  2z^2  “  A)  (X  +  z^2  -  zn^2  (A. 7) 

yields  the  eigenvalues 

^i  »  z ^  +  2z12  (A. 8a) 

and 

A2  *  z^  -  z12  (A. 8b) 

which  have  no  branch  points  if  z^  and  z^2  are  rational,  but  where 
is  of  algebraic  multiplicity  2.  The  eigenvector  corresponding  to  is 


237 


1 


[Ill 


and  the  eigenvectors  corresponding  to  X ^  are 


1 

1 

^.l1  - 

ej  211/3 
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ej47T/3 

ej2TT/3 

(A. 9a) 


(A. 9b) 


Since  two  independent  eigenvectors  corresponding  to  X^  could  be  found, 
its  geometric  and  algebraic  multiplicities  are  equal.  Also,  owing  to 
the  symmetry  of  the  network,  we  see  that  the  eigenvectors  are  the 
so-called  "symmetrical  components"  and  are  frequency  independent. 

The  impedance  matrix  of  the  last  network  example  has  properties 
similar  to  the  EFIE  operator  in  the  case  of  the  circular  loop  antenna 
and  the  spherical  antenna — both  rotationally  symmetric  bodies.  As  in 
the  case  given  above,  both  the  circular  loop  and  the  spherical  antenna 
have  multiple  eigenvalues  for  which  independent  eigenfunctions  can  be 
found.  These  eigenfunctions  are  frequency  independent,  too. 

A. 2.2  Two-Port  Example — Lumped-Element  Circuit 

The  simple  two-port  of  Fig .  A.  1  (a)  provides  an  example  explicitly 
illustrating  the  occurrence  of  branch  points  in  the  eigenvalues  of  an 
impedance  matrix  and  their  mutual  cancellation  in  the  singularity  expansion  of 
the  inverse  of  the  matrix.  The  impedance  matrix  of  the  circuit  is 
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of  the  impedance  matrix. 


m  = 


l  +  s/2  1 


1  +  2s 


(A. 10) 


with  the  eigenvalues 


X  (s)  *  1  +7  s  ±  /A(s)  , 

1  4 

2 


where 


A(s)  =  1 + 


3 

4  S 


(A. 11) 


(A. 12) 


and  with  the  eigenvectors 


M 


7-  s  i  /A(s) 
4 


(A. 13) 


,-l. 


From  (A. 11)  one  finds  that  X^(— 5/2)  =0  and  X^fO)51©,  and  hence  [Z  ] 
has  poies  at  s  =  0  and  s  = -5/2,  However,  there  are  also  branch  points 
at  s  ■  tj4/3,  where  the  eigenvalues  become  degenerate  (Fig.  A. 1(b)).  Away 
from  these  isolated  points  we  can  find  [Z  according  to  (A. 2)  as 


Z_1(s)1 


7-  s  +  /A  (  s ) 
4 


\  S  +  /Ms)  jj-  s  +  /A(s)J 


-1  1  + 


4  s  +  v/A(s_)  Zt'KisJ  pA(s)  +  ~s 

4  _ 


i 

i 

) 

4 


!  7s”  /A(s) 

4 

|  s  -  /A(s)  |j  s  -  /A(s)J 


1  +  |s  -  /A(s)  2/A(s)  ,STiT-|s 


(A. 14) 
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At  the  branch  points  the  denominators  of  (A. 14)  vanish  and  the  expansion 
is  not  valid.  However,  if  the  terms  in  (A. 14)  are  combined  together  we 
obtain 


M  ■ 


2s  +  1 


-1 


s2  +  -|  s  s2  +-|  s 


-1 


s/2  +  1 


s2  +-j  S  S2  +-J  s 


(A. 15) 


Its  only  singularities 


Obviously,  there  are  no  branch  points  in 
are  the  poles  at  s  =  0  and  at  s  =  -5/2. 

At  a  branch  point  the  root  system  can  be  used  to  expand  the  inverse 
matrix.  The  coefficients  of  this  expansion  can  be  found  witli  the  help 
of  the  biorthogonal  root  system  [Pease,  1965].  For  example,  at  the  branch 
point  s=+j4/3,  [Z]  has  a  double  eigenvalue  X  =  l+j5/3  and  a  root  system 


[I1,1]  =  [j]  ‘  [I1,2]  =  [l] 


(A. 16) 


If  the  scalar  product 


<[I1]5[I21>  -  [Ij  [i2J 


(A. 17) 


is  introduced  where  t  denotes  the  Hermitian  transpose,  a  biorthonormal 
root  system  can  be  found  as 


[vi-r:i-  iv] 


(A.  18) 


.-l, 


Using  (A.  16)  and  (A. 18)  [Z  ]  can  be  expanded  at  the  branch  point 
s  ■  +J4/3  as 
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7  "if  4  4  fl  _  1 
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r-H 
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1 
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i _ 

i:  f  5-1 

1 

»— t 

o 

fl+ll’ 

1 " 

L-i  d 

(A. 19) 


with  a  similar  conjugate  representation  at  the  other  branch  point, 
s  =  -j4/3.  Note  that  the  eigenfunction  expansion  (A. 14)  is  valid  through¬ 
out  the  complex  s-plane  except  at  the  two  isolated  branch  points. 

One  might  speculate  that  the  resistive  element  in  this  circuit 
is  the  feature  which  leads  to  the  branch  points.  Since  we  wish  to  draw 
analogies  with  scattering  from  lossless  surfaces,  this  would  indeed 
be  untimely.  However,  one  can  readily  construct  LC  circuit  examples 
as  well,  where  the  branch  points  arise.  Furthermore,  examples  are  easily 
constructed  in  which  the  branch  points  lie  in  the  right  half  plane. 


A. 2. 3  Two-Port  Example — Non-Uniform  Transmission  Line 

To  illustrate  that  the  branch  point  occurrence  in  the  eigenvalues 
is  not  limited  to  lumped  systems,  we  consider  a  section  of  lossless 
non-uniform  transmission  line  of  length  L  (Fig.  A. 2(a)),  The  impedance 
and  admittance  per  unit  length  are  assumed  to  vary  exponentially  with  the 
distance  x  along  the  line  according  to,  respectively,  z(x)  =  exp(2otx)  and 
y (x)  =  y  exp(-2ax),  where  z  =  si  and  y  =  sc  and  i,  and  c  denote,  respectively, 
the  per-unit-length  inductance  and  capacitance  of  the  line.  Viewing  the 
line  as  a  two-port  network  we  may  find  its  chain-parameter  matrix  by  the 
method  described  in  ( BertnolTi ,  1967 J.  From  the  chain-parameter 
description  the  impedance  matrix  can  be  found  as 


fZ] 


—  ctnh(l'L)  +  — 

y  y 


y  sinh(l’L) 


y  sinh(fL) 


(A.  20) 
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Figure  A. 2  (a)  Exponential  transmission  line,  (b)  branch  points  of  the 

eigenvalues  of  the  impedance  matrix  and  (c)  branch-point  trajectory  as 
a  function  of  the  taper  parameter. 


i 

i 


24  3 


where  P  =  /a2  +  zy  =  /a2  +  (s/v)  with  v=l / /lc  denoting  the  speed  of  light 
(a  deviation  from  the  generally  accepted  notation  is  made  here  because 
the  symbol  c  is  used  in  this  section  to  denote  the  per-unit-length 
capacitance  jf  the  line).  The  inverse  of  [Z]  can  be  readily  found  as 


rz-f 


Lz 


-  -a  L 

£  e _ 

z  sinh(ri) 


-a  L 


-  ctnh(TL)-  -  -  - 


z  sinh(fi) 


ctnh(rL) +- 


-2a  L 


(A. 21) 


If  the  hyperbolic  functions  in  I A. 21 )  are  expanded  in  a  power  series. 


z'1  , 


hence 


it  is  seen  that  only  even  powers  of  T  enter  the  elements  of 
the  inverse  matrix  has  no  branch  points.  Its  only  singularities  are  the 
infinite  number  of  poles  given  by 


~  =  ±j/(mr)2  +  (a L)2  ,  n-0,1,2,... 


(A. 22) 


The  eigenvalues  of  [Z]  can  be  found  to  be 


X  l  =  e 

7. 


a  L 


r 


^  ctnh(FL)cosh(aL)  -  ^  sinh(aL)  ±  /A(s) l  , 


(A. 23) 


where 

A(s) 


(— )  ctnh2 (TL)cos2 (al)  -  (^)  (y)ctnh(rL)sinh(2aL) 


+  (— )  sinh2 (ctL)  -  —  . 

y  y 


(A. 24  ) 


Obviously,  there  are  branch  points  in  the  eigenvalues  (A, 23)  which  can 
he  I  omul  from  the  equation  A(s)  =0.  A  numerical  investigation  of 
(A, 24)  yields  four  branch  points  located  symmetrically  in  the  four 
quadrants  in  the  s-plane  (Fig.  A, 2(b)).  Note  that  the  branch  points 
are  not  restricted  to  the  left  half-plane.  The  trajectory  of  the 
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f irst-quadrant  branch  point  for  different  values  of  the  parameter  a 
is  shown  in  Fig.  A. 2(c). 

We  recover  a  uniform  transmission  line — i%e.,  a  symmetric  structure — 
by  letting  cx -*•().  In  the  limit,  the  branch  points  recede  to  infinity, 
continuing  the  trend  indicated  in  Fig.  A. 2(c). 

A. 3  Scalar  Scattering  from  a  Prolate  Spheroid 

The  role  of  symmetry  in  eliminating  the  appearance  of  branch  points 
in  the  preceding  examples  suggests  that  in  the  scattering  problem  we 
should  examine  a  finite-extent  structure  which,  for  simplicity,  conforms 
to  a  separable  coordinate  system  but  which  is  not  symmetric  under  a 

general  rotation,  something  separable,  but  not  spherical. 

The  vector  Helmholtz  equation  is  not  separable  in  spheroidal  coordinates, 
but  the  scalar  Helmholtz  equation  is.  We  therefore  consider  the 
associated  scalar  scattering  problem. 

The  geometry  of  the  structure  under  consideration  and  the  coordinate 
system  are  shown  in  Fig.  A.  3.  The  surface  S  of  the  spheroid  is  defined 
by  the  radial  coordinate  u  =  Uq,  and  the  semifocal  distance  is  5 .  We 
are  interested  in  the  eigenvalues  of  the  integral  equation  in  C'  defined 
by 

G0(r,r',s)^(r’ ,s)dS'  =  F(r,s),  reS  ,  (A. 25) 

S 

-►  -+ 

where  Cg(r,r  ,s)  is  the  scalar  Green  function  defined  in  (2.4)  and 
K(r,s)  is  the  forcing  function. 

Expanding  the  Green  function  in  (A. 25)  on  the  surface  S  in  terms  of 
the  spheroidal  wave  functions  gives  \Flarmer,  1957] 
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•»  <*»  (  2  —  <S  ) 

Y  r  r  om 

2tt  i  N 

m*0  n=m  mn 


S 

mn 


(q»v)smn(q»v' ) 

mn 


x  cos[m(4>-4>')]R^  (q,u  n)R^(q,u  ),  r.r'e  S 
mn  u  mn  u 


(A. 26) 


where  6  is  the  Kronecker  delta,  q  ■  -jy£  is  the  frequency  parameter, 
and  Smn  and  R^^(i*l,4)  satisfy  the  following  differential  equations 
which  arise  in  the  separation  of  the  Helmholtz  equation  in  prolate 
spheroidal  coordinates  [ Flammev,  1957]: 
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d_ 

dv 


S 

mn 


(q)  -  q 
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S  (q.v)  =  0 
mn 


(A. 27) 


and 


d_ 

du 


f  R(1)(q 

du  mn 


,u)  -  a 


a(q)  -  q‘ 

mn 


m2 

U2  -  L 


R(^(q.u) 

mn 


0  , 


(A. 28) 


where  m  and  n  are  integers  and  a  are  the  separation  parameters.  The 

mn 

normalization  factors  appearing  in  (A. 26)  are  defined  through  the 
orthogonality  of  the  angle  (periodic)  functions 
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(q»v>dv 

mn  mn 

-i 
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nn 


,N 

mn 


and  are  found  to  be 


mn 


1) 


2  I 

r=0 


(r +  2m)! 


[<%>J 


(2r  +  2m  +  l)r! 


(A. 29) 


(A. 30) 


where  the  functions  d™°(q)  are  the  coefficients  in  the  expansion  of 

(he  angle  f  unci  Ions  S  In  terms  ul  Lite  assoc  J  a  toil  l.ogoinl  re  functions: 

mn 


^The  prime  over  the  summation  sign  indicates  that  the  summation  is 
over  (only)  even  values  when  (n-m)  is  even  and  over  (only)  odd  values 
when  (n-m)  is  odd. 
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(A. 31 ) 


1 

; 

1 

3* 

4 


( 


Smn(q*v)  =  I  dr"(q)C-r(v)  * 
r-0 

The  substitution  of  this  expansion  into  (A. 27)  yields  a  three-term 

recursion  formula  for  the  coefficients  dmn.  From  this  recursion 

r 

formula,  we  obtain  a  transcendental  equation,  the  roots  of  which  are 
the  separation  parameters  a  f Oauehi,  1970].  This  equation  has  the  form 


u(a»q)  =  0  . 

mn 


(A. 32) 


It  follows  from  (4.8)  that  the  separation  parameters  a  coupling 

inn 

(A. 27)  and  (A. 28)  are  functions  of  the  frequency  parameter  q. 

The  orthogonality  given  by  (A. 29)  is  awkward  in  identifying  (A. 26) 
with  an  eigenfunction  expansion  of  the  (scalar  counterpart  to)  form 
(2.11),  because  the  weight  in  the  integration  of  (A. 29)  is  unity,  while 
in  (A.  23)  the  surface  metric  is 


dS'  =  8.2/(u2  -  l)(u2  -  v'2)  dv  *  d  (f> 1  . 

O  O 


(A. 33) 


However  if  we  interpret 


'Kv.tf.s)  =  £2/(u2  -  1 )  (u2  -  v  ’  2 )  iKv,i>,s) 


(A. 34) 


as  Che  unknown  in  (A. 25)  and  define 

2  it  1 


<r 


n>  dv  d$  , 


(A. 35) 


0  -1 

then  our  symmetric  product  yields  the  orthogonality  integral  of  (A.2h). 
By  comparison  with  (2.11)  we  find  that 


X  (q)a  R(1)(q,u  )R(4)(q,u  ) 
mn  4/  mn  o  mn  o' 


(A. 36) 
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in  the  eigenvalues  (A. 36) 


The  radial  (non-periodic)  functions 

mn 

can  be  computed  in  terms  of  the  spherical  Bessel  functions  [Fliimmrr', 

1957].  Here,  we  are  interested  only  in  their  analytical  properties  in 

the  complex  q-plane.  Since  satisfy  the  differential  equation 

mn 

(A. 28),  their  behavior  in  the  q-plane^  can  be  inferred  from  the 

properties  of  the  associated  separation  parameters  a  .  The  branch 

mn 

points  of  a  (q)  are  also  branch  points  of  and  thus  of  * 

mn  mn  mn 

The  functions  a^Cq)  have  been  investigated  in  the  literature  [Oguchi, 

1970;  Meixner  et  at.,  1980]  and,  indeed,  it  has  been  found  that  they 

have  branch  points  in  the  q-plane.  The  branch  point  locations  can  be 

computed  by  the  procedure  outlines  below. 

The  behavior  of  the  separation  parameters  a  in  the  q-plane  is 

mn 

governed  by  the  differential  equation  [Wilton  et  at.,  1975] 


h  3mn(q)  =  -  / 
a 


(A. j7) 


.suj-ei.  results  from  (A.  32)  (the  subscript  denotes  a  partial  di t  ( erent i;«l 
with  respect  to  the  variable  indicated).  Integration  of  (A. 37)  on 
closed  contour  C  in  the  q-plane  yields 


A  a 


mn 


(A. 38) 


This  integral  gives  a  non-zero  contribution,  i.e. ,  a  is  multi-valued, 

mn 

if  the  contour  C  encloses  a  Zero  of  U  .  The  location  of  the  branch 

a 

point  of  a  coincides  with  the  location  of  the  zero  of  U  .  Thus,  the 
•  mn  a 

branch  points  can  be  found  by  solving  (A. 32)  together  with 


Following  the  convention  accepted  in  the  literature  we  use  the 
q-plane  rather  than  the  s-plane  here.  The  normalized  frequency  plane 
is  related  to  the  q-plane  by  the  trivial  mapping  Yt-jq. 
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(A. 39) 


u,(a  »q) 
a  mn 


These  simultaneous  transcendental  equations  can  be  solved  by  a  Newton- 
like  iteration  [Oguchi,  1970], 

A  partial  listing  of  branch  points  in  the  eigenvalues  for  a 
number  of  indices  m,n  is  given  in  [Oguahi,  1970].  A  more  complete 
list  is  included  in  [Meixner  et  at.,  1930].  As  an  example,  we  display 
in  Table  A.l  a  number  of  branch  points  for  m=0-2  and  n-0-4.  V»e  observe 
that  some  of  the  branch  points  for  different  eigenvalues  coincide, 
which  means  that  an  eigenvalue  for  some  m,n  can  be  obtained  t>y  an 
analytic  continuation  into  another  Riemann  sheet  in  the  q-plane  of  an 
eigenvalue  with  different  indices.  This  observation  substantiates 
the  degeneracy  interpretation  of  the  branch  point  which  is  mentioned 
in  Section  2.3. 

When  the  spheroid  is  deformed  into  a  sphere,  which  can  oe  accomplishe  ' 
by  letting  the  focal  distance  l  go  to  zero,  the  points  in  the  y-plano 
corresponding  to  branch  points  in  the  q-pl ane  recede  to  infinity— 
a  phenomenon  analogous  to  that  observed  in  the  case  of  a  non-uniform 
transmission  line (Section  A. 2. 3).  Indeed,  it  is  well  known  that  the  eigen¬ 
values  of  the  sphere  are  proportional  to  the  spherical  Hankel  functions 
of  the  second  kind,  which  are  meromorphic  in  the  finite  complex 
frequency  plane.  The  absence  of  branch  points  can  be  thus  attributed 
to  the  geometrical  symmetry  of  this  structure,  which  is  lacking  in  the 
spheroid  case. 
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Table  A.l  Branch  points  of  the  eigenvalues  in  the  prolate  spheroid 
problem. 


m  n 

Branch  Points  of  a  (Y&) 

m,n 

0  0,2 

-2.601670  +  j  1.824770 

0  2,4 

-5.807965  +  j  2.094267 

0  2,4 

-3.081362  +  j  5.217093 

0  1,3 

-2.887165  +  j  3.563644 

1  1,3 

-4.097453  +  j  1.998555 

1  2,4 

-4.492300  +  j  3.862833 

2  2,4 

-5.449457  +  j  2.736987 

r  i 


A. 4  rM  Scattering  from  an  Elliptic  Cylinder 

In  this  section  we  consider  an  electromagnetic  scattering  problem 
which  demonstrates  the  introduction  of  branch  points  because  of  its 
departure  from  rotational  symmetry:  that  of  transverse-magnetic  (TM) 
scattering  from  a  perfectly  conducting  elliptic  cylinder.  Since  the 
structure  is  infinite  in  extent,  a  branch  point  at  the  origin  of  the 
frequency  plane  appears  as  in  the  case  of  the  circular  cylinder 
[Schafer  and  Kouyoumjian ,  1975].  It  is  demonstrated  here  that  the 
elliptic  cross-section  introduces  additional  branch  points  into  the 
frequency  plane. 

The  geometry  of  the  structure  and  the  coordinate  system  used 

are  shown  in  Fig.  A. 4.  The  elliptical  surface  of  the  cylinder,  which  is 

infinite  in  the  z-direction,  is  defined  by  the  radial  coordinate  u=uq. 

The  semifocal  distance  is  l.  The  integral  equation  for  the  current 

densitv  J  (p,y)  excited  by  the  incident  electric  field  E  (i  , <)  is 
2  2 


-  f  Z  y 
4  o 


J  (p '  ,Y)H^2)  (-jy|p-p'  |  )dL'  =  E*(p,y),  DEL  , 
z  u  z 


(A. 40) 


where  'is  the  Hankel  function  of  the  second  kind  and  zeroth  order 
and  is  the  two-dimensional  free-space  Green  function. 

The  Hankel  function  can  be  expanded  in  terms  of  Mathieu  functions 
as  [Jones,  1964] 


H 


(2) 


(- JY  I P-P '  i  )  “  2  l  ce  (v,q)ce  (v’,q)Mc^1\uo,q)Mc^4)(u(i,q) 


m=0 


+  2  l  sem(v,q)sem(v' ,q)Msm  (uQ,q)Msm  (»Q.q)  . 
m*l 

p,p'  EL  .  (A. 41) 
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vs  const. 

u=  const. 


The  frequency  parameter  is  q*-(r£/2)2  and  ce  and  se  are,  respectively, 

m  m 

the  even  and  odd  angular  (periodic)  Mathieu  functions  [ McLaahlan ,  1947] 
satisfying 

I> +M-2,..82vilr^:uo 
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ffam(q)l 

+ 

and  and  Ms^^(i 

m  m 

functions  satisfying 


and  McU;  and  Ms'  ;(i  =  l,4)  are  the  radial  (non-periodic)  Mathieu 
m  m 


i!_  + 

dU2 

Mc^ (u , q ) 
m 


=  0  . 


(A. 43) 


As  in  the  spheroidal  equation  case  in  the  preceding  section,  the 
separation  constants  a^(q)  and  bm(q)  of  the  Mathieu' s  equation  satisfy 
transcendental  characteristic  equations  and  are  dependent  on  the  complex 
variable  q.  The  angular  functions  possess  the  following  orthogonality 
properties: 

2tt 


ce  (v,q)ce  (v,q)dv  *  tt  6  , 

m  n  mn 


se  (v,q)se  (v,q)dv  =  n 
J  m  n 


Sem(v,q)cen(v’q)dv  ”  0  » 


(A. 44a) 


(A»44b) 


(A. 44c) 


while  the  integration  metric  is 


dL  m  £/cosh2u  -  cos2v  dv 
o 


(A. 45) 
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By  applying  reasoning  analogous  to  that  used  in  the  spheroid  development, 
we  conclude  that  the  frequency  variation  of  the  eigenvalues 
of  the  operator  in  (A. 40)  follows 

[Kcn1>(Vq)Mcn4)(Vq)] 


X  (qW 

“o 


Ms^(u  ,q)Ms^(u  ,q) 
no  no 


(A. 46) 


Since  the  radial  functions  Mc^  and  Ms are  related  to  the  Hankel 

n  n 

functions  of  the  second  kind  [MaLachlan,  1947],  they  share  with  them 

the  branch  point  at  the  origin — the  branch  point  due  to  the  infinite 

extent  of  the  scatterer.  The  radial  functions  in  (A. 46)  manifest 

additional  branch  points,  however. 

As  in  the  spheroid  case,  the  analytic  properties  of  these  functions 

in  the  complex  q-plane  can  be  inferred  from  the  properties  of  the 

separation  parameters  a^(q)  an<*  b^(q).  These  functions  have  been 

investigated  in  the  literature  [Blanch  and  Clemm,  1969;  Meixner  et  at., 

1980]  where  it  was  shown  that  a  and  b  have  branch  points  in  the  q- 

n  n 

plane,  and  the  locations  of  these  points  were  computed  for  a  number 

of  values  of  n  by  a  scheme  similar  to  that  outlined  in  Section  A. 3. 

The  listing  of  branch  points  in  a  and  b  in  yt-plane  for  n=0-7  is 

n  n 

displayed  in  Table  A. 2. 

When  the  elliptic  cylinder  is  deformed  into  a  cylinder  with 
circular  cross-section,  the  eigenvalues  of  the  external  scattering 
problem  become  the  Hankel  functions  having  only  a  branch  point  at  the 
origin.  Thus  wo  may  associate  the  additional  branch  points  in  the 
elliptic  cylinder  with  the  lack  of  symmetry  in  this  structure. 


255 


4 


Table  A-2  Branch  points  of  the  eigenvalues  in  the  elliptic  cylinder 
problem. 


- - - - - - 1 

n 

Branch  Points  of  a^Cv?)  ' 

0,2 

-1.713919  +  j  1.713919  ! 

1,3 

-1.917584  +  j  3.376788  I 

2,4 

-2.046414  +  j  4.988464  1 

3,5 

-2.140312  +  j  6.582698  i 

4,6 

-5.739541  +  j  5.739541 

4,6 

-2.213980  +  j  8.168790  j 

5,7 

-5.895408  +  j  7.313526 

5,7 

-2.274478  +  j  9.750397 

n 


Branch  Points  of  b  (>?,) 
n 


2, A 


3,5 


4,6 


5,7 


6,8 


6,8 


7,9 


-3.722621  +  j  3.722621 
-3.968A60  +  j  5.416204 
-4.159649  +  j  7.064678 
-4.316064  +  j  8.689667 
-7.758450  +  j  7.758450 
-4.448352  +  j  10.300680 
-8.027379  +  j  9.468520 
-4.562911  +  j  11.902630 


7,9 


A. 5  Conclusions 


It  is  demonstrated  by  several  specific  examples  that  eigenvalues 
of  impedance-type  integral  operators  can  have  branch  points  in  the 
complex  frequency  plane  (s-plane).  Through  a  two-port  circuit  example, 
it  is  demonstrated  that  as  s  approaches  a  branch  point  two  eigenvalues 
become  degenerate,  the  eigenfunction  expansion  of  the  inverse  operator 
becomes  invalid,  and  the  root  vector  expansion  must  be  used  to  expand 
the  inverse  operator  at  this  point.  If  the  terms  involving  the  degenerate 
eigenvalues  are  grouped  together  a  representation  is  obtained  which 
does  not  exhibit  any  branch  points  and  simple  poles  are  the  only 
singularities  of  the  inverse  operator. 

In  the  eigenanalysis  of  scattering  from  a  prolate  spheroid  or 
elliptic  cylinder  branch  points  also  occur  and  two  or  more  eigenvalues 
become  degenerate.  In  fact,  since  some  of  the  branch  points  for 
different  eigenvalues  coincide,  one  of  the  eigenvalues  can  be  obtained 
by  the  analytic  continuation  into  another  Riemann  sheet  in  the  s-plane 
of  the  "neighboring"  eigenvalue. 

These  examples  demonstrate  that  symmetry  plays  an  important  role 
in  the  presence  or  absence  of  branch  points  in  the  eigenvalues  of 
the  EFIE  operator.  It  is  conjectured  that  branch  points  may  always 
be  present  when  sufficient  object  symmetry  is  lacking.  This  conjecture 
is  supported  by  the  fact  that  branch  points  appear  when  the  sphere  is 
deformed  into  a  spheroid  or  when  a  circular  cylinder  is  deformed  into 
an  elliptic  one.  The  analogous  phenomenon  has  been  observed  in  the 
circuit  and  transmission  line  problems.  For  example,  when  the  taper 
parameter  of  a  non-uniform  transmission  line  goes  to  zero  (the  line 
becomes  uniform),  the  branch  points  move  away  to  infinity. 


The  examples  substantiate  that  the  previously-speculated-upon 
branch  points  do  occur  in  the  scattering  theory.  In  addition,  they 
suggest  a  tie  between  the  branch  points  and  the  need  for  root  vectors 
in  a  basis  spanning  the  solution  space.  Whether  the  need  for  root 
vectors  in  the  solution  expansion  is  isolated  to  branch  points  in  the 
frequency  plane  and  what  ties  exist  between  the  branch  points  and 
the  normal  operator  condition  stated  by  Ream  [1980]  remain  open 
questions. 
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APPENDIX  B 


AN  APPLICATION  OF  THE  SEM  EQUIVALENT  CIRCUITS  TO 
THE  ANALYSIS  OF  ANTENNAS  WITH  NONLINEAR  LOADS 

In  this  appendix  we  give  some  examples  of  the  application  of  the 
SEM  equivalent  circuits  to  the  analysis  of  antennas  with  nonlinear 
loads.  Nonlinear  effects  are  important  for  antenna  systems  containing 
semiconductors,  integrated  circuits,  and  voltage  limiters  when  they 
may  be  illuminated  by  an  extremely  strong  signal,  such  as  that  produced 
by  a  lightning  strike  or  an  electromagnetic  pulse  (EMP)  [Liu  and 
Tesahe ,  1976]. 

One  approach  to  the  analysis  of  nonlinearly  loaded  energy-collecting 
structures  has  been  to  solve  a  space-time  domain  integral  equation  for 
the  current  flft lief  et  at.,  1976;  Liu  and  Tesahe ,  1976;  Schuman,  1974], 

At  each  time  step  a  system  of  nonlinear  algebraic  equations  must  be 
solved  in  this  method  by  a  Newton-like  iterative  procedure.  Even 
complex-shaped  bodies  could  be  treated  by  this  technique  if  the  procedures 
developed  recently  by  Rao  et  al.  [1981]  are  applied.  This  approach 
is  not  very  practical,  however,  since  the  complex  boundary  value 
problem  must  be  resolved  each  time  the  loading  is  changed. 

An  original  technique  for  the  analysis  of  antennas  with  nonlinear 
loads  has  been  developed  recently  by  Tesahe  and  Liu  [1975]  (see  also 
[Liu  and  Tesahe ,  1976 \Liu  et  al.,  1977].  They  use  frequency  domain 
data  of  an  unloaded  structure  and  a  numerical  Fourier  transform  to 
formulate  a  Volterra-type  integral  equation  for  the  current  in  the 
time  domain.  Subsequently,  this  equation  is  solved  by  stepping  in 
time  and  iteration  procedure. 


It  is  believed  that  the  approach  proposed  in  the  present  work, 
i,e.,  the  application  of  the  SEM-derived  equivalent  circuits  and  a 
standard  circuit  analysis  code  to  the  analysis  of  nonlinearly  loaded 


energy-collecting  structures,  can  be  an  useful  alternative  to  the 
method  of  Tesahe  and  Liu  [1975],  It  has  the  advantage  that  no  special 
purpose  computer  code  is  required.  Furthermore,  the  input  languages 
of  many  advanced  circuit  analysis  programs,  such  as  SCEPTRE  [Jensen  and 
McNamee,  1976],  allow  ou<-  to  define  almost  any  real  or  conceivable 
loading  situation.  Also,  an  EMP  with  an  arbitrary  time  history  can 
he  easily  defined. 

In  the  examples  given  below  the  nonlinear  loading  is  provided 
by  an  1N414B  high-speed  switching  diode.  The  diode  was  represented 
in  the  circuit  by  a  voltage  dependent  current  source  I = Io[exp(AV)-l ] 
with  1^  -  2.9*10  ^  A  and  A = 15  V  *  in  shunt  with  a  0.5  pF  capacitor 
modeling  the  junction  capacitance.  The  diode  parameters  were  taken 
from  [Liu  et  al 1977]. 

In  Figure  B.l  is  plotted  the  input  current  response  of  a  circular 
loop  (H  *  15)  driven  by  a  double  exponential  function  voltage  source 
(with  a  peak  voltage  of  10V)  through  a  diode  connected  in  the  forward 
direction.  Shown  is  also  the  input  current  with  the  diode  removed. 

The  topology  and  element  values  of  the  equivalent  circuit  used  are 
given  in  Table  5.2. 

In  Figures  B.2  and  B.3  are  shown  the  current  and  voltage  responses 
at  a  port  of  a  plane-wave-excited  loop  antenna  (I,!=15)  loaded  by  a  diode. 
The  incident  wave  time  history  is  a  step  function  with  a  peak  field 
strength  2nb  E* *  100V  (note  that  the  response  cannot  be  magnitude- 
normalized  as  in  the  examples  given  in  Chapters  5  and  6,  since  the 
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Figure  B.l  Input  current  of  a  circular  loop  antenna  driven  through  a  diode  by  a  double-exponential 
voltage  source.  Shown  is  also  the  input  current  with  the  diode  removed. 


Figure  B.2  Port  current  (solid)  and  voltage  (dashed)  responses  of  a  circular  loop  antenna  loaded 
by  a  diode  and  excited  by  a  step  function  plane  wave.  The  diode  is  polarized  initially  in  the 
"forward"  direction. 


Figure  B.3  Port  current  (solid)  and  voltage  (dashed)  responses  of  a  circular  loop  antenna 
loaded  by  a  diode  and  excited  by  a  step  function  plane  wave.  The  diode  is  polarized  initially 
in  the  "reverse"  direction. 


present  problem  is  nonlinear).  In  the  case  of  Figure  B.2  the  diode 
is  connected  so  that  it  is  polarized  in  the  "on"  direction  from  the 
beginning  and  has  little  irifluence  on  the  current  waveform  (of.  Figure 
5.8).  When  the  terminals  of  the  diode  are  interchanged  the  situation 
changes  drastically,  as  shown  in  Figure  B.3.  The  current  is  initially 
blocked  and  then  flows  in  the  reverse  direction.  Also  in  this  case, 
the  reverse  voltage  on  the  diode  can  reach  quite  high  levels  in  early 
times.  The  topology  and  element  values  of  the  equivalent  circuit 
used  are  defined  in  Tables  5.2  and  5.4.  Only  the  first  eight  dominant 
pole-pair  modules  were  employed.  Note  that  in  this  case  the  transfer 
networks  for  n*  1,3,5,..,  are  not  present,  as  indicated  in  Table  5.4. 

In  Figures  B.4  and  B.5  are  plotted  the  current  and  voltage  wave¬ 
forms  at  a  port  of  a  cylindrical  dipole  antenna  (Q=10.6)  center-loaded 
with  a  diode,  exicted  by  a  broadside  incident  plane  wave  with  a  double 
exponential  time  history  and  a  peak  field  strength  of  £E^ *  100V. 

In  Figure  B.4  the  diode  is  connected  so  that  it  is  initially  forward 
conducting.  In  Figure  B.5  the  terminals  of  the  diode  are  interchanged. 
As  expected  [Liu  and  Tesohe,  1976],  the  current  of  the  diode  for  tc/£>l 
oscillates  at  about  twice  the  rate  of  the  short-circuit  current  shown 
in  Figure  6.21.  This  phenomenon  occurs  because  the  diode,  when  reverse 
conducting,  acts  as  an  open  circuit.  The  antenna  thus  behaves  in 
effect  like  two  collinear  wires,  each  having  its  fundamental  natural 
resonance  at  a  frequency  roughly  twice  that  of  the  whole  structure. 

The  simplified  circuit  of  Figure  6.19  was  used  in  the  SCEPTRE  analyses, 
with  the  elements  defiend  in  Tables  6.2  and  6.12. 
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APPENDIX  C 


STATIC  CAPACITANCE  OF  THE  SPHERICAL  ANTENNA 


C.l  Introduction 

The  rapidity  of  the  convergence  of  the  eigenfunction  or  SEM 
series  for  the  driving-point  admittance  of  a  simply-connected  object 
can  be  significantly  increased  if  the  static  capacitance  of  the 
structure  is  known,  as  discussed  in  Chapter  7.  The  early  time  response 
of  an  SFM  equivalent  circuit  can  often  be  improved  by  adjusting  the 
static  capacitance  of  the  circuit  to  the  proper  value  as  demonstrated 
in  Chapters  6  and  7.  It  is  in  this  context  that  we  discuss  here  the 
computation  of  the  static  capacitance  of  the  spherical  antenna. 

In  the  first  part  of  this  appendix  approximate  closed  form  formulas 

are  derived  for  the  internal  (Cint)  and  external  (CeXt)  capacitances 

of  a  spherical  antenna  with  an  equatorial  slot.  The  derivation 

essentially  follows  that  given  by  Fvcmceschetti  [1976].  A  constant 

field  distribution  in  the  slot  region  is  assumed.  In  the  second  part 

a  numerical  technique  is  used  to  compute  the  total  capacitance  of 

the  spherical  antenna.  Results  are  presented  for  0^*90®  (equatorial 

slot)  and  9  *60°  with  the  slot  width  d/a  ”0.05  in  both  cases, 
g 


C.2  Derivation  of  Closed  Form  Expressions  for  the  Internal  and 
External  Capacitance 

In  the  equatorial  slot  case  the  expressions  (7.31)  and  (7.32) 
can  be  rewritten  as 
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V"/2,ao> 


P2n-l(no)/rr;?  1 

2n  2a 


P2n_i(n)dn 


-n. 


(C.  3) 


where  n“sina,  a  =  7r/2-0,  and  a  =d/(2a). 

o 

Following  Frcmoesahetti  [1976]  we  define 

f”(ot  )  -  i  V°°(a  )  (C.4) 

no  n  n  o 

00 

where  V  (a  )  is  the  asymptotic  expansion  of  V  (Tr/2,a  )  for  n-*-°°. 
no  no 

From  [ Frcmoesahetti s  1976] 

P2n-l(rl)  “  cos(mT)cOS[(2n-ls)a]  (C.5) 

it  follows  that 
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n  o  nTT 


n  * 


/cosa  cos[  (2n-Jj)  (a+aQ)  ]da 


-a 
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Now  the  series  (C.l)  and  (C.2)  can  be  recast  in  the  following  form 
„int  rra 
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(C.14) 


Using  (C.14)  and  (C.15)  it  can  be  shown  that 
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(C.17) 


The  substitution  of  (C.ll),  (C.16)  and  (C.17)  into  (C.8)  and  (C.9) 
yields 


and 


(C.18) 


,ext  „  2 

=  Z  c 
o 


(C.19) 


For  d/a » 0.05  these  expressions  give  Cints 70.66a  [pF]  and 

CeXt  -  104.74a  [pF],  so  that  the  total  capacitance  Ct0t * Cint  +  C6Xt s 


75.4a[pF] . 

It  is  of  interest  to  point  out  that  Sahelkunoff  [1952]  derived 
the  following  expression  for  the  capacitance  of  a  symmetrical  spherical 
antenna  driven  by  a  biconical  transmission  line: 


c  -zT  [■-©♦•■» 

O  *—  —I 


(C. 20) 


This  expression  differs  significantly  from  (C.19),  except  for  very 
narrow  slots. 


C.3  Computation  of  the  Total  Capacitance  by  the  MoM  Technique 

In  this  part  of  the  appendix  the  static  capr-itance  of  the  spherical 
antenna  is  computed  by  solving  numerically  the  Integral  equation  for 
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the  charge  density.  The  integral  equation  is  formulated  in  terms  of 
the  total  charge  density  on  the  spherical  shell,  i.e.a  it  cannot 
distinguish  between  the  charge  residing  on  the  internal  and  the 
external  side  of  the  shell.  As  a  result  the  capacitance  obtained  is  the 
total  capacitance,  which  is  a  sum  of  the  internal  and  external 
capacitances  discussed  above.  The  numerical  technique  applied  is  the 
Method  of  Moments  (MOM)  [Harrington,  1968]  with  subsectional  basis 
functions. 

The  charge  distribution  a(r)  on  the  spherical  antenna  (Figure 
7.1)  can  be  found  from  the  system  of  two  coupled  integral  equations 
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47TER 


f  q(r') 

4ttR 


dS' 


,  r  e  S1 


V2  ,  reS2 


(C.21) 


where  R =*  |r-r’|,  and  and  S2  denote  the  upper  and  lower  parts  of  the 

sphere,  respectively,  having  the  arbitrary  constant  potentials  and 

V2>  and  e  is  the  permittivity  of  the  medium.  In  order  to  find  the 

capacitance  C  of  the  structure,  the  system  (C.21)  must  be  solved  with 
12  12 

the  constraint  Q  ■  Q  * -Q  where  Q  and  Q  denote  the  total  charges 
on  Sj  and  S2»  respectively.  Then  we  have 


V  -  V 
1  2 


(C. 22) 


Exploiting  the  azimuthal  symmetry  of  the  problem,  eq.  (C.21) 
can  be  rewritten  in  the  form 
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and  K(8)  denotes  the  complete  elliptic  integral  of  the  first  kind 
[Abramo&itz  and  Stegun,  1965].  In  terms  of  the  new  function  p(0)  the 

total  charge  on  surface  is  given  by 

0  -a 

Q1  -  a  e  {  p(0)sin0d0  (C.27) 


0 


with  an  analogous  expression  for  S2> 

The  MoM  technique  is  applied  to  (C.23)  by  dividing  the  surface 
of  the  spherical  antenna  into  N  subsections  and  expanding  the  unknown 
function  p(0)  in  terms  of  the  pulse  functions  as  follows 


where 
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This  is  illustrated  by  Figure  C.l. 


The  substitution  of  (C.28)  into  (C.23)  and  matching  the  integral 


equations  at  0j  results  in  the  matrix  equation  for  the  expansion 
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where  the  index  i  takes  on  the  values  i*l,2,.,.,N,  and 

(N1  +  N2»N)  denote  the  number  of  zones  on  and  S^,  respectively. 

The  integration  is  done  numerically  by  a  Gaussian  quadrature. 

The  integrand  in  (C.31)  is  singular  for  i*j,  thus  the  integration 
for  the  self-term  cannot  be  performed  directly.  This  difficulty  is 
taken  care  of  by  subtracting  the  logarithmic  singularity  of  the  elliptic 
integral  and  integrating  the  singular  term  analytically  as  follows 
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In  order  to  gain  some  confidence  in  this  method,  the  capacitance 

1  2 

with  respect  to  infinity  of  a  solid  sphere  (aQ*0,  V^»V2“V,  Q«Q  +Q  ), 
C»Q/V,  was  computed  first.  With  N  ■  10  the  value  111.267  a  [pF]  was 
obtained,  which  is  in  an  excellent  agreement  with  the  exact  value  4irea. 
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In  order  to  find  the  capacitance  (or  "capacity"  {Adams,  1972]) 
of  the  slotted  sphere,  the  system  (B.23)  must  be  solved  together  with  the 
charge  neutrality  constraint.  However,  as  pointed  out  in  [Adams,  1972; 
Lee,  1978],  it  is  more  convenient  to  solve  the  unconstrained  system 
(C.23)  twice  for  even  and  odd  voltage  excitation  modes  (Figure  C.2) 
and  then  suitably  superpose  the  results.  It  can  be  shown  that  the 
capacity  is  then  given  by 


C  = 


qV  -  qV 

o  e  e  o 

2V(Qg+Q^) 


(C.33) 


where  the  superscripts  "e"  and  "o"  denote,  respectively,  the  even  and 
odd  excitation  modes. 

The  capacity  of  the  spherical  antenna  computed  by  this  scheme  is 

shown  in  Table  C.l  and  in  Figure  C.3  as  a  function  of  the  number  of 

zones  N.  The  extrapolated  values  of  the  normalized  capacitances  are 

165  [pF/m]  for  9  “90°  and  137  [pF/m]  for  0  *60°.  The  width  the 
8  8 

slot  was  d/a  *  0.05  in  both  cases.  The  value  165  a  [pF]  obtained  in 
the  equatorial  slot  case  agrees  favorably  with  the  value  175.4  a  [pF] 
obtained  from  the  closed  form  formulas  derived  in  the  first  part  of 
this  appendix,  considering  chat  approximations  are  made  in  both 
approaches. 
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Table  C.l  The  capacitance  per  unit  radius  of  a  spherical  antenna  with 
a  slot  of  width  d/a  "0.05  located  at  0g  *  90°  and  0g*6O#,  as  a  function 
of  the  number  of  zones  N  in  the  MoM  procedure. 


10  132.7  107.5 
20  145.9  119.4 
30  151.5  124.7 
40  154.6  127.5 
50  156.6  129.3 
60  158.0  130.6 


70 


159.0 


131.5 


d/a=0.05 


l/N 


Figure  C.3  Capacitance 
0  *  90°  and  at  9  ■  60°  as 
in  the  MoM  solution. 


of  the  spherical  antenna  with  a  slot  located  at 
a  function  of  l/N,  where  N  is  the  number  of  zones 
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